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SUMMARY 


Upon  discovery  of  Che  infrared  sensitivity  of  weakly  ionized  plasmas, 
research  was  initiated  at  establishing  the  quantum  origins  of  the  observed 
interactions.  Taking  into  account  the  obvious  inverse  bremsstrahlung 
process,  studies  were  made  to  evaluate  the  role  of  bound-bound  and  bound- 
free  transitions  under  the  influence  of  low  energy  photon  flux  corre- 
sponding to  the  intermediate  and  far-infrared  portions  of  the  spectrum. 

It  was  established  and  documented  that  in  a broad  range  of  electron 
temperatures  and  densities  of  weakly  ionized  plasmas,  the  bound-free  tran- 
sitions induced  by  Infrared  photons  affect  the  recombination  process  to  a 
significant  and,  frequently,  to  a dominant  degree.  Theoretical  results  are 
obtained  on  the  basis  of  a hydro genic  model.  Experiments  demonstrate  the 
generality  of  the  infrared  sensitivity  of  plasmas. 

A separate  phase  of  research  in  optical  quantum  electronics  dealt 
with  optical  interaction  with  surface  and  near  surface  electrons  in  photo- 
emissive  materials.  In  particular,  the  dynamics  on  picosecond  time  scales 
of  electrons  released  from  the  surface  and  near-surface  of  select  materials 
is  of  both  fundamental  and  practical  interest.  Negative  electron  affinity 
materials,  known  for  their  high  quantum  efficiency,  are  predicted  to  exhibit 
maximum  time  dispersion  of  a step  excitation.  Preliminary  studies  and  the 
development  of  tools  needed  for  an  experimental  determination  of  parameters 
on  picosecond  time  scales,  were  the  principal  goals  of  this  effort.  As  a 
byproduct,  a patent  was  filed  and  granted  on  the  more  successful  configu- 
rations of  devices  designed  to  carry  out  the  basic  surface  emission  studies. 
The  program  is  currently  well  under  way  and  results  continue  to  emerge  in 
an  area  that  was  heretofore  not  amenable  to  experimental  exploration  due  to 
the  unavailability  of  picosecond-resolving  devices. 
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EVIDENCE  OF  LOW-BINDING-ENERGY  ELECTRONS 
IN  A LOW -TEMPERATURE  HELIUM  AFTERGLOW* 

J.  P.  Kaplafka,  H.  Merkelo,  and  L.  Goldstein 
Gaaeous  Electronics  Laboratory,  University  of  Illinois,  Urbana,  Illinois 
(Received  12  July  1968) 

Moderate,  pulsed  microwave  heating  of  the  free  electrons  In  a weakly  ionized  helium 
afterglow  plasma  at  9.5  Torr  at  a gas  temperature  of  77"K  produces  an  unusual  modula- 
tion of  the  visible  afterglow  light.  Observations  of  a second  (probing)  microwave  signal 
indicate  an  accompanying  significant  increase  in  free-electron  density  under  the  assump- 
tion that  there  is  no  significant  anomalous  phase  shift  or  absorption  of  the  probing  signal. 


Observations  have  been  made  of  the  modulation 
of  both  visible  afterglow  radiation  and  microwave 
transmission  resulting  from  perturbations  of  the 
free-electron  gas  temperature  imposed  during 
various  intervals  in  the  decay  of  a weakly  ionized 
helium  gas  at  a pressure  of  9.5  Torr  maintained 
at  a gas  temperature  near  77°K.  The  perturba- 
tion of  the  free-electron  temperature  Te  is  pro- 
duced by  the  introduction  of  a second  microwave 
signal  of  a different  frequency  and  necessarily 
higher  power  than  the  probing  microwave  signal. 
This  “heating”  microwave  signal  is  «200  Msec  in 
duration.  Under  the  tentative  assumption  that 
there  is  no  significant  anomalous  phase  shift  or 
absorption  of  the  transmitted  probing  microwave 
signal,  interpretation  of  the  modulation  on  phase 
shift  and  absorption  implies  a significant,  unex- 
pected increase  in  free-electron  density  ne  (on 
the  order  of  unperturbed  ne)  in  the  course  of 
moderate  microwave  heating  and  a subsequent 
fast  (“50-100  nsec)  return  to  unperturbed  ne  up- 
on removal  of  heating. 

The  discharge  tube  is  Pyrex,  rectangular,  ap- 
proximately 24  cm  in  length,  and  makes  a close 
fit  within  an  .Y-band  microwave  waveguide.  A 2- 
Msec,  high-voltage  (~7-kV)  pulse  applied  between 
electrodes  located  outside  the  waveguide  produc- 
es breakdown  of  the  cataphoretically  pure  helium. 
At  200  Msec  after  breakdown,  ne  is  found  to  be 
“10ll/cc,  falling  to  lO’/cc  in  approximately  4 
msec.  Spectrally  resolved  observations  of  visi- 
ble light  intensities  in  the  early  afterglow  (within 
800  Msec)  were  made  with  a Jarrell-Ash,  Model 
82-000,  0.5-m  monochromator  and  an  RCA  photo- 
tube (S-ll  photocathode).  The  decay  of  atomic 
He  lines  is  markedly  different  than  that  of  di- 
atomic molecular  He  bands.  The  lines  are  very 
strong  in  the  breakdown  but  decay  rapidly  in  in- 
tensity. By  approximately  300  Msec  the  sum  of 
atomic  line  contributions  to  the  signal  from  the 
photomultiplier  can  be  neglected  in  comparison 
with  the  sum  of  molecular  band  contributions. 


The  most  prominent  six  molecular  bands  (5735  A 
[3*Au-2*ng],  4725  A [4  ‘Su+-2  *%],  4650  A 
(3sng-2*Zu+],  4547  A [4,ru+-2sng],  -4470  A, 
and  4439  A [4snu-2sng])1’*  can  be  followed  with 
15 -A  resolution  to  beyond  800  Msec.  These  six 
bands  appear  to  decay  at  the  same  rate.  Beyond 
800  Msec  the  observations  of  light  intensity  are 
not  spectrally  resolved;  i.e.,  they  are  observa- 
tions of  the  signal  from  the  photomultiplier  with 
total  visible  radiation  incident  on  the  photocath- 
ode. 

When  a short  pulse  of  microwave  radiation  of 
Insufficient  power  to  cause  significant  ionization 
of  ground-state  or  metastable -state  He  atoms  is 
incident  on  an  He  afterglow  plasma  of  gas  tem- 
perature near  300  °K  at  times  on  the  order  of  a 
millisecond  after  the  discharge,  the  visible  af- 
terglow light  intensity  is  partially  quenched  dur- 
ing application  of  the  pulse.3’4  Generally  such  a 
pulse  of  constant  power  will  maintain  a constant 
free-electron  temperature  Te  (after  the  time 
necessary  to  establish  this  Te)  during  its  appli- 
cation even  in  the  event  that  ne  decays  appreciab- 
ly during  its  application.  If  the  pulse  length  is 
short  (i.e.,  no  significant  change  of  ne  by  recom- 
bination or  diffusion  losses  occurs  during  appli- 
cation), then  the  afterglow  light  is  observed  to 
maintain  a constant  level  until  the  removal  of  the 
microwaves.  At  this  time  the  light  intensity  re- 
turns to  approximately  the  same  intensity  as  be- 
fore introduction  of  the  microwaves,  in  the  time 
required  for  Te  to  decay  by  means  of  elastic  col- 
lisions between  the  free  electrons  and  helium  at- 
oms.4 The  partial  quenching  of  visible  afterglow 
light  intensity  has  been  interpreted  as  due  to  a 
reduction  of  the  recombination  coefficient  by  in- 
crease in  Te  which  results  in  subsequent  reduc- 
tion of  the  light  from  recombination. 

The  upper  traces  of  Fig.  1 illustrate  the  very 
different  effect  produced  on  the  visible  afterglow 
light  of  a 9. 5 -Torr,  helium  afterglow  plasma  at 
a gas  temperature  of  77 °K.  The  electron  temper- 


870 


I 


•i 


2 


p 


BJMMWUWHIWW  JAWiWW^tlMlWW^ aiWSIIl  ""IM- 


B trnmmm 


I 


Volumi  21,  Numbbb  U PHYSICAL  REVIEW  LETTERS  30  SirriMBcm  1968 


. FIG.  1.  Effects  of  moderate  microwave  heating  upon 
the  visible  afterglow  light  (upper  traces)  and  upon  a 
transmitted  probing  microwave  signal  Gower  traces) 
for  a He  afterglow  plasma  at  a gas  pressure  of  9.5 
Torr  and  gas  temperature  of  77°  K. 


ature  increases  rapidly  (~3  psec)  at  the  time  of 
application  of  the  microwave  pulse  tu  accompa- 
nied by  an  expected  rapid  partial  quenching  of  the 
visible  afterglow  light.  However,  visible  after- 
glow light  subsequently  increases  and  reaches  a 
nearly  constant  amplitude  within  approximately 
150  Msec.  This  amplitude  is  nearly  coincident 
with  the  unperturbed  amplitude  when  the  heating 
is  relatively  small  (this  is  the  case  in  Fig.  1)  and 
progressively  less  than  the  unperturbed  ampli- 
tude as  heating  power  is  increased.  At  the  time 
of  removal  of  the  heating  microwaves,  visible  af- 
terglow light  intensity  increases,  v/ithin  approxi- 
mately 2 psec,  to  an  amplitude  that  can  be  sev- 
eral times  the  unperturbed  and  then  decays  to  the 
unperturbed  amplitude  in  approximately  103-200 
Msec.  Such  a transient  increase  in  afterglow 
light  upon  removal  of  heating  microwaves  was 
observed  in  He  afterglows  at  gas  temperatures 
of  approximately  77  and  4.2°K  by  Goldan,  Berlan- 
de,  and  Goldstein  and  termed  the  “afterpulse.  ”* 
The  lower  traces  of  Fig.  1 illustrate  amplitude 
modulation  produced  on  the  low -power  probing 
microwave  signal  (applied  over  a longer  interval 
of  time  than  the  microwave  heating  pulse)  due  to 
application  of  the  heating  pulse.  The  rapid  in- 
crease in  the  absorption  at  f,  is  in  response  to  a 
rapid  increase  in  Te,  inferred  from  the  increase 
in  the  elastic  electron -neutral  collision  frequen- 
cy for  momentum  transfer  as  determined  from 


FIG.  2.  Typical  plot  of  electron  density  at  150  Msec 
of  microwave  heating  versus  electron  temperature  dur- 
ing heating.  Each  has  been  normalized  to  the  respec- 
tive unperturbed  value  for  this  particular  afterglow 
time.  (T^0~  170°K  for  this  plot,  gas  temperature  of 
77°K  and  pressure  of  9.5  Torr.) 


absorption  and  phase -shift  measurements  on  the 
transmitted  probing  microwave  signal.  The  sub- 
sequent more  gradual  increase  in  absorption  dur- 
ing application  of  the  heating  microwaves  is 
found  to  result  from  an  increase  in  ne-Te  re- 
maining approximately  constant.  In  the  case  of 
Fig.  1 the  increase  in  Te  by  heating  is  about 
100%  and  the  modulations  of  the  probing  micro- 
wave  signal  and  the  visible  afterglow  light  are 
typical  of  this  relative  increase  in  Te  throughout 
the  region  where  most  observations  were  made, 
i.e.,  between  1 and  3 msec. 

Figure  2 is  a typical  plot  of  ne,  as  inferred 
from  the  absorption  and  phase-shift  measure- 
ments, at  150  Msec  of  microwave  heating  versus 
Te  during  heating,  both  normalized  to  their  re- 
spective unperturbed  values.  It  should  be  noted 
that  Teg  at  1 msec  of  the  afterglow  (i.e.,  unper- 
turbed Te  at  1 msec)  is  still  falling  and  well 
above  77°K  (~190  and  ~150°K  at  2 msec).®  Ambi- 
polar  diffusion  enhanced  by  increasing  Te  must 
be  significant  in  determining  the  shape  of  the  lat- 
ter portions  of  the  curve,  but  nonetheless  ne  can 
be  increased  by  a factor  of  almost  3 during  mi- 
crowave heating,  and  merely  doubling  Te  gives 
rise  to  an  almost  50%  increase  in  ne. 

When  the  heating  microwaves  are  suddenly  re- 
moved at  tit  Te  falls  rapidly  to  a value  near  un- 
perturbed Te,  accompanied  by  the  almost  discon- 
tinuous fall  in  absorption  seen  in  Fig.  1.  The 
more  gradual  change  in  absorption  after  this  fall 
results  mainly  from  a decrease  in  ne  to  unper- 
turbed ne. 
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The  modulation  of  the  unresolved  visible  light 
intensity  is  assumed  to  involve  no  spectral  redis- 
tribution of  the  observed  radiation  because  the 
six  major  molecular  bands  listed  earlier  show 
apparently  identical  modulations,  relative  to 
their  own  unperturbed  intensities,  for  early  af- 
terglow times  where  they  can  be  spectrally  re- 
solved. 

The  preceding  remarks  were  based  on  the  as- 
sumption that  the  free -electron  density  and  elec- 
tron temperature  measurements  are  reliable. 
However,  the  reliablility  of  determining  free- 
electron  density  and  electron  temperature  from 
phase-shift  and  absorption  measurements  on  the 
probing  microwave  signal  rests  with  the  assump- 
tion that  the  phase  shifts  and  absorptions  are  re- 
lated to  the  free  electrons  undergoing  elastic 
electron-atom  collisions.  The  following  are 
some  additional  possibilities: 

(1)  Inelastic  electron  collisions.  However,  at 
a given  time  in  the  afterglow  a significant  effect 
on  phase  shift  or  absorption  due  to  free  electrons 
undergoing  inelastic  collisions  requires  a popu- 
lation of  “target”  states  several  order  of  magni- 
tude larger  than  the  free-electron  density  or  in- 
elastic collision  cross  sections  several  orders 
of  magnitude  larger  than  the  elastic  electron- 
ground -state -atom  collision  cross  section  for 
momentum  transfer. 

(2)  Resonant  absorption  of  probing  microwaves 
(i.e.,  absorption  by  means  of  photon -induced 
transitions  between  bound  states).  The  probing 
microwave  photon  energy  is  small  (~3X10_5  eV). 
Therefore  it  is  necessary  to  consider  states 
which  are  situated  close  in  energy.  It  is  difficult 
to  explain  the  modulation  of  visible  molecular  af- 
terglow light  if  states  responsible  for  a resonant 
microwave  absorption  are  solely  states  of  the  at- 
om. Moreover,  it  appears  that  if  the  probing  mi- 
crowave signal  is  significantly  absorbed  as  the 
result  of  probing  microwave  photons  inducing 
transitions  between  molecular  rotational  levels, 
the  density  of  molecules  in  such  rotational  levels 
would  be  nearly  as  large  as  the  density  of  helium 
atoms  initially  present  in  the  gas  (~10“/cc). 

(3)  Particle  redistribution  or  oscillation.  The 
possibility  that  the  observed  perturbations  are 
due  to  particle  redistribution  was  examined  in 
several  ways:  (a)  by  spatially  resolved  observa- 
tions of  afterglow  light  across  the  width  and 
along  the  length  of  the  discharge  tube,  (b)  by 
heating  in  the  TE^,  microwave  mode  rather  than 
TEI0  and  observing  visible  afterglow  light,  and 
(c)  by  observing  visible  afterglow  light  emitted 


by  a discharge  tube  situated  outside  the  micro- 
wave  waveguide  and  irradiated  by  means  of  a mi- 
crowave horn.  There  was  no  evidence  in  these 
cases  to  support  the  possibility  of  spatial  redis- 
tribution. 

Tentative  acceptance  of  the  reliability  of  the 
measurements  of  ne  and  Te  leads  to  certain  im- 
plications: The  rise  of  afterglow  light  (after  the 
initial  partial  quenching)  in  the  course  of  heating 
could  be  explained  on  the  basis  of  increasing  rate 
of  recombination,  caused  by  an  increasing  ng. 

The  afterpulse  of  light  could  be  understood  as  a 
response  to  the  sudden  increase  in  recombina- 
tion coefficient,  the  result  of  a rapid  decrease  in 
Tg.  It  should  be  noted  that  since  it  is  difficult  to 
determine  at  what  time  ne  actually  returns  to  un- 
perturbed ne,  it  should  not  be  inferred  from  Fig. 

1 that  the  light  in  the  tail  of  the  afterpulse  contin- 
ues above  its  unperturbed  level  beyond  this  time. 

The  functional  dependence  of  ne  on  time  t after 
removal  of  the  microwave  heating  before  ne  re- 
turns to  normal  cannot  be  explained  by  invoking 
a simple  recombination  law.  Assuming  that  dne/ 
dtccng*7,  slopes  of  plots  of  In (dne/dt)  vs  In ne  give 
inconsistent  and,  at  times,  inordinately  high  val- 
ues of  a (~3  «r  <~5),  Te  being  essentially  con- 
stant. Moreover  dne/dt  within  this  period  is 
much  greater  than  can  be  expected  on  the  basis 
of  the  decay  rates  of  ne  observed  in  the  unper- 
turbed afterglow.  The  additional  observations 
that  ne  returns  to  unperturbed  ne  and  that  dne/dt 
at  the  time  of  application  of  the  heating  micro- 
waves  increases  monotonically  with  Te  (i.e., 
monotonically  with  microwave  heating  power) 
suggest  the  following  interpretation:  There  is  a 
population  [ AT]  of  weakly  bound  electrons  in  a 
near-collisional  equilibrium  with  the  free-elec- 
tron gas  at  a given  afterglow  time.  Microwave 
heating  of  the  free-electron  gas  has  the  effect  of 
disturbing  this  equilibrium  by  increasing  the  de- 
struction rate  of  [x]  and  perhaps  reducing  the 
rate  of  formation.  The  electron  density  increas- 
es during  heating  until  a new  equilibrium  popula- 
tion of  state  or  states  X is  established  at  the  en- 
hanced Te  maintained  during  heating.  After  the 
microwave  heating  pulse  is  removed,  Te  falls  to 
a value  near  unperturbed  Te  and  [x]  is  subse- 
quently re-established  at  the  value  previous  to 
heating.  ' 

At  present  additional  efforts  are  being  consid- 
ered to  confirm  further  the  existence  of  a signifi- 
cant population  of  weakly  bound  electrons. 
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QUENCHING  OF  HELIUM  AFTERGLOW  BY  INFRARED  PHOTONS 

J.  P.  Kaplafka , H.  Merkelo,  and  L.  Goldstein 
Gaseous  Electronics  Laboratory 
University  of  Illinois 
Urbana,  Illinois  61801 

(Received  3 June  1969;  in  final  form  16  July  1969) 


The  visible  afterglow  radiation  from  low-pressure  weakly  ionized  helium  exhibits  a partial  quench- 
ing due  to  absorption  of  infrared  photons  (0.11-0.14  eV).  The  helium  sample  is  placed  in  the  cavity 
of  a Q- switched  0O2  laser  and  studied  at  gas  temperatures  In  the  range  from  77  to  300°K.  The  par- 
tial quenching  appears  to  originate  from  photon-induced  change  of  population  of  a state  or  states  of 
a neutral  or  ionic  helium  molecule. 


I 


A pulsed  helium  discharge  tube  was  placed  in  the 
cavity  of  a Q-switched  C02  laser.  A partial  quench- 
ing of  the  visible  afterglow  radiation  from  the 
weakly  ionized  helium— as  large  as  50%  at  a gas 
temperature  near  77“K— is  observed  during  the 
brief  (“- 10  usee)  time  interval  of  lasering.  Gas 
pressure  was  varied  in  a range  from  a few  tenths 
of  a Torr  to  approximately  15  Torr,  but  the  quench- 
ing Is  not  observable  below  ~0.4  Torr  at  a gas 
temperature  of  77°K  or  below  ~2  Torr  at  room 
temperature.  (Below  these  pressures  the  visible 
afterglow  radiation  consists  primarily  of  atomic 
lines  rather  than  molecular  bands.)  It  is  present 
throughout  all  but  the  very  early  afterglow  period 
and  appears  to  originate  from  photon- induced 
change  of  population  of  a state  or  states  of  a neu- 
tral or  ionic  helium  molecule. 

A major  portion  of  the  helium  discharge  tube 
(~50  cm  in  length  and  10  mm  i.d.)  is  contained 
within  a styrofoam  box  to  allow  cooling  of  the  he- 
lium gas  in  the  region  where  observations  are 
made  of  the  visible  afterglow  radiation  emitted  fol- 
lowing application  of  a 2-  psec  several-thousand- 
volt  discharge  pulse.  The  helium  is  cataphoretically 
cleaned.  The  visible  afterglow  radiation  is  trans- 
mitted to  the  input  of  an  RCA  photomultiplier  having 
an  S-ll  photocathode.  <?  switching  of  the  CO,  laser 
Is  accomplished  by  means  of  rotating  one  of  the 
two  spherical,  cavity  mirrors  at  a rate  of  60  rps. 
Lasering  occurs  every  16.6  msec  in  the  form  of 
aeveral  short  pulses  spanning  an  interval  of  approx- 
imately 10  psec.  The  time  interval  during  which 
lasering  takes  place  can  be  positioned  to  any  time 
during  the  helium  afterglow.  The  onset  of  quench- 
ing is  essentially  coincident  in  time  with  the  onset 
of  lasering.  The  measured  rise  times  of  the 
quenched  light  signals  are  found  to  be  shorter  than 
0,5  psec  in  all  cases. 

With  the  helium  tube  in  a liquid-nitrogen  bath 
(77*10  the  percentage  of  quenching  at  most  pres- 
sures can  be  conveniently  obtained  up  to  3 msec 
after  application  of  the  2-psec  helium-discharge 
voltage  pulse.  Within  the  precision  of  measurement 
of  the  change  in  visible  radiation  produced  during 
lasering,  the  percentage  of  quenching  varies  little, 
after  reaching  a maximum  in  the  very  early  after- 
glow, in  comparison  to  the  changes  of  parameters 


of  the  electron  gas.  Based  on  microwave  sensing 
measurements  on  systems  under  similar  conditions 
of  gas  temperature  (i.e.,  near  77°K),  gas  pressure, 
breakdown,  and  discharge  tube  geometry,  it  is  ex- 
pected that  the  free-electron  density  (~10lo/cc) 
changes  by  at  least  a factor  of  10  and  the  free- 
electron  temperature  by  at  least  a factor  of  3. 

In  order  to  verify  directly  that  the  partial 
quenching  is  not  due  to  an  increase  of  temperature 
of  the  electron  gas,  the  helium  tube  was  inserted 
through  a section  of  square  cross-section  micro- 
wave  guide.  Sensing  microwaves  at  a frequency  of 
8.6  GHz  were  transmitted  through  this  guide  section 
at  a low  power  level  (so  as  not  to  disturb  the  elec- 
tron-gas temperature)  and  observations  were  made 
of  the  transmitted  signal  during  the  interval  of  la- 
sering. (An  increase  in  temperature  of  the  electron 
gas  will  cause  an  increase  in  the  absorption  of  the 
microwave  signal  by  the  He  plasma.)  There  was  no 
clearly  discernible  change  in  absorption  of  the 
sensing  microwaves  during  the  interval  of  lasering. 

In  order  to  verify  the  sensitivity  of  this  method  a 
proven  means  of  changing  electron  temperature 
was  Introduced  for  the  purpose  of  comparison.  In 
order  to  selectively  heat  only  the  free-electron 
gas,  a second  microwave  source  operating  at  a fre- 
quency of  10.2  GHz  was  pulsed-on  for  a short  time. 

The  signal  was  transmitted  through  the  guide  sec- 
tion containing  the  helium  tube,  and  its  power  was 
adjusted  so  that  It  produced  the  same  percentage  of 
quenching  of  the  visible  afterglow  radiation,  at  the 
afterglow  time  of  comparison,  as  was  observed 
during  lasering.  The  change  in  transmission  (in- 
creased absorption)  of  the  low  power  (sensing)  mi- 
crowave signal  was  clearly  visible,  as  illustrated 
in  Fig.  1.  It  is  interesting  to  note  that  the  helium 
afterglow  has  this  capability  of  detecting  and  dis- 
tinguishing between  the  incident  infrared  and  mi- 
crowave radiations.  A partial  quenching  by  the 
infrared  irradiation  is  seen  even  during  simultane- 
ous application  of  a pulsed  microwave  signal  of 
sufficient  power  to  very  substantially  increase  the 
electron  temperature. 

Since  the  partial  quenching  by  the  Infrared  pho- 
tons with  helium  is  not  due  to  a heating  of  the  elec- 
tron gas,  it  appears  that  the  lasering  produces  a 
partial  destruction  of  the  population  of  an  excited 
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Fig.  1.  Trace  1 is  visible  helium  afterglow 
radiation  between  ~840  and  940  pse c after 
discharge  (i.e.,  10  psec/cm)  increasing 
downward.  The  radiation  is  shown  partially 
quenched  by  infrared  lasering  and  just  adja- 
cent to  this  it  is  partially  quenched  by  a 
square  pulse  of  microwave  power.  Trace  2 
is  the  transmitted  low  power  (sensing)  mi- 
crowave signal  clearly  perturbed  (increased 
absorption)  during  microwave  heating  and 
virtually  undisturbed  during  lasering.  (Note 
that  the  lower  trace  is  displaced  slightly  to 
the  right  on  the  screen.)  The  gas  pressure 
was  8 Torr  and  gas  temperature  ~ 77‘K. 


specie(s)  existing  somewhere  in  the  recombination 
cascade  (even  a state  of  the  ion  itself)  that  leads  to 
the  visible  afterglow  radiation. 

It  was  previously  observed  by  the  authors,  at 
gas  temperature  near  77'K  and  pressures  of  sev- 
eral Torr,  that  with  similar  breakdown  and  in  the 
same  discharge  tube  geometry  the  visible  afterglow 
radiation  after  very  early  afterglow  times  is  com- 
pletely dominated  by  band  radiation  from  the  He2 
molecule.1  The  He2  ion  has  been  identified  in  room- 
temperature  helium  afterglows  by  means  of  mass 
spectrometry  and  recombination  via  this  ion  is  as- 
sumed responsible  for  the  observed  He2  band  radi- 
ation in  room-temperature  afterglows.  Recently, 
however,  heavier  helium  ions  have  been  identified 
at  lower  gas  temperatures.  The  existence  of  an  HeS 
ion  over  a range  of  gas  temperatures  from  76  to 
200CK  and  pressures  from  1 to  30  Torr  was  in- 
ferred by  means  of  drift-tube  ion- mobility  experi- 
ments.2,3 In  the  course  of  these  studies,  Patterson3 
made  the  assessment  that  at  77  K the  rate  of  for- 
mation of  Ile2  from  He2  is  faster  than  the  three- 
body  formation  of  He2  from  He*.  In  this  regard  a"d 


as  a direct  proof  of  its  existence,  Ferguson  ct  al* 
found  the  Hej  ion,  by  means  of  mass  spectrometry, 
in  a flowing  helium  afterglow  at  a gas  temperature 
of  82"K.  In  fact  Hej  was  the  dominant  ion  they  ob- 
served. We  have  since  been  informed  that  deVries 
and  Oskam  have  also  identified  the  Hej  ion  by  mass 
spectrometry,  in  helium  afterglows  at  gas  temper- 
atures near  77‘K,  and  in  addition  discovered  an  Hej 
ion  in  these  same  afterglows.5  Recombination  of 
these  necessarily  weakly  bound  heavier  ions  is 
likely  to  be  a source  of  excited  states  of  He2  and 
thereby  responsible  for  the  visible  band  radiation 
that  we  observe  at  the  lower  gas  temperatures. 

When  the  gas  temperature  is  maintained  near 
77°K  the  percentage  of  quenching  is  largest  near 
0.7  Torr  and  decreases  with  increasing  pressure. 

Preliminary  observations  at  a gas  temperature 
near  77°K  and  a pressure  of  several  Torr  using 
spectral  resolution  with  the  present  system,  at 
early  afterglow  times  where  some  intensities  are 
adequate  for  resolution  and  where  the  rapidly  de- 
caying 5876  A (3 3D-  2 3P)  atomic  line  can  still  be 
observed  « ~300  psec),  indicate  that  the  partial 
quenching  during  lasering  is  present  and  the  per- 
centage of  quenching  constant  on  the  strongest  mo- 
lecular bands  (e.g.,  4439  A [43n„  - 23n<f  ],  4725  A 
[412j2-  2'nJ)  and  apparently  absent  on  the  strong- 
est atomic  line  (5876  A). 

The  visible  afterglow  quenching  by  the  infrared 
photons  of  the  C02  laser  appears  related  to  recom- 
bination of  the  molecular  ions  present  and  is  due 
to  photon -induced  changes  in  population  of  a state 
or  states  of  a neutral  molecule  formed  in  recom- 
bination or  a change  in  population  of  a state  of  the 
ion  that  contributes  significantly  to  recombination. 
In  either  case  the  response  to  the  infrared  irradi- 
tion  can  be  expected  to  be  very  rapid. 

We  wish  to  thank  Professor  H.  J.  Oskam  for  his 
preprint  regarding  mass  analysis  of  the  low  gas 
temperature  helium  afterglow.  We  also  wish  to 
acknowledge  the  able  technical  assistance  of  W. 
Johnson,  K.  Kuehl,  E.  Boose,  and  J.  James. 
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Th*  purpose  of  Chis  Celk  is  to  report  briefly  on  a preliminary  study  of  the  interaction  of  infrared 
radiations  with  gaseous  plasmas  and  on  the  possible  applications  of  the  resulting  effects.  Particular 
emphasis  is  laid  on  the  use  of  gaseous  plasmas  for  detection  of  pulsed  infrared  radiation,  in  the  9-11 
micron  wavelength  range  of  CO^  lasers,  in  the  presence  and/or  absence  of  pulsed  microwaves. 

1.  Generalities  on  the  plasma  media  considered.  Interaction  with  high  frequency  electromagnetic  waves. 


For  simplicity  we  shall  first  consider  weakly  ionized,  low  pressure  (<10  Torr)  noble  gases  and/or 
some  of  their  mixtures. 

The  plasmas  are  produced  by  electric  fields,  either  pulsed  or  continuous,  between  electrodes.  The 
configurations  and  locations  of  these  electrodes  are  compatible  with  the  propagation  of  Infrared  and 
microwave  radiations  through  the  plasma.  The  gases  investigated  are  housed  in  containers  that  are  trans- 
parent to  the  visible  radiations  they  emit  when  excited. 

He  shall  restrict  generality  for  the  purpose  of  this  summary  and  consider  only  pulsed  Infrared  in 
the  9 to  11  micron  wavelength  range  (from  CO^  lasers)  and  microwaves  in  the  3-10  cm  wavelength  range. 

1-1,  A weakly  ionized  gas  can  be  generally  considered  as  a mixture  of  the  gases  of  neutral  atoms  and  the 
gases  of  the  charged  particles  l.e.,  the  free  electrons,  e,  and  the  positive  ions,  1.  These  are  in 
continuous  Interaction  with  each  other.  An  equation  of  state  for  such  a mixture  can  be  written,  to  a 


good  approximation,  in  a form  expressing  the  law  that  the  total  gas  pressure,  p^,  in  the  volume  containing 
this  gas  is  the  sum  of  the  partial  pressures  of  each  constituent: 

Pt  - P.  + Pi  ♦ Pn 

Or,  according  to  the  kinetic  theory  of  gases, 

Pt  * % ^ I,  + k Ij  + H k Ij 

where  (n  , T ),  (n. , T.)  and  (N,  T ) represent  the  density  and  temperatures  of  the  electrons,  ions  and 

• « 1 1 g 

neutral  gas  atoms  respectively.  A certain  fraction  of  the  neutral  atoms  is  in  various  excited,  radiating 
or  nonradiating,  states*. 


Strictly  speaking  we  should  also  include  a photon  gas  component.  However,  in  view  of  circumstances 
here  considered,  we  may  neglect  it  with  no  great  loss  of  understanding;  the  plasmas  being  transparent  for 
■oat  of  their  own  radiations. 
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EFFECTS  OF  INFRARED  RADIATION  ON  CASEOUS  PLASMAS 
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Any  phenomenon  observed  result*  from  the  Interaction  between  the  constituents.  Since  wa  excite  these 
gases  by  electrical  means,  tha  electrons  and  Iona  gain  and  accumulate  kinetic  energy  from  the  field. 

There  are  a great  number  of  collisions  and  these  tend  to  randomize  their  motions.  Thus  their  temperatures 
are  higher  than  that  of  the  neutral  gas.  The  Ions,  having  the  mass  of  the  atoms,  are  in  excellent  thermal 
contact  with  the  neutral  gas  atoms  and  reach  a temperature  only  very  slightly  higher  than  that  of  the 
neutral  component.  The  electron  gas,  however,  because  of  poor  thermal  contact  with  the  heavy  particles, 
reaches  temperatures  which  may  be  several  orders  of  magnitude  higher  than  gas  temperature.  Thus  the 
electrons  are  responsible  for  both  ionization  and  excitation  of  the  neutral  gas  atoms,  which  take  place 
either  directly  from  ground  state  or  from  some  already  excited  state. 

As  long  as  energy  Is  supplied  to  the  electrons  from  an  external  field,  this  situation  prevails  and 
T^  » T^  <v  Tg.  The  electron  temperature  rise,  or  even  better,  the  mean  energy  of  the  electron  gas  is 
limited  by  geometrical  conditions.  For  these  low  pressure  cases  the  geometrical  conditions  determine 
the  loss  rates  for  charged  particles.  Steady  state  conditions  correspond  to  balance  between  charged 
particle  production  and  losses.  As  long  as  ionization  of  the  gas  la  taking  place  from  the  ground  stats  of 
the  atoms,  the  plasma  Is  said  to  be  active. 

1-2.  The  Inactive  plasma:  the  afterglow. 

When  the  source  of  excitation  and  maintenance  of  the  plasma  Is  removed,  the  electron  temperature 
decreases  toward  the  gas  temperature  chiefly  by  thermal  contact  with  the  neutral  gas.  Once  the  electron 
temperature  decays  sufficiently  that  there  remain  very  few  electrons  to  produce  further  Ionization — even 
from  excited  (metastable)  atoms — the  plasma  becomes  inactive.  Such  a plasma  decaya  in  a length  of  time 
determined  by  the  Initial  conditions  and  the  loss  rates  of  the  remaining  charges. 

The  charged  particles  are  lost  by  volume  recombination  and  surface  recombination.  The  excited  atoms 
decay  rapidly  (<  10  6 sec)  by  radiation  and  those  In  nonradiating  states  by  other  processes.  Those 
charges  which  do  not  undergo  volume  recombination  reach  the  bounding  surfaces  by  diffusion  (a  process 
highly  dependent  on  geometrical  conditions). 

In  what  follows  we  shall  be  dealing  with  plasmas  In  low  pressure,  monatomic  gases  at  ordinary 
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temperatures.  The  degree  of  Ionization  l.e.  — , can  be  as  small  as  10  to  10  . The  charge  density 
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range  la  ilO  to  10  /cc  with  gas  densities  of  10  to  IO'1  atoms/cc. 

On  the  other  hand,  from  the  point  of  view  of  electron  temperature,  we  shall  discuss  two  kinds  of 
such  plasmas: 

T 

a)  where  — » 1 as  In  continuously  excited  active  plasmas 
T. 

and  b)  where  y~  x 1 •*  1"  the  case  of  decaying  Inactive  afterglowlng  plasmas. 

Investigations  of  volume  recombination  processes  In  Inactive  helium  plasmas,  In  a range  of  gas  temp- 
erature from  77*K  to  300*K,  led  us  to  Irradiate  such  plasmas  with  low  energy  photons  (hv  < 0.15  eV). 


These  were  conveniently  provided  by  COj  lasers. 

1- 3 . Effect  of  CO.  loser  radiation  on  inactive  plasmas  produced  in  low  pressure  noble  gases.  Tha 

underlying  thought  for  using  the  'ilOum  infrared  quanta  was  essentially  the  following: 
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When  an  electron  recombines  with  an  atonic,  positive  Ion,  the  noutrallzed  Ion  nay  find  Itself  on  any 
of  the  possible  excited  states  of  the  atom.  Cascading  toward  the  ground  state,  auch  an  atoa  ealts  light 
quanta — some  of  which  are  In  the  visible  part  of  the  spectrua.  The  total  light  thus  ealtted  by  all  atoaa 
constitutes  the  afterglow,  which  con  be  spectrally  resolved. 

In  the  aost  probable  process  for  such  a recombination  (which  Involves  two  electrons  with  one  positive 
Ion)  there  Is  high  probability  for  the  captured  electron  to  find  Itself  Initially  on  a very  high  lying 
excited  state  of  the  atoa. 

Such  a state  Is  characterized  by  a high  principal  quantua  number  (n  > 8, 9, 10, etc. ) . The  energy 
level  of  such  an  excited  state,  E^,  is  quite  close  to  the  Ionization  energy  of  the  atoa,  E^.  If  thare 
are  highly  excited  atoms,  in  adequate  concentrations  in  levels  where  hu  (of  the  Infrared)  > E^-E^,  photo- 
lonlzatlon  can  occur  and  as  a result  the  visible  light  output  of  the  plasma  should  be  quenched.  This 
follows  from  the  fact  that  if  destroyed,  the  atoms  excited  Into  those  high  lying  states  are  no  longer 
available  to  feed  the  populations  of  lower  energy  states  from  which  visible  light  originates.  Also, 
since  new  electrons  are  liberated,  the  free  electron  density  should  be  Increased. 

Wa  have,  therefore,  two  Independent  ways  of  testing  this  hypothesis: 

1)  Optical— by  observing  light  quenching.  « 

2)  Electrical — by  measuring  a change  in  the  electrical  conductivity  of  the  plaaaa. 

2.  generality  of  the  phenomenon. 

In  as  auch  as  any  atomic  specie  has  energy  levels  very  near  the  continuum.  It  Is  to  be  expected 
that  the  phenomenon  of  visible  light  quenching  by  the  infrared  photons  will  occur  In  a variety  of  gases. 
While  our  first  experiments  (see  Fig.  1)  were  performed  with  helium  plasaas  and  Infrared  from  a Q-swltched 
C02  laser,  we  have  been  able  to  readily  show  that  afterglows  in  all  the  noble  gases  at  low  pressures  show 
the  same  effect.  Figure  2 illustrates  the  effect. 

It  can  be  readily  understood  that  In  view  of  the  relatively  low  photolonlzatlon  cross  sections, 
rather  high  Intensity  infrared  radiation  may  be  needed  to  observe  the  effects.  The  Q-swltched,  CO^  laser 
Is  operated  at  about  1 kW  of  ln-cavity  peak  power. 

If  high  Intensity  beams  are  used  then.  In  addition  to  single  photon  effects,  multiple  photon  effects 
may  also  occur.  Resonance  excitation  from  a atate  of  lower  to  a state  of  higher  quantua  number  may  occur 
with  no  photolonlzatlon.  Moreover,  at  high  power  the  Infrared  may  Interact  directly  with  the  free-electroa 
gas  of  the  plasma  and  this  may  be  partly  responsible  for  the  phenomena  observed. 

2-1.  Effect  of  CO ^ laser  radiation  on  active  discharge  plasaas. 

Pursuing  the  hypothesis  that  the  Infrared  photons  interact  with  atoms  that  are  excited  to  statea 

cloae  to  the  Ionization  atate,  there  la  no  reason  to  limit  the  plasmas  to  conditions  of  low  electron 

temperature  and  density  ss  in  afterglow  plasmas.  Excited  states  close  to  the  Ionization  limit  muat  be 

populated  contlnuoualy  from  ground  atate  atoms  in  dc  or  ac  discharge  plasmas  In  the  low  pressure  noble 

gaaes.  However,  these  excited  atoms  are  not  In  the  same  concentrations  as  In  the  afterglow  plaamaa  where 

high  principal  quantua  numbered  states  are  populated  In  recombination  of  slow  electrons  with  positive  Ions, 

and  not  froa  lower  states.  Here  again  the  Infrared  radiation  effect  should  be  essentially  Independent  of 
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Che  nature  of  the  monatomic  gas  used. 

In  order  to  test  spaclally  distinct  regions  of  the  plasma  that  differ  In  electron  density,  electron 
temperature  and  associated  excited  atom  populations,  a cold  cathode  discharge  tube  of  the  type  shown  In 
Figure  1 was  used.  Two  essentially  different  regions  are  considered:  a)  the  cathode  region  and  2)  the 

positive  column  region  (nearly  or  fully  developed) . 

Figures  4,  S and  6 Illustrate  the  generality  of  the  effect  of  Infrared  by  showing  quenching  produced 
on  dc  excited  noble  gases.  It  is  seen  that  the  Q-swltched  CO^  laser  radiation  produces  photoeffects  In 
the  visible  region  of  the  spectrum  which  are  readily  monitored  with  photomultipliers.  Thus  there  Is 
possible  application  to  infrared- Image  transformation.  The  plasma  response  Is  rapid.  Figure  7A  and  7B 
show  comparisons  between  the  infrared  signal  as  detected  by  a Ce:Au  Infrared  detector  (which  requires  a 
low  temperature,  77*K,  bath  for  operation)  and  as  seen  as  perturbation  of  the  light  from  a plasma  at  room 
temperature. 

In  view  of  the  Interaction  between  Che  electrons  and  Che  atoms  and  Ions  of  the  host  gas,  any  light 
Intensity  variation  In  the  plasma  will  be  accompanied  by  a perturbation  in  the  electron  gas — directly  or 
Indirectly  Induced  by  the  infrared.  Perturbation  In  the  electron  gas  Implies  change  in  either  Its 
density  or  its  temperature  or  In  both.  In  a dc  discharge,  current  variation  Is  observed  during  Infrared 
Irradiation.  It  occurs  In  the  form  of  an  increase  in  current.  This  current  modulation  is  as  readily 
observed  and  measured  as  the  photoeffect.  Figure  8 illustrates  that  it  is  as  rapid  and  as  faithful.  We 
have  to  keep  in  mind,  however,  the  fact  that  in  comparing  these,  the  current  Is  an  integrated  effect  over 
the  whole  discharge,  whereas  the  visible  light  is  sampled  only  from  a small  volume  of  the  plasma,  not 
all  of  which  Is  traversed  by  the  Infrared  beam.  In  any  event,  these  are  two  independent  ways  of  monitor- 
ing the  CO^-laser  pulses.  It  is  to  be  noted  that  the  cross  section  of  the  infrared  beam  can  be  smaller 
or  larger  than  that  of  the  plasma.  Gaseous  plasmas  are  highly  flexible  media.  They  can  be  dimensioned 
and/or  tailored  to  many  specialized  tasks. 

One  such  task  Is  tbs  use  of  gaseous  plasmas  for  the  detection  of  both  simultaneous  or  nonslmultaneous 
pulsed  microwave  and  CO^- laser  Infrared  radiations. 

2-2 , Pulsed  microwave  and  Infrared  radiations  In  gaseous  plasmas. 

It  Is  now  well  understood  that  microwaves,  which  can  propagate  In  gaseous  plasmas.  Interact  directly 
and  only  with  the  free  electron  gas  of  the  plasma.  The  resulting  physical  phenomena  can  be  used  for 
detection  and  monitoring  of  microwave  pulses. 

Figure  9 shows  the  effects  observed  on  such  plasmas  when  pulsed  microwaves  and  Infrared  radiations 
are  propagated  either  simultaneously  or  one  following  the  other. 

A particularly  Interesting  feature  In  the  pulsed  light  outputs  from  the  plasma  Is  the  fact  that  the 
Infrared  pulse  quenches  the  light  output  regardless  of  the  sign  of  the  light  Intensity  variation  produced 
by  the  microwaves.  Thus,  the  plasma  possesses  a discriminatory  feature  by  which  It  can  distinguish 
between  pulses  of  infrared  and  microwave  frequencies. 
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2-3.  Sensitivity  as  a detector  of  Infrared  radiation  in  the  9 to  11  micron  range. 

Thu*  far  no  atteapts  have  been  made  Co  determine  or  optimize  the  minimum  detectable  Infrared  power. 

2 

Vary  preliminary  estimates  Indicate  that  In  tubes  of  the  order  of  1 cm  cross  section,  and  for  the  plasma 
condltlonsdlscussed,  dc  current  perturbations  provide  a sensitivity  of  the  order  of  30  microvolts  per 
watt  Incident. 

It  has  been  shorn  that  Infrared  radiation  In  Che  10  micron  range  Interacts  with  gaseous  plasmas  of 
monatomic  gases.  The  results  of  this  Interaction  provide  two  Independent  methods  (optical  and  electrical) 
for  the  detection  of  pulsed  C02  laser  radiation.  The  detection  response  time  Is  relatively  rapid 
(<  lusec).  No  restrictions  are  Imposed  on  the  temperature  of  the  gas  hosting  the  plasma.  The  plasma 
medium  can  be  dimensioned  and  shaped  for  particular  tasks.  The  plasmas  described  can  be  used  for 
monitoring  both  pulsed  Infrared  and  microwave  radiations  even  when  these  are  present  in  the  plasma 
simultaneously.  The  study  pursuing  the  Identification  and  quantitative  understanding  of  the  basic 
processes  Involved  In  the  observed  effects  is  being  continued. 

It  is  a pleasure  to  acknowledge  the  Important  contributions  of  Mr.  J.  H.  Hammond  to  this  work. 
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Figure  1.  Simplified  Schematic  Cross-Section  of  Apparatus 


Figure  Z.  Eiomple  of  the  quenching  of  the  visible  otter- 
glow  of  a helium  plasma  by  a single  infrared 
pulse;  Z ^isec/Wiv. 
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Figure  4.  Quenching  of  visible  emission  by  IR  in  o dc-eicitedHe 
discharge)  5/xsec/div,  5mv/div. 

Observations  made  near  the  cathode  neck 
(Discharge  085  To rr  helium,  14mA  in  16mm  dia  tube, 
hollow  cathode ) 
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Figure  5.  Ouenching  ot  visible  emission  by  IR  in  o dc-excitedNe 
dischorge  ; 20^isec/div , 5mv/div 
Observations  made  near  the  cathode  neck. 

(Discharge  < 1.3  Torr  neon,  3.1mA  in  16mm  dio  tube, 
hollow  colhode.) 


Figure  6.  Ouenching  of  visible  emission  by  |R  in  o dc-excited  Xe 
dischorge,  IO/xsec /div.  5mv/div 
Observations  made  neor  the  cathode  neck. 
(Discharge  l7Torr  xenon,  75mA  in  16mm  dio  tube, 
hollow  cathode  ) 


Figure  7 A 

Ouenching  of  visible  emission  (bottom  troee • 5/i.sec/div, 
lOmv/div)  by  IR  (top  troce:  Ge:Au  detector,  S^isec/div, 
20mv/div)  in  a dc-exciled  He  discharge.  Observation 
mode  near  the  cathode  neck.  (Discharge.  0.9Torr  helium, 
14mA  in  16mm  dia  lube,  hollow  cathode.) 


Figure  78 

(Some  os  Figure  7A,  except  2^isec/div.) 


11?;  * r . ; 

& 

t - 

— 

| - . 

\xpwv*  X'  .b  * 

£ 

Light  Current 


lO^xsec/div  S/i  sec/div 

Figure  8. 


Current  perturbation  (top  trace-  IO^i,A/div)  and  quenching  of  plasma 
emission  (bottom  troce  •- lO/xv/divlby  IR  in  a dc-e»eited  He  plosmo. 
Observations  made  near  the  cathode  neck.  (Discharge- l.3Torr  helium, 
7.3mA  in  16mm  dia  tube,  hollow  cathode.) 


Figure  9. 

Illustration  of  the  independence  of  IR  and  microwave  interactions  with  a dc -excited 
plasma.  Upper  trace:  discharge  current;  lower  trace:  plasma  light  emission. 
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Quenching  of  Radiation  from  Ionized  Gases  by  Infrared  Photons* 

J.P.  Kaplafka,  H.  Merkelo,  and  L.  Goldstein 
Gaseous  Electronics  Laboratory , University  of  Illinois,  Urbana , Illinois  61801 
(Received  28  June  1971) 

Visible  afterglow  radiations  from  low-pressure  weakly  ionized  inert  gases  are  partially 
quenched  during  irradiation  of  the  gases  by  the  infrared  produced  with  COj  lasers  (photon 
energies  ~ 0. 1 eV).  The  most  plausible  interpretation  of  the  effect  appears  to  involve  photo- 
ionization or  photoexcitation  of  neutral  atoms  in  excited  states  populated  during  the  process 
of  recombination. 


The  original  observation  of  a partial  quenching  of 
visible  emission  from  a plasma,  induced  by  the  in- 
frared radiation  of  a COj  laser  (operating  near  10 
pm),  was  reported  for  a weakly  ionized  helium 
afterglow  plasma. 1 The  quenching  (as  large  as  50^) 
occurred  on  visible  band  radiations  from  the  di- 
atomic helium  molecule  and  was  interpreted  as 
originating  from  photon-produced  change  of  popula- 
tions of  states  of  a neutral  or  ionic  helium  molecule. 
However,  the  quenching  has  since  been  observed 
with  afterglow  plasmas  of  Ne,  Ar,  Kr,  and  Xe  and 
at  pressures  where  atomic  ion  recombination  is 
known  to  dominate  the  volume  loss  rate.  Addition- 
ally, quenching  of  atomic  lines  of  helium  has  since 
been  observed.  Therefore,  it  is  apparent  that  the 
effect  is  of  much  greater  generality.  Thus  far.  the 
most  plausible  interpretation  is  based  on  photoioni- 
zation or  photoexcitation  of  neutral  species  from 
energy  states  near  the  continuum. 

Quenching  of  visible  afterglow  emission  from  room- 
temperature  Ne,  Ar,  and  Xe  has  been  observed  down 
to  pressures  on  the  order  of  0. 1 Torr  but  has  not 


FIG.  1.  Infrared-induced  partial  quenching  of  total  visible 
radiation  from  a Ne  afterglow  plasma  at  a gas  pressure  of 
1 Torr.  The  time  base  is  1 psec/div  and  the  time  in  the  af- 
terglow is  approximately  600  usee. 


been  observed  for  pressures  beyond  a few  Torr.  As 
in  the  case  of  He,  quenching  in  Ne  and  Ar  was  ob- 
served everywhere  within  the  range  where  the  gas 
temperature  was  varied  (from  77  to  300  °K).  Xe 
was  tested  in  the  range  from  195  to  300  °K  and  ex- 
hibits quenching  throughout  this  range.  The  source 
of  the  infrared  has  been  either  a mechanically 
Q- switched  or  a transversely  excited,  pulsed  C02 
laser.  Figure  1 is  for  the  purpose  of  illustration 
and  shows  a partial  quenching  of  the  total  visible 
radiation  from  a Ne  afterglow  plasma  at  a gas 
pressure  of  1 Torr  on  a l-psec/div  time  base.  Al- 
though most  observations  were  made  with  afterglow 
plasmas,  similar  observations  have  been  made  re- 
cently in  dc  discharges  in  He,  Ne,  Ar,  and  Xe.  But 
in  the  case  of  the  dc-discharge  plasmas,  there  is 
less  evidence  to  support  an  interpretation. 

It  is  expected  that  the  free-electron  densities,  in  the 
afterglow  plasmas  studied,  lie  within  the  range  of 
about  1012  to  loVcm3  during  the  afterglow  times 
when  quenching  has  been  observed.  There  is  always 
a certain  "blind  time"  in  the  early  afterglows  where 
quenching  is  unobservable  or  very  small.  This  blind 
time  extends  approximately  to  a time  where  the 
visible  afterglow  emission  generally  assumes  a 
relative  maximum— a time  generally  somewhere  in 
a range  from  about  50  to  100  Msec  after  the  2- Msec 
active  discharge. 

There  is  evidence  that  the  quenching  is  not  primarily 
the  result  of  heating  of  the  free-electron  gas  by  the 
infrared.  On  the  other  hand,  such  heating  of  the 
electron  gas  by  electromagnetic  radiation  is  easily 
produced  at  microwave  frequencies  and  does  result 
in  quenching.  This  occurs  because  increasing  the 
electron  temperature  reduces  the  recombination  co- 
efficient and  thus  quenches  recombination  emission. 
A short-duration  exposure  to  microwave  radiation  at 
any  time  of  the  afterglow  results  In  an  increase  in 
electron  temperature  and  a decrease  toward  normal 
after  the  end  of  irradiation.  The  time  characteristic 
of  this  decay  toward  normal  is  reflected  in  the  re- 
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FIG.  2.  Theoretical  percentage  of  quenching  of  radiation 
originating  from  states  of  low  quantum  number  vs  the  power 
density  of  the  incident  infrared  radiation  at  two  different 
temperatures  of  the  free-electron  gas.  The  free-electron 
density  is  taken  to  be  lO^/cm3. 

covery  of  the  afterglow  light  and  depends  on  the 
relative  mass  of  electron  and  atom  and  on  the  colli- 
sion rates  between  electrons  and  atoms  or  ions.  It 
is  found  experimentally  that  under  many  conditions, 
when  power  levels  are  adjusted  to  produce  com- 
parable amounts  of  quenching,  the  recovery  of  the 
afterglow  light  after  infrared  irradiation  is  very 
much  faster  than  the  recovery  after  microwave  ir- 
radiation. This  has  been  most  easily  observed  with 
afterglow  plasmas  in  Ne  and  Ar.  As  a result,  the 
cause  of  the  quenching  by  infrared  appears  to  differ 
from  the  cause  of  quenching  by  microwaves. 

The  more  intense  of  the  visible  emissions  originate 
from  energy  states  of  low  principal  quantum  num- 
ber, e.g.,  when  employing  spectral  resolution  while 
observing  quenching  by  infrared,  the  spectral  lines 
that  were  studied  originate  from  states  of  principal 
quantum  number  3,  4,  or  5.  In  many  cases  it  is 
clear  that  populations  in  such  levels  cannot  be  al- 
tered directly  by  means  of  photoexcitation  produced 
by  absorption  of  single  infrared  photons  from  a C02 
laser,  because  the  separation  between  energy  levels 
is  greater  than  the  photon  energies,  i.e.,  greater 
than  ~0. 1 eV.  An  appreciable  amount  of  multi- 
photon photoexcitation  does  not  appear  to  be  a strong 
possibility  at  the  infrared  power  densities  employed. 

Photon  absorption  by  atoms  in  more  highly  excited 
states  is  possible  since  energy  levels  arc  less 
widely  separated.  The  quenching  phenomenon  in  af- 
terglows of  Nc,  Ar,  Kr,  and  Xe  is  favored  at  lower 
pressures  where  "collisional- radiative”  recombina- 
tion* is  believed  to  describe  the  volume  recombina- 
tion process.  In  this  type  of  recombination  there  i6 
generally  little  direct  formation  of  excited  states  of 
lower  quantum  number.  Instead,  the  states  of  lower 


quantum  number  arc  principally  populated  by  cas- 
cade (both  radiative  and  collisionally  induced  cas- 
cade) from  excited  states  closer  to  the  continuum. 
Therefore,  reduction  of  populations  of  more  highly 
excited  states  can  result  in  reduced  populations  in 
lower  levels  and  less  emission  from  these  levels. 

In  order  for  photoexcitation  to  occur  at  a significant 
rate  there  must  be  close  coincidence  between  the 
infrared-photon  energies  of  the  laser  and  energy 
separations  in  the  atoms.  However,  photoionization 
does  not  require  such  coincidences  and  actually  re- 
sults in  the  creation  of  a positive  ion  and  free  elec- 
tron so  that  there  can  be  no  subsequent  cascade  of 
the  atom  until  the  ion  again  undergoes  neutralization. 

In  order  to  investigate  the  credibility  of  the  photo- 
ionization hypothesis,  the  approach  of  Bates,  King- 
ston, and  McWhirter  in  examining  recombination2 
was  adapted  to  a recombining  system  of  atomic  hy- 
drogen ions  and  electrons  in  the  presence  of  infra- 
red radiation.  It  was  postulated  that  all  the  laser 
energy  appears  at  10.  6 p and  loss  terms  were  in- 
cluded to  account  for  the  rates  of  photoionization 
from  excited  states  sufficiently  close  to  the  con- 
tinuum to  be  photoionized,  i.  e. , those  of  principal 
quantum  number,  « = 11,  12,  13,  etc.  Then,  by  com- 
paring the  population  of  a given  state  of  low  quantum 
number,  which  was  predicted  with  the  infrared  in- 
cident at  a given  intensity,  to  the  population  pre- 
dicted with  no  infrared  incident,  the  percentage  of 
quenching  expected  on  all  emissions  originating  from 
the  given  state  of  low  quantum  number  was  obtained. 
By  using  the  form  of  photoionization  cross  section 
for  excited  states  of  hydrogen  derived  by  Menzel  and 
Pekeris,  * estimations  were  obtained  of  the  percent- 
age of  quenching  due  to  photoionization  in  terms  of 
kW/cm2  of  the  incident  infrared  beam.  Some  of  the 
results  are  shown  in  Fig.  2.  Actually,  only  an  upper 
limit  on  the  expected  percentage  of  quenching  is  ob- 
tained in  this  manner  because  the  calculation  gives 
predictions  valid  after  equilibrium  has  been  estab- 
lished during  infrared  irradiations.  It  is  expected 
that  such  is  not  the  case  when  pulses  of  infrared 
power  occur  for  only  short  duration.  It  should  be 
noted  that  under  the  conditions  indicated  in  Fig.  2 
the  estimated  percentage  quenchings  on  emissions 
originating  from  levels  of  n = 3,  4,  5 were  nearly 
the  same  so  only  one  curve  is  drawn.  Quenching  is 
observed  at  power  densities  contained  within  the 
range  on  Fig.  2,  although  quantitative  agreement  is 
not  implied.  It  is  also  of  interest  to  note  that  no 
lines  in  an  afterglow  have  yet  been  observed  which 
show  an  enhancement  rather  than  a quenching. 

Figure  2 suggests  a decrease  in  the  percentage  of 
quenching  as  the  temperature  of  the  free  electrons 
increases.  A decrease  is  observed  when  the  tem- 
perature of  the  free-electron  gas  is  deliberately  in- 
creased for  a short  period  of  the  afterglow,  cover- 
ing the  time  when  the  infrared  is  incident.  The 
method  used  to  selectively  increase  the  electron 
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temperature  has  been  application  of  microwave  ra- 
diation. Moreover,  the  largest  percentages  of 
quenching  for  a given  level  of  infrared  have  gener- 
ally been  observed  with  the  plasma  in  a liquid-nitro- 
gen bath. 

The  observed  disappearance  of  the  quenching  in  af- 
terglows of  Ne,  Ar,  and  Xe  with  increase  in  pres- 
sure can  be  due  to  the  fact  that  the  dissociative  re- 
combination of  the  diatomic  molecular  ions — favored 
as  pressure  increases — does  not  result  in  the  for- 
mation of  atoms  in  excited  states  sufficiently  close 
to  the  continuum  for  photoionization  or  significant 
photoexcitation  to  occur,  if  the  majority  of  such  ions 
are  vibrationally  unexcited.  He  is  the  exception  in 
the  sense  that  quenching  does  occur  at  pressures 
where  recombination  with  a diatomic  molecular  ion 
dominates  recombination  loss.  But  the  observation 
of  quenching  on  visible  bands  of  He2  is  not  surpris- 
ing because  the  very  existence  of  bands  indicates 
that  the  recombination  is  not  directly  dissociative. 

It  would  be  difficult  to  ascertain  whether  or  not  pho- 
toionization takes  place  by  seeking  to  record  an  in- 


crease in  density  of  free  electrons  during  infrared 
irradiation.  The  difficulty  is  that  estimates,  using 
the  recombination  model  referred  to  earlier,  in- 
dicate the  sum  of  populations  in  excited  levels  sub- 
ject to  appreciable  photoionization  is  generally  a few 
orders  of  magnitude  less  than  the  corresponding  den- 
sity of  free  electrons.  Moreover,  a small  increase 
in  free- electron  density  could  occur  even  in  the 
presence  of  photoexcitation  alone,  by  means  of  ion- 
izing collisions  with  free  electrons.  Such  collisions 
normally  play  an  important  part  in  determining  the 
populations  of  highly  excited  states. 

Identification  of  the  origin  of  the  quenching  phenome- 
non requires  further  evidence. 
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This  text  comprises  an  analysis  of  the  atomic  processes  occurring  in 
atomic  hydrogen  and  helium  as  they  pertain  to  the  population  densities  of 
excited  atomic  energy  levels.  Three-body  collisional  ionization  and 
recombination,  collisional  excitation  and  de-excitation  between  discrete 
levels,  radiative  recombination,  dnd  spontaneous  transitions  are  considered 
along  with  the  effects  of  an  external  source  of  infrared  photon  flux. 

Detailed  calculations  of  the  population  densities  performed  both  with 
and  \</ithout  the  infrared  radiation  indicate  a decrease  of  densities  due  to 
photoionization  of  highly  excited  states  by  the  infrared.  From  the  calculated 
densities,  quantitative  predictions  are  made  for  the  percentage  quenching  of 
visible  light  emission  which  originates  from  low-lying  states  in  hydrogen 
and  helium  afterglow  plasmas  when  subjected  to  infrared  radiation  at  a 
number  of  wavelengths . 

Infrared  wavelengths  in  the  range  1.46  um  to  32.8  ym  with  a constant 
2 

intensity  of  1.0  kW/cm  are  considered.  Calculations  for  hydrogen  at  a 
temperature  of  300°  K and  a free  electron  density  of  10'°/cm^  are  performed. 
Helium  calculations  were  done  for  temperatures  of  300°  K,  500°  K,  and 
1000°  I<  with  electron  densities  of  10^/cm^,  10^/cm^,  10^/cm^  and  10^/cm^. 


INTRODUCTION 
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There  is  experimental  evidence  of  the  effects  of  infrared  radiation  at 
10.6  pm  on  gaseous  plasmas  [1],  [2],  [3]  and  [A],  One  such  effect  is  the 
quenching  of  visible  afterglow  light  in  plasmas  subjected  to  an  infrared  laser 
pulse,  indicating  a reduction  of  population  densities  on  lower  excited  atomic 
energy  levels.  At  present,  the  hypothesis  is  that  the  quenching  is  indirectly 
due  to  photoionization  of  highly  excited  states  although  the  possibility 
also  exists  that  the  quenching  may,  to  some  degree,  be  due  to  infrared  heating 
of  the  free  electrons.  Since  an  increase  in  electron  temperature  decreases 
the  recombination  coefficient  and  thus  produces  light  quenching,  quantitative 
predictions  of  quenching  due  to  photoionization  were  needed. 

The  present  text  is  a theoretical  analysis  of  the  atomic  processes 
determining  the  various  population  densities,  and  the  change  in  densities 
due  to  a steady  state  presence  of  infrared  radiation  at  a variety  of  wave- 
lengths. From  the  population  densities  calculated,  percentage  quenching  of 
visible  light  can  be  calculated.  This  is  done  for  atomic  hydrogen  and 
helium. 

The  reaction  rates  governing  the  population  density  of  excited  levels 
for  helium  are  derived  from  those  predicted  for  hydrogen  and  the  calculation 
of  excited  level  populations  follows  the  procedure  presented  in  [5]  with 
modification  made  for  additional  loss  mechanisms  such  as  collislonal 
de-excitation  at  the  walls  of  the  discharge  tube  and  photoionization  by 
infrared  photons. 

8 3 11  3 

Free  electron  densities  in  the  range  10  /cm  - 10  /cm  and  electron 
o o 

temperatures  from  300  K - 1000  K are  considered. 

Explanations  of  the  basic  procedure  for  analysis  of  hydrogen  and  helium 
are  identical  except  where  noted. 
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I.  HYDROGEN  AND  HELIUM  ATOMS  AND  THEIR  AFTERGLOWS 


The  hydrogen  atom  consists  of  one  bound  electron  orbiting  a nucleus  of 
charge  +1.  The  energy  level  diagram  is  shown  in  Figure  1.  Using  the  free 
electron  energy  as  a zero  reference  point,  the  ground  state  energy  is 
-109678  cm  \ with  the  first  excited  level  at  -27419  cm  The  measured 
energy  levels  are  very  closely  approximated  by  the  Bohr  formula  in  which 
degenerate  levels  of  a principal  quantum  number  arc  taken  together. 

The  Bohr  formula  is : 


ENERGY 


109678  -1 

2 — cm 

P 


where  p = principal  quantum  number. 


(1.1) 


Since  the  energy  spread  between  degenerate  levels  is  much  smaller  than  the 
difference  in  energy  between  levels  of  consecutive  principal  quantum  numbers, 
these  levels  are  grouped  together  and  are  considered  as  one  level.  The 
population  of  that  level  is  then  assumed  to  be  equally  distributed  among 
the  degenerate  states. 

The  helium  atom  has  two  bound  electrons  surrounding  a positive  nuclear 
charge  of  +2.  Figure  2a  illustrates  the  ground  state  condition  in  which 
each  electron  is  at  -318066  cm  Figure  2b  shows  the  effects  of  exciting 
one  of  the  electrons,  called  here  e' . As  e'  is  excited  from  the  ground  state, 
e forms  an  orbit  closer  to  the  nucleus  at  -437827  cm  ^ . The  energy  lost  by  e 
is  gained  by  e' , which  then  requires  a total  energy  of  only  -J98305  cm  ^ to 
be  Ionized.  Thus,  since  e remains  in  the  new  ground  state  (actually  it 
moves  slightly  lower  in  energy  as  e'  increases  in  energy)  until  e'  is  fully 
excited  to  the  ionization  energy,  e'  can  be  considered  to  be  part  of  a one 
electron  type  of  atom  with  a ground  state  energy  of  -198305  cm  \ 
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As  e'  becomes  increasingly  excited,  the  interaction  between  the  two 
electrons  decreases  and  the  shielding  of  one  of  the  protons  by  e makes  it 
appear  that  e'  is  part  of  a hydrogen  atom.  Comparison  between  tabulated 
helium  energy  levels  [6]  and  those  for  hydrogen  calculated  by  Equation  (1.1) 
for  p = 4 show  an  error  of  less  than  0.1%.  This  error  decreases  with 
increasing  p.  Furthermore,  the  energy  spread  between  degenerate  levels 
decreases  with  increasing  p,  allowing  the  approximation  to  be  made  that  for 
levels  p > 4,  helium  can  be  considered  basically  hydrogen-like  in  energy 
structure.  Belov;  p = 4,  the  energy  spread  between  degenerate  levels  is 
increasing,  necessitating  the  grouping  of  these  levels  by  angular  momentum 
quantum  number.  Thus,  for  p less  than  four  in  helium  the  levels  are  treated 
as  approximately  hydrogen-like  but  treated  as  separate  levels.  For  p > 4 
the  levels  are  considered  to  be  essentially  hydrogen-like  with  the 
degenerate  levels  grouped  together  by  principal  quantum  number. 

Both  hydrogen  and  helium  afterglow  plasmas  are  treated  in  this  text. 

The  afterglow  plasma  consists  of  free  electrons,  ions,  excited  atoms,  and, 
predominantly,  neutral  atoms  in  the  ground  state,  i.e.,  those  constituents 
remaining  after  the  removal  of  the  accelerating  electric  field  in  the 
discharge.  Once  the  field  is  removed,  the  electron  temperature  decays 
toward  the  gas  temperature  by  thermal  contact  with  the  neutral  gas  until  the 
two  temperatures  are  approximately  equal. 

Principle  recombination  mechanisms  for  the  free  electrons  and  ions  for 


atomic  hydrogen  and  helium  are 


+ * 

X + e + e-*-X  + e 


X + e ■+■  X + hv  . 
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Reaction  (1.2a),  which  is  a three-body  collisional  recombination  into  an 
excited  level, predominates  in  the  early  part  of  the  afterglow  where  the  free 
electron  and  ion  densities  are  high.  The  reaction  of  (1.2b),  in  which  the 
excess  energy  of  the  combining  electron  is  given  off  as  radiation,  predominates 
in  the  later  afterglow.  Thus,  both  reactions  must  be  considered  simultaneously 
as  presented  in  [1],  As  previously  mentioned,  electron  densities  of 
10  /cm  - 10  /cm  and  temperatures  of  300°  K - 1000°  K are  investigated. 
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II.  MKTIIOD  OF  CALCULATION  OF  DENSITIES 

The  calculation  of  light  quenching  requires  precise  predictions  of 
population  densities  of  excited  states.  The  method  used  follows  that  given 
in  [5]  in  presenting  the  collisional-radiativc  recombination  theory  for 
hydrogen  and  hydrogen-like  ions.  The  following  method  is  written  for  singly- 
charged  positive  ions.  It  is  also  restricted  to  situations  in  which  the 
free  electron  density  is  equal  to  the  ion  density.  Furthermore,  in  the 
following,  p and  q refer  to  principal  quantum  numbers  for  hydrogen  and  level 
numbers  in  the  case  of  helium  in  which  levels  2 and  3 are  broken  into  their 
degenerate  levels. 

Let  K(c,p)  be  the  rate  coefficient  for  the  three-body  recombination 
process 

N+  + e + e -*•  N (p)  + e , (2.1) 

where  N*  = the  hydrogen  or  helium  ion 

N(p)  = excited  neutral  atom  in  state  p 
e ■ free  electron. 

For  the  reverse  process  of  ionization,  the  rate  coefficient  K(p,c) 
describes 

N(p)  + e-*-N+  + e + e . (2.2) 

K(p,q)  is  the  rate  coefficient  for  collisional  excitation  or  de-excitation 
as  described  by 

N(p)  + e - N (q)  + e , (2.3) 

and  0(p)  is  the  rate  coefficient  for  radiative  recombination  into  level  p 
described  by 
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N + e -*  N (p)  + hv  . 


(2.4) 


A(p,q)  is  the  spontaneous  transition  rate  for  the  process 


N(p)  -*■  N(q)  + hv  p > q . 


(2.5) 


A uniform  distribution  among  degenerate  states  of  a level  number  is  assumed. 

Assuming  electronic  transitions  due  to  heavy  particle  collisions  are 
negligible,  the  plasma  optically  thin,  and  ignoring  boundary  effects,  the  rate 
of  change  of  population  of  level  p can  be  written  as 


dn(p) 


* -n(p)  [n(c)K(p)  + a(p)]  + n(c)  7 K(q,p) 

q?p 


+ l n(q)A(q,p)  + n(c)  [K(c,p)  + 3 (p) ] , 

q>p 


(2.6) 


where  n(p)  = population  density  of  level  p 

n(c)  * free  electron  density 

K(p)  * K(p,c)  + l K(p,q) 

q^p 

a(p)  * l A(p,q)  . 

q<p 


Denoting  by  n (p)  the  number  density  in  level  number  p in  Saha  equilibrium, 
£• 

letting 


p(p)  = n(p)/nE(p)  , 


(2.7) 


and  using  the  relations  required  for  detailed  balancing 


nE(q)K(q,p)  **  nE(p)K(p,q) 
n(c)K(c,p)  « nE(p)K(p,c)  , 


T 
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Equation  (2.6)  can  be  put  in  the  form 

■ -P(p)[n(c)K(p)  + a(p) ] + £ p(q)n(c)K(p,q) 

E q^p 

n (q)  2 . 

+ £ P(q)  TTT^T  A(q.p)  + n(c)K(p,c)  + S-gf  3(P)  : (2.9) 

q>p  Evp/  Ev*^ 

an  infinite  set  of  coupled  differential  equations. 

For  most  levels  and  typical  values  of  n(c)  and  T, 

n£(p)  <<  n(c)  p ^ 1 

and 

n(p)  <<  n(c)  p ^ 1 . 

Provided  that  the  temperature  is  not  high  enough  to  permit  substantial 
electron  pumping  from  the  ground  state,  the  total  population  density  of 
excited  states  is  much  less  than  the  free  electron  density.  Under  these 
conditions  and  using  the  fact  that  radiative  and  collisional  processes  are 
very  rapid,  the  excited  state  population  densities  quickly  reach  an  equili- 
brium value  before  n(c),  the  free  electron  density,  can  change  appreciably. 
Thus  a given  n(c)  and  T uniquely  determine  the  excited  state  populations. 

This  simplifies  the  problem  because  n(p)  for  all  excited  levels  in 
Equation  (2.9)  can  be  set  to  zero,  resulting  in  a set  of  coupled  linear 
equations  Instead  of  differential  equations.  Furthermore,  when  p is  large 
enough,  collisional  processes  dominate  and  n(p)  can  be  assumed  to  be  equal 
to  the  Saha  density.  Designating  this  critical  level  s,  it  is  then  only 
necessary  to  work  with  s-1  equations,  since  n(l)  is  assumed  known,  thus 
avoiding  an  infinite  matrix.  The  initial  choice  for  s is  somewhat  arbitrary 
and  determines  the  accuracy  of  calculations.  It  is  first  assumed,  then  the 
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equations  are  solved  and  checked  to  sec  that  p( s),  which  is  n(s)/n  (s),  is 

E 

sufficiently  close  to  unity.  The  value  of  s = 28  was  chosen  for  hydrogen  and 
helium.  Adding  also  a term  for  de-excitation  into  the  ground  state  by 
collision  with  the  walls  of  the  discharge  tube  and  loss  due  to  photoioniza- 
tion due  to  infrared  photons,  the  final  form  of  Equation  (2.9)  to  be  solved 


p(p)ln(c)K(p)  + a(p)  + D(p)  + IR(p)] 


_ nF(q) 

- I p(q)n(c)K(p,q)  - £ p (q)  — — A(q , p) 

I^P  <1>P  "E(P) 


-n(c)[K(p,o)  +K(p,c)]  +a^r£i+  I r^Kq.P) 


nE(P) 


<,ts 


(2.10) 


for  p f 1 


p < s * 28  , 


where  K(p,o) 


I K(p,q) 
q>s 

de-excitation  rate  with  walls 


IR(p)  ■ IR  photoionization  loss  rate. 


The  percentage  of  light  quenching  that  occurs  is  calculated  for  the 
levels  from  which  the  spontaneous  emission  is  in  the  visible  wavelength 
range.  For  hydrogen,  these  are  levels  2 through  6,  for  helium,  level 
numbers  11  through  14  which  correspond  to  principal  quantum  numbers  3-6. 
Since  the  light  from  a level  p is  given  by 


I « n(p)A(p,q)hv 


pq  • 


(2.11) 


where  I ■ intensity 


n(p)  » density  of  level  p 


T 


11 


p,q  • levels,  where  q < p 
A(p,q)  ■ spontaneous  transition  rate 


h ■ Planck's  constant 


v • frequency  of  emitted  light, 

P<1 


the  percentage  quenching  of  light  is  proportional  to  the  change  in  population 
density  and  is  given  by 


n(p)|  - n(p) 


x 100  = 


PERCENT 

QUENCHING 


(2.12) 


The  equations  are  solved  by  the  Gauss-Seidel  method  for  solving 
simultaneous  linear  equations.  Consider  an  example  of  three  equations. 


Allpl  + A12p2  + A13P3  ” C1 


A21P1  + A22P2  + A23P3  “ C2 


(2.13) 


A31P1  + A32p2  + A33P3  ~ °3 


arranged  in  the  following  form: 


C1  ' A12P2  “ A13P3 


(2.14a) 


C2  " A21P1  " A23P3 


(2.14b) 


C3  “ A31pl  “ A32p2 


(2.14c) 


It  is  assumed  that  ■ Pj  “ Pj  ■ 0,  where  the  subscript  refers  to  the 


variable  number  and  the  superscript  refers  to  the  number  of  iterations. 


These  are  put  into  (2.14a)  giving  pj,  then  and  p^  are  put  into  (2.14b) 
giving  p*,  and  finally  p*  and  p*  are  put  into  (2.14c)  giving  p*. 

This  process  is  repeated  until 


MK  , max  |pf>  - pf-^l  < IQ'6  , 


i.e. , the  maximum  value  of  the  difference  between  successive  iterations  of 
the  value  of  p for  all  p's  is  less  than  10~6. 
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III.  ATOMIC  COEFFICIENTS 


The  Saha  equilibrium  density  of  excited  states  is  that  density  which  is 
predicted  for  a closed  system  in  thermodynamic  equilibrium.  From  Fowler  [7], 
the  Saha  density  for  an  atom  in  state  p is: 


nE(P> 


n (c)gQ(p) 


h 

2TrmkT 


2 \3/2  h) 

kT  eXpkT, 


(3.1) 


where  n,,(p)  = Saha  density  of  state  p 


n(c)  * free  electron  density 
ggCp)  = degeneracy  of  atom  in  state  p 

g^  * degeneracy  of  ion  energy  level 


h * Planck's  constant 


m * electron  mass 


T * equilibrium  temperature 
k = Boltzmann's  constant 


I * energy  needed  to  ionize  the  atom  from 
^ state  p 

w - 3.1415  . . . 


Hydrogen  has  a degeneracy  of  2p  for  an  atom  with  its  electron  in  level  p, 
and  the  degeneracy  of  the  positive  ion  is  1.  The  ionization  energy  from 
level  p in  inverse  centimeters  is 


„ _ 109678  -1 

I — cm 

P 


(3.2) 


Substituting  these  values  and  the  constants  into  Equation  (3.1)  results  in 


n£(p)  - 4.2  x 10'16  x n2(c)  x p2  x T~3/2  cxp[157859/p2T]  . (3.3) 
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Thus,  given  an  electron  density  and  temperature,  the  Saha  density  for  the 
various  levels  arc  computed. 

The  singly  charged  Hc+  ion  has  a degeneracy  of  2.  The  atom  with  one 

electron  in  the  ground  state  and  the  other  in  state  p has  a degeneracy  of 

2 

4p  if  p is  taken  as  the  group  of  degenerate  states  with  principal  quantum 
number  p.  For  those  degenerate  levels  which  must  be  considered  separately, 
the  degeneracy  is  2J  + 1,  where  J is  the  total  angular  momentun  quantum 
number.  Thus  the  formula  for  computing  the  Saha  density  for  helium  is 

nE(p)  = 1.05  x lO-16  x n2(c)  x gQ(p)  x T-3^2  exp[I(p)/kT]  . (3.4) 


The  values  of  gg(p)  and  I(p)  are  stored  in  the  array  EHE(35,4)  and  are  used 
in  the  program  for  calculating  the  Saha  densities  for  helium  - SAEQP,  which 
is  presented  in  the  appendix.  The  calculated  densities  are  fed  into  the 
main  program  at  execution  time. 

The  compilation  of  spontaneous  transition  probabilities  for  hydrogen 
is  relatively  straightforward.  For  transitions  between  levels  less  than 
twenty-one,  the  tabulated  values  (8]  were  used.  Transitions  from  levels 
greater  than  twenty  were  calculated  by  use  of  the  formula  presented  in  [9] 
for  the  asymptotic  oscillator  strength 


f(p,q) 


_i\”3  _i  _i  _i 

2 3 3 

q / P q *p 


where  p,q 


g 


P 


the  two  levels 

the  degeneracy  of  level  p 


(3.5) 


The  spontaneous  transition  probability  is  then  given  by  [8]  as 
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A(q»p)  - f(Ptq)  , 

* 8q 

where  8p»8q  " degeneracy  of  levels  p,q 

f(p,q)  ■ oscillator  strength 

X ■ wavelength  of  emitted  light. 


With  the  use  of  Equation  (3.5),  and  the  appropriate  conversion  from  X to 
energy.  Equation  (3.6)  becomes 


A(q,p)  “ 


1.29  3q (E  - E )2 

q p 

, 2 M 
(q  - p ) 


q > 20  , 


where  (E  - E ) = energy  difference  of  levels  q and  p 
q P in  cm'1. 


Helium  presents  a difficult  problem  due  to  the  scarcity  of  published 
values  for  transitions  between  highly  excited  states.  For  transitions  to 
level  numbers  eleven  and  below  from  higher  levels  up  to  fourteen  or  twenty, 
published  values  [8]  and  [10]  were  used.  Transitions  to  these  levels  from 
higher  excited  states  were  calculated  with  the  use  of  a form  for  the 
asymptotic  oscillator  strength  as  presented  in  [11] 


f(Piq)  * — • 


This  was  then  substituted  into  Equation  (3.G)  and  normalized  for  each  level 
to  the  rate  for  the  highest  level  published. 

This  resulted  in  the  following  equations: 


A(q,p)  - C x g x (E  - E)  /q' 

P P q P 
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C 


A* (q,p)q 

»P<E,  - y: 


(3.10) 


where  Cp  - normalizing  constant  for  level  p 
A' (q,p)  * rate  from  the  highest  of  known  levels 


and  gp. 


E and  E are  as  previously  defined. 
<1  P 


The  calculated  values  of  C are  given  in  Figure  3. 

P 

The  spontaneous  transition  rates  between  higher  excited  levels  were 
assumed  to  be  approximately  the  same  as  those  for  hydrogen.  Use  was  made 
of  the  published  and  asymptotic  values  previously  mentioned  for  hydrogen 
atoms.  The  rates  for  helium  were  tabulated  and  calculated  and  then  stored 
by  use  of  the  program  SPONA.  This  program  punches  the  rates  onto  cards  for 
use  by  the  main  program  and  is  presented  in  the  Appendix. 

The  radiative  recombination  coefficient  describes  the  rate  at  which 
free  electrons  are  captured  into  excited  atomic  states.  Bates  [5]  gives 
the  general  formula  for  hydrogen  as 


w f/2  exp  (I  /kT)  V 

f5(a,T)  = — a 2 1/2  3’/"2  f K ’ 

i U V ' (mkT)  ' J 


(3.11) 


where  w , wJ 
a’  i 


w_  f/2  exp(Ia/kT) 

(mkT) 

degeneracy  of  the  atomic  state  a,  and  the  ion 


I 

a 


K 


with  hv 


k 


ionization  energy  from  state  a 

00 

/ (hv)2oal(v)  expj^jrj.  d(hv)  , 

*a 

energy  transferred  to  electron 
photoionization  cross  section  for  state  a 

Boltzmann’s  constant. 
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l 


■ 


: 


m ■ electron  mass 
h ■ Planck's  constant 
T ■ temperature 
c - speed  of  light. 


Seaton  [12 J gives  the  photolonlzatlon  cross  section  as 

-6  2 

2 aira  n 2 _3 

o ..(v)  * ZT~  (1  + n e)  8i  1 (n>e)  . 

nl  3/3  11 

with  n = principal  quantum  number  in  hydrogen 
o ■ fine  structure  constant 
Sq  * Bohr  radius 

g^3  * Kramers-Gaunt  factor 

e * excess  energy  of  ejected  electron. 


(3.12) 


Using  (3.11)  and  (3.12)  for  hydrogen,  where  w /w.  = 2n  results  in 

3 X 


B(n,T) 


DX 


1/2 


n n n 


x S (X)  * 


(3.13) 


18 


with  the  effective  quantum  number  resulting  In 


1/2 

Dg*  ' x S<X) 

B(p,T)  - P 

4"eff 


(3.14) 


where  now,  p » level  number 


gp  - degeneracy  of  level  number  p 


■ effective  quantum  number  for  level  p. 


Seaton  [12]  gives  an  asymptotic  expansion  for  g 


11 


gn(n,e)  = 1 + 0.1728  n*2/3(U  + 1)_2/3(U  - 1) 

- 0.0496  n'4/3(U  + 1)_4/3(U2  + 4U/3  + 1)  + . . . (3.15) 


Substituting  (3.15)  into  S (X)  results  in 

n 


§ 


S (X)  = S(0)(x  ) + X_1/3S(1)(x  ) + X_2/3S(2)(x  ) + . 


(3.16) 


where  S^(x)  = exp(x)  E^x) 
S(I)(x)  = 0.1728  Xx(x) 
S(2)(x)  = -0.0496  x2/3  X2(x) 


and 


“>  -v 
e 


Ei(x) " / V dv 


X (x)  - / (U  + 1)'5/3(U  - 1)  e_xU  dU 
0 


<30 

X- (x)  - / (U  + 1)"7/3(U2  + 4U/3  + 1)  exp(-xU)  dU  . 
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For  0.02  < xQ  < 20,  values  for  xS^\  xS^\  and  xS^  were  interpolated 

from  Seaton's  table  of  calculated  values.  Values  of  x >20  were  calculated 

n 

by  his  approximate  formulas 


xS 


xS 


xS 


(0) 


(1) 


(2) 


1 + 1/x  + 2!/x2  + 3!/x3  + 4!/x4 


-0.1728  x 


1/3 


-0.0496  x 


2/3 


N 

M 


70 


(3x)" 

32 


800  + 11440 


(3x) 
448 


(3x) 


(3x)‘ 


(3x) 


- 

i 


(3.17) 


The  values  of  8(p,T)  for  hydrogen  were  calculated  in  the  program  BETA  as 
shown  In  the  Appendix  and  then  by  comparison  with  (3.14),  converted  to  give 
those  for  helium. 

Gryzinski  [13]  gives  as  the  cross  section  for  collisional  ionization  of 
an  atom: 


Q(p,c)  = / fi(vo)Q'(U)  dv  , (3.18) 

0 e e 

where  Q(p,c)  * cross  section  for  electron  in  orbit  p 

f*(v  ) » velocity  distribution  of  the  orbital 
electron 

Q' (U)  * Gryzinski’s  derived  formula 
and,  fA(ve)  « 26[ve  - (2E1/m)1/2]  (3.19) 

where  E^  *>  kinetic  energy  of  orbital  electron 
m **  electron  mass. 

He  gives  Q' (U)  as 

-14 

Q'(U)  = -,-56-  j-  10 g[E./IJ,  E./U)  , (3.20) 

U 
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where  U - ionization  energy  from  the  level  p 

E,  • kinetic  energy  of  orbital  electron 
in  p 

" energy  of  incident  electron 


I 


and, 


«tVV 


2 \3/2, 


,2.2 

'V1  + v2; 


t ) / \? 

'f,  E U E,  /U  P")  * 

J2  _1  + _R  i _ _1  _ _H  l for 

l3E2  E2\  E2 J \hj  f " + h <- 


IVj  + v2 


21  ins  e2  e2 


ill 


(3.21) 


In  plasmas,  is  not  a single  value  but  is  a distribution  which  is 


assumed  to  be: 


f (e)de 


2ir 


(*kT) 


3/2 


1/2 

e exp(-e/kT)  de 


(3.22) 


which  is  the  Maxwell-Boltzmann  distribution  of  electron  energies  for  tempera- 
ture T.  Simplification  of  equations  results  when  the  following  change  of 
variables  is  performed 


X - E2/Up  . 


(3.23) 


With  the  additional  formula  for  conversion  of  the  cross  section  into  rate  of 
ionization 
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f i 
I ' 


K(p,c)  - / Q(p,c)v  f(c  )de  (3.24) 

0 Z 

where  K(p,c)  - rate  of  ionization  from  level  p 
Q(p,c)  ■ Equation  (3.18) 

■ velocity  of  incident  electron 

■ Equation  (3.22)  , 

Equations  (3.18)  through  (3.23)  result  in  the  following  equation  as  presented 
in  [5]  for  the  rate  of  ionization  from  level  p; 

* _xu  /kT 

K(p,c)  = 5.45  T'3/2  / I e P dx  , (3.25) 

1 

where  (x  - l)3/2/(x  + 1)3/2  x < 2 

(5x  - 6)x1/2/3(x  + 1)3/2  x > 2 
T ■ electron  temperature 
11^  * ionization  energy  of  level  p 

k ■ Boltzmann's  constant  . 

For  helium  the  same  formula  was  used  with  the  p referring  to  level 

number  and  using  the  appropriate  U . 

P 

Collisional  excitation  between  discrete  levels  is  calculated  in  a 
similar  manner.  In  this  case,  the  cross  section  for  excitation  from  level  p 
to  level  q is  given  by 

m 

Q(p.q)  - / f(ve>  Q'(Uq+1.Uq)  dve  (3.26) 

where  Q' (U  ,U  ) 3 Q' (U  ) - Qf (U  ) . 

q+1  q q q+1 

t 
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The  ocher  symbols  are  Che  same  as  for  collisions!  ionization.  In  this 

instance,  however,  the  classical  cross  section  given  by  (3.26)  must  be 

-3 

modified  by  multiplying  it  by  q to  give  a quantum  related  cross  section 
as  shown  by  McCarroll  [ 14 3 > since  excitation  between  discrete  levels  is 
alien  to  classical  mechanics. 

With  this  modification,  the  use  of  Equations  (3.19)  through  (3.23), 
and  the  Bohr  formula  for  hydrogenic  energy  levels.  Equation  (1.1), 

Equation  (3.26)  becomes  that  given  in  [5]  for  the  collisional  excitation 
rate  between  levels  p,q: 


_2  _2, 


\ _ 10.9  p r _ . f- 15 7 890 (g  - p )y1  , 

K(p»q)  - 0/9  2 2 / ^ exPi  2 2 if 

T ' q(q  - P ) 1 1 P q T J 


(3.27) 


{2  - a + y(l  + 4a) } (y  - l)1/2(a  + 1)1/2 


- 1/2.  T x3/2 
3a  (a  + y) 


y < 1 + a 


(-3a  + y(3  + 4a)}y 
3 (a  + y)3/2 


1/2 


y > 1 + a 


Jl.,-2  _2. 

a ■ q / (q  - P ) • 


where  in  hydrogen. 


p - p 

principal  quantum  numbers  (3.28) 

q - q 


and  in  the  helium  approximation, 

p,q  ■ level  numbers 

p,q  ■ effective  principal  quantum  numbers 

corresponding  to  p,q.  (3.29) 

Collisional  Ionization  and  excitation  rates.  Equations  (3.25)  and  (3.27), 
are  evaluated  by  the  Gaussian  integration  method.  The  integrals,  by  change 
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of  variable,  are  put  in  the  form 

1 

/ f(x)dx  . (3.30) 

-1 


They  are  approximated  by  the  sum 


l a f (x  ) , (3.31) 

i=l 


where  the  x/s  are  the  zeros  of  the  Legendre  polynomial  of  order  n.  The 
a^'s  satisify 


a^  + . . 


+ a = / dx  = 2 
n -i 


Vi + 


+ a x = f xdx  = 0 
n n ^ 


n-l  , 

Vi  + • 


n-l 

+ a x 
n n 


= / 


n-l , 
x dx 


-1 


n even 


— n odd 


(3.32) 


The  tabulated  values  of  x^  and  a^  for  n = 16  are  used  in  the  program  for 
calculating  K(p,c)  and  K(p,q).  This  program,  found  in  the  Appendix,  stores 
the  resulting  rates  in  COI.IZ  and  COLEX  which  are  punched  onto  cards  and  read 
into  the  main  program. 

Collislonal  de-excitation  rates  are  calculated  by 

nE(q)K(q,p)  = nE(p)K(p,q)  , (3.33) 


47 


as  required  by  detailed  balancing,  where  n(q),n  (p)  = Saha  density  of  levels 


The  loss  from  level  p due  to  photoionizatlon  is  given  by 


where  n(p)  = density  of  level  p 


o(p,v)  = photoionization  cross  section 


The  equations  for  I and  a(p,v)  are  given  by 


where  P = power  density  of  infrared 


X = wavelength  of  infrared  radiation 


h = Planck's  constant 


Hydrogen 


where  g = gaunt  factor  - 1 


e = electron  charge 


Rydberg 


1 = ionization  energy  for  level  p 


v = frequency  of  infrared  radiation 


as  given  in  [15]  for  hydrogen  and  hydrogen-like  ions.  Equations  (3.35)  and 


(3.36),  for  a constant  power  density  of  1 kW/cm 


simplify  to 


25 


_ „ „ 5.25  s 1031  x A4  , -lv 

Io(A,p)  * ; — = loss  rate  (sec  ) 


(3.37) 


The  value  of  the  lowest  level  photoionized  is  given  by 


where  U = ionization  energy  for  hydrogen 
ground  state 

X * wavelength  of  radiation. 


(3.38) 


Equation  (3.38)  was  solved  for  the  particular  wavelength  necessary  to  ionize 

levels  4,  5,  6,  ...  , 20.  For  example,  a X slightly  smaller  than  15.39  urn 

will  ionize  levels  13  and  above.  A wavelength  slightly  larger  will  only 

ionize  levels  14  and  above.  In  this  way  calculations  were  made  for  X from 

1.46  pm  to  32.88  pm.  The  wavelengths  and  corresponding  p are  stored  in 

the  main  program.  Since  p for  X > 1.46  pm  is  > 4,  which  has  approximately 

min  ~ — 

the  same  ionization  energy  in  both  hydrogen  and  helium,  the  same  values  are 
used  for  both. 
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NORMALIZATION  CONSTANTS 
FOR  HELIUM  A(q,p) 
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level  number 
P 

actual 

level 

c 

P 

1 

Is2 

46.5 

2 

2s  3S 

9.54 

3 

2s  1S 

51.9 

4 

2p  3P 

16.6 

5 

2p  1P 

44.1 

6 

3 

3s  S 

1.27 

7 

3s  1S 

40.2 

8 

3p  3P 

17.7 

9 

3d  3D 

9.41 

10 

3d  XD 

27.76 

11 

3p  XP 

53.5 

Figure  3. 
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IV.  ARRAYS 

The  array  EUE(35,4)  is  used  in  Che  calculation  of  most  of  the  atomic 
coefficients  for  helium.  It  contains,  for  35  levels:  the  energy  difference 

between  the  level  and  the  ground  state,  ionization  energies,  an  "effective" 
principal  quantum  number,  and  the  degeneracy  of  the  level.  A diagram  of  the 
array  is  shown  in  Figure  4.  The  rows  denote  the  level  number;  the  correspon- 
dence between  the  level  number  and  the  actual  level  in  helium  is  shown  in 
Figure  5.  As  shown  in  that  figure,  levels  one  through  eleven  are  the  first 
eleven  degenerate  levels,  with  twelve  through  thirty-five  grouping  the 
degenerate  levels  together.  Level  12,  for  example,  contains  the  4s,  4p, 

4d,  and  4f  energy  levels. 

Column  one  of  EHE  contains  the  energy  difference  between  the  level  and 
the  apparent  ground  level  of  198305  cm  \ The  first  eleven  levels  were 
obtained  from  tabulated  values  [6];  the  rest  were  calculated  using  the 
hydrogen  approximation 

EHE (p,l)  = 198305  - 109678/ (p  - 8)2  , (4.1) 

where  p denotes  level  number,  not  quantum  number. 

Column  two  holds  the  energy  needed  to  ionize  from  that  level  number 
and  was  calculated  by 


EHE (p , 2 ) = 198305  - EHE(p.l)  . (4.2) 

Column  three  contains  an  "effective"  principal  quantum  number.  This 
makes  the  calculation  of  some  coefficients  which  depend  on  energy  level  easier 
to  calculate,  and  is  the  principal  quantum  number  that  an  energy  level  would 
have  in  a hydrogen  atom 
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ElIE(p, 3)  = [l09678/EHE(p,2)]1/2  , (4.3) 

where  109678  cm  1 = Ionization  energy  of  hydrogen. 

Column  four  contains  the  degeneracy  of  the  level  number.  For  one 
through  eleven  the  degeneracy  is 

EHE (p,4)  = 2Jp  + 1 1 < p < 11  , (4.4) 

where  J is  the  total  angular  momentum  quantum  number  for  the  level.  Eleven 
and  up  are  groups  of  degenerate  levels,  and  for  helium, 

EHE(p,4)  = 4 (p  - 8)2  , (4.5) 

since  a level  number  greater  than  eleven  corresponds  to  a principal  quantum 
number  eight  less.  Array  EHE  is  produced  in  the  program  for  calculating 
spontaneous  transition  rates  and  is  presented  in  the  Appendix. 

Figure  6 shows  the  calculated  energy  levels  for  atomic  hydrogen  using 
the  Bohr  formula 


E (p)  = 109678/p2  . (4.7) 

The  values  of  the  Saha  densities  are  stored  in  the  array  SAEQP (35 ,4, 7) . 

The  first  subscript  refers  to  level  number  and  the  second  to  temperatures  of 
300°  K,  400°  K,  500°  K and  1000°  K.  The  third  refers  to  free  electron 
densities  of  10^,  5 x 10^,  10^,  5 x 10^,  10^,  5 x 10^  and  10^  per  cubic 
centimeter. 

The  rates  of  colHsional  ionization,  K(p,c),  are  stored  in  array  C0I.IZ (30,4) . 
This  array  stores  the  rates  for  the  first  thirty  levels,  and  the  four 


temperatures. 
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The  array  C0LEX(30, 35,4)  contains  the  rates  of  collisional  excitation, 
K(p,q),  where  q > p.  The  first  subscript  refers  to  the  lower  level,  the 
second  to  the  upper  level,  and  the  third  to  the  four  temperatures. 

GSIRL(30,2)  contains  the  infrared  loss  term  for  the  equations  for  a 
possible  use  of  up  to  thirty  levels.  The  second  subscript  refers  to  the  two 
cases:  1 for  no  incident  infrared  in  which  the  loss  term  is  zero,  and  2 for 

the  case  with  infrared  radiation.  This  array  is  calculated  in  the  main 
program. 

SPONA(35,35)  stores  the  spontaneous  transition  rates,  A(p,q),  between  the 
first  thirty-five  levels. 

The  array  BETA(29,4)  holds  the  radiative  recombination  rates  for  twenty- 
nine  levels  for  the  four  temperatures. 

ELECD(7)  and  TEMP (4),  respectively,  hold  the  values  of  the  seven  electron 
densities  and  four  temperatures. 

RLAMB(17,2),  in  the  first  column,  stores  seventeen  wavelengths  of 
infrared  radiation,  and  in  column  two,  the  minimum  quantum  level  effected  by 
the  radiation.  This  information  is  used  to  calculate  the  infrared  loss  terms 
in  the  rate  equations. 

COLDEX(35, 30)  stores  the  values  of  collisional  de-excitation  rates, 

K(p,q),  where  p > q.  These  values  are  calculated  from  COLEX  and  SAEQP  by 
Equation  (3.33). 

Arrays  are  also  used  to  calculate  and  store  the  coefficients  of 

Equation  (2.10).  AE(30)  contains  a(p)  = £ A(p,q),  where  q is,  as  before, 

q<p 

the  level  above  which  the  atoms  are  assumed  to  be  in  Saha  equilibrium. 

RKAPA(30)  holds  K(p)  = K(p,c)  + } K(p,q),  which  is  the  depopulation  rate 

qrp 

for  level  p by  collisional  processes. 

nF(q) 

ASUM(30)  contains  V — A(q,p),  and  SIC.MK(30)  stores  K(p,o)  2 > K(p,q). 

q>s  nE(p)  q>s 
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These  atomic  parameters  are  then  arranged  to  form  the  coefficients  of  the 
rate  equations  and  stored  in  GS(29,29,2)  with  the  last  subscript  denoting  one 
set  with  no  infrared  and  one  with  infrared.  The  constant  terms  [the  right-hand 
side  of  Equation  (2.10)]  are  stored  in  CSC(29,2).  RHO(29,2)  holds  the  values 
of  p and  QUEN(30)  stores  percent  light  quenching  on  each  level  for  printout. 
IITER(2)  and  EEPS(2)  store  the  number  of  iterations  performed  and  the  final 

g 

value  of  M in  Equation  (2.15). 


EIIE(35,4)  ARRAY 
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ATOMIC  HYDROGEN  ENERGY  LEVELS 
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V.  RESULTS  AND  CONCLUSIONS 

Figures  7-10  are  the  tabulated  results  of  the  computer  analysis  for 
atomic  hydrogen  and  helium.  Only  one  case  of  atomic  hydrogen  was  run: 

T ■ 300°  K and  n(c)  = 10^/cm\  it  was  run  primarily  to  demonstrate  the 
feasibility  of  the  method  of  analysis  and  the  programs.  For  atomic  hcLium,  the 

following  cases  were  calculated:  T = 300°  K,  500°  K,  1000°  K and 

8 3 9 3 10  3 11  3 

n(c)  ■ 10  /cm  , 10  /cm  , 10  /cm  and  10  /cm  . As  can  be  seen  from  the 

resultant  values  of  p(s)  which  should  asymptotically  approach  1.0,  the  cases 
8 3 9 3 

of  n(c)  = 10  /cm  and  10  /cm  are  probably  less  accurate  than  the  other  cases 

since  p(s)  for  these  cases  was  approximately  0.4  and  0.8,  respectively.  This 

indicates  that  more  levels  need  to  be  considered  in  the  equations.  In  these 

figures,  the  minus  designates  the  results  for  X slightly  lower  than  that 

indicated  and  the  plus  slightly  higher. 

The  tabulated  results  are  shown  graphically  in  Figures  11  - 23.  As  seen 

from  the  graphs,  the  percentage  light  quenching  increases  with  increasing 

4 

wavelength,  since  the  loss  term  due  to  ionization  is  ^ X , until  the  photon 
energy  decreases  to  a value  which  photoionizes  one  less  level.  At  this  point 
there  is  a sharp  drop.  This  is  repeated,  with  the  average  quenching  increasing, 
until  the  effect  of  the  decreasing  number  of  levels  photoionized  begins  to  pre- 
dominate over  the  increasing  value  of  cross  section,  and  the  average  level  of 
quenching  begins  to  decrease.  The  results  obtained  provide  a basis  for  further 
quantitative  analysis  of  light  quenching  due  to  photoionization  of  highly  exci- 
ted states  and  also  for  analysis  of  the  possible  role  of  electron  heating. 

The  present  analysis  provided  insight  into  the  atomic  processes  occurring 
and  also  pointed  out  important  considerations  for  future  work.  First,  the  solu- 
tion of  the  simultaneous  equations  must  be  performed  in  double  precision 
arithmetic  due  to  the  large  number  of  calculations  involved  in  the  Iteration 
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Second,  more  equations  must  be  used  for  the  lower  electron  densities  so 
that  p(s)  approaches  1.0  more  closely.  This  insures  that  the  upper  levels 
are  very  nearly  in  Saha  equilibrium  with  the  free  electrons,  one  of  the  basic 
assumptions  of  the  method. 

Third,  care  must  be  taken  when  manipulating  the  equations  to  be  solved. 
Operations  which  are  permitted  when  solving  equations  exactly  such  as  multiplying 
each  term  in  an  equation  by  a constant  greatly  influence  the  accuracy  of  the 
results  and  the  efficiency  of  the  program  when  solving  the  equations  by  an 
iterative  method.  This  becomes  increasingly  important  as  the  number  of 
equations  increases. 

Finally,  more  cases  of  electron  density,  temperature  and  wavelength  must 
be  analyzed  along  with  different  values  of  infrared  power  density  to  obtain 
a more  complete  understanding  of  possible  trends  in  the  results. 
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HELIUM  PERCENT  QUENCHING 


5.82  .16 


p(s)  - 0.35  p(s)  = 0.  74  p(s)  = 0.92  p(s)  = 0.97 


n(c)  - 10 
Um  - + 


1.46  .04  .02  .04  .02  .04  .02  .04  .03 


2.28  .12  .05  .13  .06  .14  .07  .16  .09 


3.28  .22  .02  .26  .05  .31  .10  .39  .18 


4.46  .06  .05  .18  .17  .35  .34  .63  .60 


.50  .49  .99  .97  1.71  1.58 


7.38  .39  .38  1.25  1.241  2.45  2.40  3.94  3.32 


9.11  .84  .83  2.74  2.71  5.28  5.15  7.29  5.42 


10.6  1.44  1.44  4.66  4.66  8.79  8.79  9.41  9.41 


11.0  1.64  1.62  5.29  5.24  9.95  9.46  10.72  6.96 


13.11  2.87  2.85  9.16  9. 


15.39  4.49  4.44  14.03  13.74  23.56  19.68  13.70  7.55 


17.85  6.27  6.19  19.15 


20.49  7.95  7.82  23.82  22.6!  31.53  22.94  11.17  5.70 


23.32  9.34  9.13  27.48  25.42  31.70  21.37  9.19  4.49 


26.32  10.30  9.92  29.73  26. If  29.33  17.80  7.06  3.20 


29.51  10.79  10.01  30.14  23.9'  24.56  12. 4S  4.95  1.94 


32. 88  10.64  8.58  27.71  16.5 


Figure  8 
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HELIUM  PERCENT  QUENCHING 
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HELIUM  PERCENT  QUENCHING 
1000°  K 
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PERCENT  QUENCHING  VS 
WAVELENGTH  IN  HYDROGEN 
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Figure  11. 


PERCENT  QUENCHING  VS 
WAVELENGTH  IN  HELIUM 

TEMP=300  °K  ELECTRON  0ENSI7Y=106  /CM3  INFRRREO  JNTENSITY=1KW/CM2 
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Figure  12. 
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Figure  14. 
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PERCENT  QUENCHING  VS 
NRVELENGTH  IN  HELIUM 

TEMP=300  °K  ELECTRON  0ENSm=10n /CM3  INFRARED  JNTENSI7Y=IKW/CM2 
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Figure  15. 


PERCENT  QUENCHING  VS 
WRVELENGTH  in  helium 

7EMP=500  °K  ELECTRON  0ENSI7Y-108  /CM3  INFRARED  INTENSITY=1KW/CM2 
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Figure  16. 


PERCENT  QUENCHING  VS 
WAVELENGTH  IN  HELIUM 

7EMP=S00  °K  ELECTRON  0ENSm=109  /CM3  JNFRflREO  INTENSJ7Y-1KW/CM2 
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Figure  17. 
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PERCENT  QUENCHING  VS 
WAVELENGTH  IN  HELIUM 
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Figure  18. 
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PERCENT  QUENCHING  VS 
WAVELENGTH  IN  HELIUM 

TEMP-500  °K  ELECTRON  0ENSiTY=10n /CM3  INFRRREO  INTENSITY-tKW/CM2 


WAVELENGTH  - 10'b 

Figure  19. 
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PERCENT  QUENCHING  VS 
WAVELENGTH  IN  HELIUM 

TEMP=100QPK  ELECTRON  0ENSJTY=108  /cm3  JNFRflREO  INTENSITY=1KW/CM2 
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Figure  20. 
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PERCENT  QUENCHING  VS 
WRVELENGTH  IN  HELIUM 

TEMP-lOQQPK  ELECTRON  OENSITY^IO10 /CM3  INFRARED  INTENSITY=IKW/CMZ 


Figure  22. 
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PERCENT  QUENCHING  VS 
WAVELENGTH  IN  HELIUM  1 

TEMP*tOOO°K  ELECTRON  0ENSnY=10u /CM3  INPRRREO  INTENSITY=IKW/CM2 


WAVELENGTH  - 10'D  METERS 


Figure  23. 
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The  following  programs  are  those  used  in  calculating  the  rate  coefficients 
and  percent  quenching  in  hydrogen  and  helium. 

The  first  two  programs  are  the  main  programs  for  calculating  densities 
and  quenchings  for  hydrogen  and  helium.  The  third  calculates  the  radiative 
recombination  rates  for  hydrogen  and  the  fourth  takes  these  hydrogen  results 
and  converts  them  to  helium  rates.  The  fifth  program  calculates  collisional 
excitation  and  ionization  rates  for  helium.  The  sixth  program  calculates  the 
Saha  densities  for  helium,  and  the  seventh  program  calculates  the  spontaneous 
transition  rates. 
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calculation  of  excited  state  POPULATIONS  ANO  PEPCENT  OUENCHINO 
FOR  HYOROGEN 

CONSTANT  INTENSITY  l XW  / CK**2 

VARIABLE  NUMBER  OF  LEVELS 

NOLEV  IS  THE  NUMBER  OF  LEVELS  USED  MAXIMUM  LEVEL=29 


01  KENS  I ON  SAEOPIsQ)  ,CCLlZ(30i  .COLEXl  30.'>0)  ,0S  IRL I 30 , 2 ) , AE < 30) , 
ISPOriAl  35.35)  »RETA(2»30).°:<A3A(33I  . S I tV'R  ( 3 .3 ) , a SUM  t 3 3 1 , GS  ( 29 , 20 , 2 ) . 
XGSCI 29.2) »RHG( 29.2) .OUENI 30 J , 1 1 TER  1 2) . EEP S ( 2 ) , ELECUI A ) , TEN? ( 2 1 , 

IRL  AM  BUS,  2 ) .COLOcXI  AO,  30)  .EI35J 

OCUBLE  PRECISION  T , P , SUM , 3, .R , El , E2 . FI . F 2.  C ( 1 6 J , X ( 16  ) , OSORT, OSXP 

DATA  S AEQ?  »COLI 7 «COLEX, CCI.OEX, GS I RL  » AS. SP3NA , BETA ,R  < APA , S 1 CMX , ASUM 
A, GS, CSC,. RHO.OUEN.tE PS/  5769*0. 0/. 1 1 TER/  2*0/ . EPS/ 1 . OS-6 / , MAX/ 350/, 
5ELECO/5.0L3, 1 . OC-9. 5 .0E9 , l .06 10/  , T E M P / 3 00  • 0 » 9 00  « 0/ » 

6f  PL  AMR  ( I . l)  , f ~ 1 . 1 6 J / I .96E-6.2.23S-5,  3 . 2SE -6 , 4 .46E-6 , 5. 82E-6 , 
7T.33E-6.9.11E-6,  U.E-E,  < 3.  I IE-i  . I 5 .39E-6  , l 7 .35E-6  , 20. 495-6, 

823. 32 E- 6.  26.32E-6, 29.5 XE-S . 32.  S8E-5/ , ( F I.A M3 < l , 2 > . I = l . X 6) /4.  , 5.  , 
96. ,7. , G . »9. ,10. ,11. ,12. ,X3.,19.,X5.«X6.,X7.,X8.,19./ 

IT  = X =300  2=900 

10=1=5. E8  2=1. E9  3=5. E9  9=1. ElO 

IL=l=l .06  2=3.39  3=5.0  9=10.6  5=28.0 

READ  IN  TEMP ,OEN, LAMBDA  QESIRFO 

RCAO.IT.IO.IL 
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10*4 
I T * l 
REAO.E 
C 

c 

C READ  IN  VALUfcS  OF  HCTA(2,25) 

C 

READ.IflETAI  l, I). 1=2, 19) 

PE£0, (HETAI2, [) ,1=2,191 
READ, (BETA  (1,1),  1=2').  26  ) 

REAO. ( BETA!  I , I ) , 1=26, JO) 

C READ  IN  VALUES  Of;  SPONA(20.20) 

00  l 1=1.19 
11=1*1 

1 REAO.ISPCNAIJ.I)  ,J=II,20) 

C CALCULATE  HIGH  VALUES  OF  SPGNNA 
DO  501  J=2l,35 
JJ=J-1 

DO  501  K = 1 , J J 

501  SPONAI J,lC)s<K«*3l*J«( ( E ( J )-E I K ) )*«2)/(  ( (J**2)-(K*»2) ) **  3 ) 

C 

C 

Cl  l)=. 0 2715245 04 1175400 
C 1 2 ) = .062253  52393364800 
C(  3 ) = . 095  1505  l 1602'.93D0 
C(4)=. 12462097125553400 
C( 5)*. I 49 595 9 8 3 8165/7 00 
C( 6)  = .169  1 5 6 51 9395003 00 
C( 7) =.132603 41504 49 2400 
Cl «)=. 1694504 1 945506900 
CI9)=. 10 94 506104550 69 DO 
C<  lu)  * . 1 8260  J4- 1 5 04  492400 
Cl  11)  = . 169156619396 003D0 
C 112 )=. 14 95959308 1657 7D0 
Cl  13)  = . 1246289  71 255  5 3->D0 
Cl  I*)  = .39515351)68249300 
Cl  15  ) = . 05225  362393364(100 
Cl  It  ) = . 027152656411 75600 
XI  1)=-. 9894019 >499165300 
X I 2 ) =- . 944 5 75023073  23  300 
X(  3)=-.  365 1-312 02 33733200 
X ( 4)  = -.  755404406  35500  >00 
X ( 5 1 *- .61  787624440264400 
XI 5) =-.4560167  7 7o5  722  700 
XI 7)=-. 23160355077925900 
XI 8) =-.09 50 l 2 50983763 700 
XI9I=. 095012  5''9337637  DO 
XI )0)=. 23160365077925900 
X(  HI*. 4530 1. o 77 765722  700 
X(  12)  = . 61  / 8 7 6 2 4 4 4 02  7>4  00 
X'  13)  = . 7554  0'-'.  33  3 55  00  300 
XI  14|«.j;,56  3120238783200 
X(  15)  = .9445  76  32  3.3  73  23  300 
XI  16)  =.  989-.C09349  91 65000 
C 
C 

C CALCULATION  OF  SAHA  E0UILRR1UM 
00  2 <=2.60 

2 SAL  I'M  K)  =4.2E-;t  *(ELFr|)(10)*»2)  = K**2)*(FXP(lK7390./((x*«2)* 
1TEMP( IT) ) ) J/1TEMP1 I T)»“l.5l 

C 
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C CALCULATION  OF  K1P.C) 

T=TEMP( IT) 

DO  4 J=*  1 1 30 
P = J 

$UM=0.0 
DO  3 K=  l • 16 
8=C(K) 

R = XI  K ) 

El  = l .0D5D*DEXP(-7m45 .00MP  + 3.U0)  / 1 P+*2.D0*T)  1 
£2  = 1. 0050+08  XP  I -7,31560.  000/  ( P**2*T  * ( l .0U0-P  ) ) I 
F 1=2. 0*1. 4 1 42D0< DS0R r ( 1 1 « * l . 03 1 / ( 3 + 5 . DO ) ) * * 3 . DO ) /3 . DO 
F2=  14.00*1  7 . D 0 + 3 • 0 0 *R  1/(3  .DO*  ( 1.00-R  I ♦ * 2*1)50*  7 ( ( 5 . DO-R ) **3)  ) 

3 SUM=SUM+8*<F1*E1+F2*t2> 

FINAL=5.45U0*SUM/U)SQRT (T**3» *1 .0350) 

4 COL  I 7 ( J 1 =F I MAL 


CALCULATE  K(P.O) 

DO  5 J=l • 30 
P = J 

JS=J»l 

DO  5 K=J$,40 
SUM  = 0 . 0 

A=1K**2)/(K**2-P**2) 

DO  6 L = l • 16 
B = C 1 L 1 

B = X ( U 1 

El  = 0ExP(-70945. 000*19  H.O  I /<P**2*'f  1-157390.000/1  P**2*A*T  1 ) *1.0  030 
E2=DEX PI-315780. 000*1 1 . 0 + A 1 / ( P**2* A*T* ( L. 0-R I) ) * l . 0030 
F 1 = 1 4 . 0*4  **2  +7.0*A+o.0+4.0*A=*2=R*A*R)*0S1HT  ( 1 R + 1 .01*14+1  .01  1*A/ 
41 6. O+DSORT  t I 3 • 0*A  + 2 • 0+  A *R 1**3) ) 

F2=l  8. *4**2 +1  l .*A+6.+3.*A*R)*DS09T  ( ( 2.00*1  1 . 00  + A ) ) **3.00) 

1/1 3.00*»S0RT<  < 3 .00* A + 2. DO- A*R 1**3.001*1 l.-R)**2) 

6 SUM  = 5UV*  + S*(  Fl*El+F2*E2) 

F 1NAL  = 10.9D0*P**2*A*SUM/1 OS  OR T ( T + * 3 1 *K**3 I *1.00-3  0 
5 COLE  X 1 J,K)=F1NAL 


CALCULATE  KtQ.Pl 

CO  8 J = 1 1 39 

K-41-J 

DO  7 L=1.30 

IF!  L.GE.'OOO  TO  3 

7 COLOEX IK,L)=SAEQP(L)*CCLEX!L,K) /SAEOP(K) 

8 CONTINUE 


CALCULATE  A1P) 
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DSUM*0.0 

JJ=J*1 

00  15  K=JJ,29 

15  USUH  = DSUM»COLEX(  J.K» 

14  RKAPAI  J1=R«APA(  JJj-OSOM 

on  16  J=2.30 

D$UM=0.0 

JJ=J-l 

00  17  K=1.JJ 

17  OSlJM  = OSUMt-COlOEX{  J.K) 

16  RKAPAI  J) =RKAPA(J 1 +0  SUM 


NOLEV =2 8 
N01-N0LEVH 
N02=NCLEV-l 
N03=M0LEV-2 

IF(N0LEV.LT.30)G0  TO  801 

PR  I NT  802 

STOP 

802  FORMAT f l* ,35X, • LEVELS  EXCEED  29*) 

801  NO 1 =NOL£ V* 1 

CALCULATE  X ( P i S IGMA ] 

00  18  J=l,NCLEV 

DSUM=0.0 

00  19  K = N01 > 29 

19  DSUM=DSUM«-COLEX<  J.K) 

18  SIGMMJ)  = OSUM 

CALCULATE  ASUM 
00  20  1=2,28 
0SU',  = 0.0 
00  502  J = NO 1,29 

502  OSUM=OSUMK  SAEQP( Jl *SPONA( J, M/SAEQPII  ) ) 

20  ASUM( I)=OSUM 


CALCULATE  GSC 
DO  24  1=1, NC2 
N1  = IH 

24  GSC  ( I , 1 > =ELECO<  ID  I*  < S I GW  N1  1 ♦ COL  IZINlllM  IELECOI  ID  ) «=»2  >* 
1EETAI  IT.Nl  l/SAEOPINl)  ) « ASUM(Nl) 

DO  27  1=1, N02 
27  GSC (1,2)  =GSC(I,l) 


CALCULATE  IR  LOSS  TERM 
00  1315  <MN|=  l , 2 

00  416  I L = L , 16 
DO  10  1=1.30 
GSIHU  1,21=0.0 

10  GS IP  L ( I, l 1=0.0 
W = P LAN3 ( IL, l 1 

1 P M = P.  L A M 3 i IL.21 
DO  9 J=  l P'1, 30 

9 GSIP.LU.21  =(5.252L31S*(^**4)/(J**51 


CALCULATE  GS  COEFFICIENTS 
709  DO  21  L»1.N02 
Nl*L»l 
K=  ID 


on  o 
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c 


= ei  rcrm >*»kapai ni ) 

-LLtCOfK  I’»RKAPA(N1J 


*AE(NI»*GSIPLIN1,2» 
* A£ I N 1 I 


=-ELCCD(K ) *C0L0EX(N2,N1) 
=-CLEC0(K.)*C0LUEX(N2  ,Nl) 


GSIL.L.2 J 

21  GSIL.L.l) 

00  22  L’l.NOi 
Nl  = L«-l 

00  22  M=N1,N02 
N2  =M* 1 
GS(M.L.L) 

22  G S ( M , L • 2 ) 

00  23  L=l,N03 
N l =L+ 1 

00  23  M=Nl,N02 

N2=M*1 

00  23  1=1,2 

i ? GS(L.M,n=-<ELEC0tlDJ*C0LEXINl,N2)  ) - t S AEQP  ( N2  ) ♦SPO'lA ( N2 ,N  1) 
l/SAEOPINl) ) 

GAUSS  SICOEL  METHOD  TO  SOLVE  EONS. 

DO  25  1=1,2 
I TER=l 
99  8 I G=0.0 

00  100  L= 1 , N02 

102  CSUM=0.0 
IFIL.EQ.l)  GO  TO  105 
LAST=L-l 

00  106  M= 1 , L A ST 

106  OSUM=  OSIJM  *G3  ( L , M , I ) + R HO  I M,  I ) 
tF(L.EO.N02)GO  TO  103 

105  INITL  = LH 

00  107  N= IN  I TL  »N02 

107  OSUM=OSUM»GS ( L , N» I ) *PHO(N,I) 

103  T T E Mp=  ( GSCIL.I  )-DSllM)/GS(L,L.I  ) 


IFIABSITTEMO-RHGIL, I H.GT ,3IG)QIG=A8S( TTEHP-RHOIL,  1)1 
100  RHOIL . I )=TTEMP 
EEPSI I )=8IG 
UTERI  IJ  = ITER 
IF! BIG.LT.EPS JGO  TO  25 
110  IF! ITER. GE. MAX) GO  TO  25 
ITER  = ITERU 
GO  TO  99 
25  CONTINUE 


CALCULATE  PERCENT  QUENCHING 

DO  26  L=  1 « NCJ2 

LL=L*1 

26  OUF.NILI  ) = (RHO(L. ll-RHOIL. 2)  l/RHOIL,  11*100. 

PRINT409, NOLEV 

409  FORMAT (• 1 • ,55X, 'RESULTS' , 10X, 12 ,2X, • LEVELS*  ) 

PRINTS  10, TEMPI  IT) .ELECOI 1 0 i , RLAM3I I L , 1 1 

410  FORMAT  I •-*  , • TEMPERA TURE  = ' ,F  /.0,25X, 'ELECTRON  DENSITY’ • , E12. 1.25X, 
l 'WAVELENGTH  OF  IR=',£10.3) 

PRINT 411 , 1 1 T P R I 1) , I ITERI2) 

411  FORMAT ( , ' I TS9AT1 ONS  NO  IR='»I3«6QX.' ITERATIONS  l KW  IP=',I3) 

PR  I MT4 12, EEP  S I l » .EEPS 12) 

412  FORMAT ( • ', 'EPSILON  NO  I R ' , E 10 . 2 , 56 X . • f PS  I LON- I R • , E 10. 2 1 
PRINT6U 

511  FORMAT!'-') 

00  tOl  1=1, NC2 

n = i*i 

A1  = SAE0P(  II  )*Rllpf  I,  l) 

A 2 = S A E 0 ' I II )*RHOI 1.2) 

601  PRINT602.I I » A 1 , I I ,A2 

|! 
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602  FORMAT!*  *, ‘DENSITY  (*,I2,*J=  • ,E16.7,60X, *IR  DENSITY  (•, 

1 1 2 • • 1 - '.61 6.M 
00  <.11  l=2,NULEV 
11=1-1 

613  PRINT616.I.PH0! 1 1 . I > , I . PHO ( II , 2 ) 

616  FORMAT ( • •, •'UH! ,,!2,' )=  • ,616. 7, SOX, 'PHO  13  t*,I2,,)=  *,C16.7) 
PR INT610 

618  F 0 0 ’ A T I • 1 ' ) 

PR  I NT  6 15 

615  FORMAT ('-• , 50X, • PERCENT  OUENCHI NG*  J 
03  616  1=2, NOLEV 

616  P R I N T <*  l r , I , OIJCNI  I ) 

617  FORNATI *0* .65*. 'QUENCHING  IN  LEVEL  • , I 2 ,6X , F6.2) 

03  1313  KR=  1 . 16 

1313  RLAMB(K3,2)=RLAM0(K0,2)«-1. 

1315  RUTT=2. 

STOP 

END 


non 
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C CALCULATION  OF  EXCITED  STATE  POPULATIONS  AND  PERCENT  9'JENCHING  FOR  HELIU« 

C CONSTANT  INTENSITY-  1 K'<;  / CM**2 
C VARIABLE  NUMBER  OF  LEVELS 

C NOLEV  IS  THE  NUK3EP  OF  LEVELS  USED  MAXIMUM  LEVEL=29 

IMPLICIT  RE  AL*3  ( A-H.O-/. ) 

01  MENS  I ON  SACO?  05 .4 , 7)  ,COL  1 7 ( 3:1, 4 ) , CGL  EX  ( 30 , 3 5,  4 ) , OS  I»L  (.10 1 2 I , 

1AE( 30) , SPONA ( 35 .35)  ,PETA<  2 ),4) .RKSPA ( 3 Y) , SIGM\( 39) , A SUM! 30)  , 

20 S(  29,29.2)  .GSC(29.2)  ,RHi,I29,2)  .OJENI  3-3),  I ITEA<  2)  ,EE°S{  2)  . 

3CLECDI 7)  .TEMP (A) .RLAMIM  IT. 2) . EHEI35.4 ) , COLOEX ( 35 , 30 ) 
data  EPS/UD-S/.MAX/35  0/,  I ITE"/2'‘D  / .ELcCD/l.l)a,5.DB.  I .09,5.09, 

11  .010,5.010.  I.  Dll/,  . ( HV'BI  I . I ) , 1 = 1 , 17)/ 

21.460-6, 2. 28 0-6,3.230-4, 4.460-6, 5. 82 0-6, 7.330-6, 9. 110-6, l C. 60-6, 
311.0-6, 13. no- 6, 15. 390- 6, 17. 350-6. 20. 4 9 0-6, 23. 320-6, 26. 32 0-6, 
429.510-6,32.330-6/, (RLAM3I I ,2), I =1.17 1/4.09, 5.09,6. 03 ,7. 00,8. DO, 

59.00, 10.00,  11.  DO.  11.00.  12.  DO,  13.00,  l <-.1)0,  15.00,  16.00,  17.00, 13.00, 
519.00/. 

6C0LUEX/ 1050*0. 00/, R HO /53» 0.00/,  TEMP/ 300.0 0,400. 00, 500.00, 1000.00/ 

00  290  1=1.4 

200  REAO.IBETAl J.I) ,J=l,29) 

00  201  K=l,7 

00  201  J=l,4 

201  READ,  ( S4E3P ( I . J,K) , 1 = 1.35) 

DO  202  1=1.4 

202  REAO,1COLU(J,I),J  = 1,30) 

REAO.COLEX 

REAO.EHE 
RE AO, SPONA 

READ, (SAE0PI2. 1,1), 1=1. 7» 


CALCULATE  A(P) 

00  12  1=2,30 
11=1-1 
0SUM=1 .04 
DO  II  J=l.II 

11  0SU«=0S0M»SP0NA(1,J1 

12  AEtn-OSLlM 
00  1317  l T = 4 » 4 
00  1317  10=7. 7 

C CALCULATE  C (’  L E'  E X ( 3 5 , 3 0 • 4 , 7 ) , < ( P , 9 ) P.GT.O 

00  l<-‘2  1=1,35 
03  1442  J=1 , 30 
1442  COLOEXI I, J)*0. 

KK  = I T 
JJ’ll) 

DO  8 J=l , 34 
K=36-J 
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DO  7 I.  = l . 30 
IFIL.GE.KIGO  TO  8 

7 COLOEXIK.L)  = $AtOP(L.KK,J.))*:OLEXIL,K,KK)/SAeQP(K,KK,JJI 

8 CONTINUE 

I F ( IT. ME. 1)00  TO  1503 
00  1501  1=3,35 

1501  COLOEXI r.2)=C0L0EXI 1.21*1.020 
1500  BUT  = 1 . 

C 

C 

C CALCULATE  KIP) 

DO  13  1=1,30 

13  RKAP4(I)--C0LtZ(I,IT> 

00  14  J=l,23 
0SUH*0.0 

JJ  = .IH 

00  15  K = J J » 35 

15  OSUM  =0SUNi-C0LEX(  J.K  , I T) 

14  RKAPAI J)=RKAPAIJ)+DSUM 
00  16  J = 2 » 30 
OSUNrO.O 

JJ=J-1 

DO  17  K * l . J J 

17  USUM  = OSUM*-COLOEX(J,K> 

16  RKAPAI J)=RKAPA( J) ♦OSU* 

C 

c 

N0LEV=23 
NOl  = NC'LEVH 
N02=N0LEV-1 
N03=N0LEY-2 

IFIMOLEV.LT. 30)00  TO  801 

PR  I NT  802 

STOP 

C02  FORMAT  I • 1* ,35X, • LEVELS  EXCEEO  20‘) 

801  N01=NfJLEVU 
C CALCULATE  KIP, SIGMA) 

DO  18  J=l, NOLEV 
D SUM =0.0 
DO  19  K=NOl , 35 

19  OSUNsDSU’^COLEXI J.K.IT) 

18  SIGNKI J)  =OSUM 

C CALCULATE  ASUM 

00  20  1=2.23 
0SUM=0 . 0 
DO  r02  J = NO  t , 3 5 

502  OSJMmOSUMK S4E0PI J. IT , ID) ♦SPQNA I J , I I /SAEOPI I . IT. 10) ) 

20  ASUMI  U=OSUM 

IF!  IT.EO.  1)  ASUMI  2) ---ASUMI  2)*  1.0-20 

c 

c 

C CALCULATE  GSC 

K=  10 

DO  24  1=1, N02 

n i = r «•  1 

24  GSC  I 1 , 1 ) = ( EL  T CD  I ID  ) * I S 1 GWX  INI)  «-C  1L  I l ! N 1 . I T 1 I «■  ( ( TLECDI  I D ) **2  ) * 
lDETACIl,  IT  )/ SALOP  I Ml  , IT  , 10)  i * AvJ’II  M ) ) * L.U45 
00  27  1=1, N02 

IT  I IT  .EO.  I > “Sr  I 1 . I > = 1 ‘LECOI  ! D)  MS!  O^KI?  ) 'COLT  M 2.  l))»l  .045)  ♦ 

II  I EL  ECU  I 10)  *»2  ) *'UT  A<  z.  1 ) * 1 .026/S  AE  <)P(  2 , l,  ID)  ) ♦ I 4SU*'I  2 I *1.045  I 
2 7 GSC  I 1 ,2)  =r,SC(  1,1) 


62 


C 

C 

C CALCULATE  GS  COEFFICIENTS 
00  22  L=1,N03 
Nl=L*l 

00  2 2 M-’-fl.NOZ 
N2  = M*-l 

GSIM.L.l  ) =-ELFC0IK)®C0LDEX(N2,Nl)* 1.045 

22  0SIM.L.2)  = -ELEC0(K)*CGL0t:X<N2.Nl)*l.045 

on  23  L = 1 » N03 
N l = L *1 

00  23  M=Nl,N02 
N2  = MU 

DO  23  1=1.2 

GSIL.M,  I ) =-lELF:CD{  t D)  *COLEX ( Nl , N2 , IT)  )-<S4E0P(N2,  IT,  1 0 ) *SP0N4 { N2 ,N 

11 ) /SAE0PIN1 . I T. 10) ) 

23  GSIL.M,  I )=GS(L.M, I ) *1.045 

1 F ( IT.NE.DGO  TO  1502 
00  1503  1=2, N02 
11=1*1 

1503  GStl»I»l)=GSll.I«2)=-(FLECD(IO)*C0LEXI2.II»l)*l.O45)-((SAEQPlII« 

1 l . 10)*SP0NA( II ,2>/SAE9P<2, l. 10)1*1.025) 

1502  BUT =2. 

C GS  E GSC  OIVIOED  BY  10**45 

00  63  Kl=l,2 
00  63  K2=? , N02 

68  GSC(K2,Kl)=GSC(K2,M)*1.0-45 

00  69  Kl  = l,2 
DO  69  K2  = 2 . NQ2 
00  69  K3=2.N02 
IFIK3.EO.K2 JGO  TO  69 
GSIK3,K2,Kl)=GS(K3,K2,KD*  1.0-45 

69  CONTINUE 
C 

C 

C CALCULATE  IP  LOSS  TERM 

DO  1315  KMN=l,2 
DO  416  IL=l.l7 
DO  10  1=1.30 
GS  I °L ( 1 ,21=0.0 
10  GS I P L ( 1,11=0.0 
W=RLAM8( IL.l) 

IP«*RLAMRl 1 L ,2 ) *8. 

00  9 J=!PM,30 
JJ=J-3 

9 CSIRLtJ.21  =( 5.252031 ) * I W* *4 ) /( J J>*  5 1 
709  00  21  L=1.N02 
Nl=L*l 

GS(L.L.2)=  (ELECOIK )*RkAPAIN1>  *AE( N l ) *GS I RL ( N1 , 2 ) ) 

21  GS( L , L , 1 ) = (FLECO(K)*KK APA(Nl)  +AEINI)  ) 

C GAUSS  S I FUEL  METHOD  TC  SOLVE  EONS. 

00  25  1=1.2 

IF( IL.NE.I.AND.I.FO.IIGO  TO  25 
DO  1504  J 1 = l , N02 

1504  R HO ( J 1 , 1 ) =0. 

ITER* 1 

99  0 1 G=0. 0 

00  100  L =2 . N02 
102  OSUM^O.D 

I F ( L.E0.2 ) GO  TO  105 
LAST  = L-1 
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00  106  M=2,LAST 

106  l)SUM  = USUM  + GS(L,M,  I ) *SH0!M.  I) 

IF(L.EQ.N02)G0  TO  103 
105  I NI TL  = L +1 

00  107  N=tNTTI..N0? 

10  7 DSIJU=0St)M+G3  ( I. , N . I ) + RHQ(  N , I I 
i ' 103  TTfcMP=(  GSCIL.I  ) -USOV  ) /GS I L , L , I) 

IF!OAUS(T  reNP-r.HU!l,  I ) ) .GT. ‘UG  HH  3 = DABS!  TTEMP-RHOIUI  > ) 

100  RHOIl, I ) = T TEMP 
EE  DS ( l ) = 3 1 G 
UTER(  1 )=ITER 
IFIB1G.LT. EPS1C0  TO  25 
110  IF!  ITEM  .GE.MAX1G0  TO  25 

1 TEP*I  TERU 
GO  TO  90 

25  CO NT  I MU 6 

DEN1  = 0ET12  =0.00 
00  6 7 K l - 2 , N02 

DENl  = 0fcMUGS(  l.Kl , U *RH0!K1,  1) 

67  DEN2=DEN2+GS I 1 • K l • 2 ) * RHO ( K l , 2 ) 

A A 1 - 1 • DO 

1 F ( IT.NE.  l)AAl  = 1.0-20 

DEN  1* I GSC ( 1 , ll-OtUl )*l.D-65*AAl*S4E0P<  2. 1 T , I 0 ) *1 . 020/GS ( 1,1,1) 
0EM2=(GSCl  1.21-0EN2  >* l.D-65 «AAl*SAEQP ! 2 . 1 T , I 0 )*1 .020/GS ( 1 , l , 2 ) 


CALCULATE  PERCENT  O'JENCHING 
QUEN ! 2 ) = ( DE  N l— 0EN2 ) * 1 .02/ OEN 1 
00  26  L = 2«  NC2 

LL=L+1 

26  OUEN(LL)  = (RH01L. I 1*  l.D6-RH0!L,2 )*l  .06)/ (RHOIL , l 1*1.061 
PRINT  6 09, NUL  E 7 

609  FORMAT  ( ' l • , 5 5X  , ' R E S‘JL TS  * . 10X  , 1 2 , 2 X , • LEVEL  S ' J 
PR I NT6 10. TEMPI  IT) , ELECD! I 0) , RLAMB!  IL .1 ) 

610  FORMAT (•-«, * TEMPERATURE** ,F 7.0,25X,  'ELECTRON  DENS  I TV= * , E 1 2 . 1 , 2 5X, 
l* WAVELENGTH  OF  IR=',E10.3) 

PRINT  411,  I ITER ( l ) . I ITEP (2) 

611  FORMAT! ITERATIONS  MO  I R = ' , I 3 , 60X, • I TERAT I ONS  1 KW  I a = • ,131 
PR  I NT6  12 • EE°S ( l ) . EE°S(2 ) 

612  FORMAT!'  ', 'EPSILON  NO  IR» , E 10. 2 , 56X , *E PS  I LON-I R • , E 10.2 » 

PRINT511 

511  FORMAT!*-*) 

00  COl  1=1, N02 
11=1+1 

A 1 = SAEOP!  II,  IT,  ID)  *P.HO!  1,1) 

A2  = SAEQP( II • IT, 10)mRH0( 1,2) 

IF! I .EQ. 1 ) A1 =DEN 1 
IF!  I .60. I ) A2=OEN2 

601  PRIMT602.II  ,A1,II,A2 

602  FORMAT! • • , 'DENSITY  , • . 1 2 , • ) = • , El  6 . 7 , 6 OX , • I R DENSITY 

112.'  )=  ' ,E16. 7) 

DO  613  I =2 , MOLE V 

1 1 = 1-1 

613  PR  J N T 6 1 6 , 1 # 3 HO I 1 1 # 1 ) , I , RHP! 11,2) 

616  FOR"  AT  ( * ', 'RHO!  ' , I 2,  ' ) = ' , El  o.  7 , 5 OX  , • PHO  IR  !».I2»'>-  *,F:16.7) 
PR1NT61 0 
618  FORMA T ( * l • ) 

PRINT615 

615  Ftl°"AT!*-'  , B 0 X • * P C R C C NT  0' IF  VC  "I  MG'  ) 

A VE  = ( OMEN  ( 1 1 I ♦QUEM!  1 2 I +OUEN!  I 3 1 +OJLMI  16)  1/6. 

PRINT, AVE 
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DO  416  1^2. NOLEV 

416  PRI.NT41  7,  ( ,OUEN<  I 1 

417  FORM A r ( • 0 1 » 45X i ' QUENCH I NG  IN  LEVEL 
DO  1313  K0=l,l7 

1313  ft  L A '-1  B ( K B « 2 ) =PL  AMrt  ( KB,  2 ) H . 

1315  BUTT=2. 

DO  1316  KB  = 1 1 1 7 

1316  RLAM!3(  K0 , 2 ) = R L AM B < KB , 2 I -2 . 

1317  CONTINUE 
STOP 
END 


C RADIATIVE  RECOMBINATION  OF  HYDRO OEM  I C IONS 
DIMENSION  T A *5 1 A , 9 7 ) .ELI  AO  1 1 BETA (AO) 

READ, < TAB! 1 . I i ,TA8<  2,1)  ,TA3<  3, f ) .TAGIR, I)  , l -1  ,97) 

PRINT  9,  ( I ,TAR,<  l , I ) .TAB  (2, 1 ) ,TA:i<3,  I ) ,TAB<  A,  I ) , 1=1,97) 
9 FORMAT  I • • ,5X,I2,2X,F10.5,2<»F10.5,2X,F10.5,2X,F10.5) 

; REAO.T.N 

READ, ( ELI  1 1 ,1  = 1, N) 

0-5.I97E-14 
C=I57890./T 
C3=C*»I-l./3. ) 

CC=SQRT(CI 
DO  l I = 1 . N 
X=C/CL 1 1 ) **2 
Z = 3 . * X 

IFIX.LT..02IG0  TO  10 
IFI X.GT.20. 100  TO  20 
IFI  X.IE..OMGO  TO  31 
1F(X.LE..I)G0  TO  32 
IF(X.LE..2)G0  TO  33 
IFI X.LE..41GQ  TO  34 
IFIX.LE.l. }G0  TO  35 
IFIX.LE.2. )G0  TO  36 
IFIX.LE.4. JGO  TO  37 
ir-(  X.LE.IO.  1G0  TO  38 


IF( X.IT.20. ) GO  TO  39 
XSO  = . 95AA 
XSI=-.AIA7 
• XS2=-  . 350 
GO  TO  2 

10  XSO=X*t  XP<  X)  *(  X-ALOGI  X 1-0.5  T 12) 

XSl=0.t629*  X*{  l.  «-A.*X)-1.0368*X**l  . 3 33333*  ( l.  .875*X) 
X$2=-0.0672«X*(  l.«-3.  = X)«-0.  l',30*X**  l.  664  65  7*  ( l . U.8*X) 

GO  TO  2 

20  XS9=l.-l./X*-2./|X*X)-6./(X*X*X)*2A./(X*X*X*X) 

XS I’-O. I 728*X**0. 3333333*1  l .-8. /2 ♦ 70. / 1 7* 7 )- 800. / 1 X*X*X > nmO./IZ 

*+i*i*n  i 

XS2=-O.OA96*X**0.6666667*( l.-3./Zv32./<  7* 7 » -AA8 ./ < 7*7*7 J ) 

GO  TO  2 

31  01=500. 

J=l.*0»*<X-.02) 

GO  TO  3 

32  01=200. 

J=12.*0I*<X-.0A5) 

GO  TO  3 

33  01=100. 

J=2A.+Dt*(X-,ll) 

GO  TO  3 

3A  01=50. 

J = 33.+DI‘MX-.20) 

GO  TO  3 

35  01=20. 

J=4A.*9I*1X-.A5) 

GO  TO  3 

36  01=10. 

J=56.*DI*(X-1.1) 

GO  TO  3 

37  01=5. 

J = 65.*DIMX-2.9J 
GO  TO  3 

38  01=2. 

J=  76.  ♦ DIMX-A.5) 

GO  TO  3 

39  OT  = 1 . 

J=8U.*0I*tX-ll.) 

3 P=0I*1 X-TAB1 1 . J) 1 

XS0=TAB(2,.||*P*(TA3(2»J*1)“TA(?(2«JJ1 
XS t=TA3{ 3, J I »P*1TA8( 3. J»1 » - T AP ( 3, J)  1 
XS2=rAri(A,J)+P*(TA3(',,JHi-rA6(A,J)) 

2 8tTA<! >=0*CC«( XSO+C  1*X$ 1 »C 3*C3« XS2 ) /CL ( I) 
l CONTINUE 

PRINT  At ( I « 3ETA I ! > , I =1 ,N) 

A FORMAT  (•  « , 10X, ‘BETA! * ,12, ‘ ) = ’.E17. 7) 
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CALCULATION  OP  R ADI 4T I VE  RECOMBINATION  COEFF-  HEL.UM 
BE TA(  291 

READ  IN  I3ETA(291  ,DIE(  35,'») 

01  "ENSIGN  8ETA{29),EUE<33,4)  , T(4) 

READ, T 
READ.EHE 

00  l 1*1,* 

RE  AO, BETA 
00  2 J = 1 , 2 9 

2 BETA! J ) =BETA(J )P£HC IJ ,4)/(4*EHEI J, 3 )*EHE<  J,3) J 
PRIMTA  ,T( I ) 

* FORMAT  <*  1 •,  10X,  «Te'MP.=  * , F 7. 1 1 
00  5 L L = 1 , 2 9 

5 PR INT6  » LL • BETA ( LL  ) 

6 FORMAT ( • ■ , 10X,  ' BETA! • , 12, ' ) = *,616.7) 

PUNCH, BETA 

DO  3 K=L,29 

3 f)ETA(K)=0.0 
l CONTINUE 

PUNCH, EH E 
STOP 
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CALCULATION  OF  CCLISIONAL  EXCITATION  ANO  IONIZATION  HELIUM 
COL IZI 30) ,C3LEXI30,35),T14J,£HEI35,4> 

DIMENSION  C0LIZI30J.C0LCX(30,35),  EHE135.4) 

DOUBLE  PRECISION  P .SUM , 8, R , El , E2 , F 1 , F2 , C 11 6 ) , X ( 16 J , DSQPT, DE XP , 0, 

IT  < 4 I 
READ,! 

REAO.EHC- 

EHEf 9.3»  = S0PT(  1 09673. /1 22 10.) 
cHEI 10,3) =SORT< 1096 78./ 12206.) 

Cl l >=.02715245941175400 
Cl  2 > = .062253  52  39  1864800 
Cl  3>  = .C9515  85U682493D0 
C 1 4 > = . 1246239712555  3400 
C! 5 >*.1495959088165 7700 
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C 16 >=.1691565193950  )JDO 
Cl /)«.  18260341504492400 
Cl  31  = .18945  )61  0455  06900 
Ci9)*.in94506io45506900 
C(  10)  = .lrt240<4:*O44424D0 
Cl  ll>  = . 16915651939500300 
C(  12 >«.  14959598:581  6 5/700 

Cl  m *.1246239/125553400 
Cl  14>=. 09515651 168249300 

C(  15>*. 04225352393364800 
Cl  161*. 02715245 9411  75400 
XU  1 =-.93940  093499165000 
XI 21 =-.9445 7502307323300 
XI  J)*-. 36563120238703200 
XI 4)=-. 75540440035500300 
X ( 5 ) =- . 61 737624440264400 
XI t 1 =- .45 301 67  7/65/22  700 
XI 7)=-. 20160355977925900 
XI 0) =-.09501250983763700 
XI  9)  = .095012509 *5376  3700 
XI 10) =.2dl6035507 7925900 
XI  111*. 45301677765722700 
X l l 2 1 =. 6 1 78 7624440264400 
XI  13)  = . 75640440335500300 
XI  14!  = . 66563120233783200 
XI 15 ) = • 9445  7S0230  73  233D0 
XI  16)  = . 98940093499165000 
C CALCULATE  K(P,C> 

DO  4 L*i «4 
00  3 LL=l*30 
3 COL  I 21 LL ) =0.0 
00  l J=2«  30 
P*EHEI J,3) 

SUM*0.0 
DO  2 K=1.16 
B*CtK) 

R=X(K) 


61  = 1 . 0050*08 XP 1-7 8945. DO*<R  + 3.DO)/(P**2.00*ITIL)))) 

62* 1 .0030*0EXP l -43 1 569.00  / 1 P**2*T  l L ) *1 1.00-3  > > ) 

61  = 2.*  1.41 4209*0SOPT< U R* l . DO ) / U + 5 . 00) ) ** 3. 00 ) /3.D0 
F2--=  16.00* ( 7. 00+ 3.00*31/ l 3.00*1 l.OO-R > **2*0$0RT  t ( 5 .OO-R)  *=*3> ) 
2 SUK=SUM*8*lrl*£l+F2*E2> 

FINAL=5.4500*5UM/(0SQ8T(T (L ) **3 ) *1.0501 
1 COL  I 2 « J)=FINAL 
PRIMTtO.TlL) 

10  FORM  ATI* lUlOX,* TEMP.*  *,67.1) 

00  11  NRT*1 ,30 

11  4R1MT12.MRT.C0L12INRT) 

12  FORMAT { • ••10X,'CCUZl*tI2*')a  SEia.7) 

4 PUNCH. COLU 

CALCULATE  KIP.O) 

00  7 N= l « 4 
DO  8 NN=l,30 
00  3 NM»  1 • 35 
8 COLE XI MN.NM) =0. 

DO  5 J-»2, 30 
P = fcHH  J.3) 

JS  = JU 

00  5 x=.)5,35 
0*EHE ( <»  3 1 
SU1*0.0 
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A*  10**21/ (0**2”P**2 ) 

00  6 1=1.16 

B=C(U 

rt=xiu 

EI=0EXP(-7R9A«;.00*ta*l.  )/ » P **2*T  ( 9 1 1-15  73<»0.(>0  / ( P**2*  A*T(  N III  * 
11.0030 

E2*DEXP(-3157ilO.DO*Il.*AJ/(P**2*A*T(N)*l  l.-RI)  1*1.030 
F|r  U.*A**?*7.*A«-6.*a.*A*«2*R*A*R)*US0RTI  (»*■!.  )*(AU.  I)  *A/ 
AI6.*0S0RT( (3.*A*2.*AftR)**3l I 
F2«(8.*A**2* ll.*A+6.»3.*A*R)*0S0RT(  (2.00*1 l . 0')»A  ) 1 **3.001 
l/<  3.0T*0SC3T( (3.00*A*2.00-A*B )*«3.D0l*( l.-RI**2» 

6 SUM=>SUM.H*(F1*E1*F2*E2I 

FINAL* I0.900*P**2*A*SUu/< DSQRT(  TIN )**3)*K**3| 

IFIFINAL.LT. l.t -A 0)FtMAL=0. 

F IN»L*FINAL* l .E-30 
5 COLEXI J.K)=FINAL 
PRIMTIO.T(N) 

DO  13  NOC=l .30 
NOC  l =NOC* l 
DO  13  MAC *NOC 1.35 
13  PR l NT1 A, NOC. MAC. CCL EX (NOC, MAC  I 

1*  FORMAT ( • '.10X  .'COLEXI • ,12,' ,12. • )»  '.E16.7I 

7 PUNCH, CGLEX 
STOP 

END 


C SAHA  EQUILIBRIUM  DENS I TIES- HE HUM 

C SAEOP (2  I SET  TO  0.  S INC E> l . E 75-CORR ECTEO  IN  M.P. 

C SAEOPI35I-P6AI)  IN  EHE 

DIMENSION  SAEOP  I 35)  . EHE( 35.4 J . T ! 4)  . ELECOI  7) 

DATA  T/300.  ,A00.,500.,  1000. / , ELECD/ 1 .E8 , 5 . E8. I.E9,5.E9,UE10, 
15.EI0,  l.El 1/.F/4. 2E-16/ 

READ. EHE 
AND*9. 99999 
00  l K*l. 7 
00  l J = l.4 
00  8 KK*1.35 
8 SAEOP (KK)=0« 

00  2 1*2.35 

IFU.E0.1.AND.I.E0.2I  00  TO  2 

SAEOP!  I)*F*(ELECD(K  )*-»2)*(EXP(EHEU«2)/I.695*TI  JDII* 

II EHE! I ,41/4.)/<SQRT!T(J)**3) ) 

2 CONTINUE 
PRINT3.T! J) .ELECOtK) 

3 FORMAT ( ' 1 ' » 10X » * TEMP=  » ,F 7. I , 1 OX, • ELEC.  OEN.  * *,E12  • A) 

00  4 LL*1 .35 
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C CAlCUt  ATP  SPf'N  TANEOUS  T R /.MS  I T | ON  PPtJP.  - MILIUM 
DIMENSION  =nr(  39, 4)  , SP-»NA  ( v; , 35)  ,C  ( 1 1 ) 

VFAH, ICMcI I , 1 ) • ( = 1,  1 1 ) 

PFAO  , ( FME  ( I ,3)  ,1=1, ID 
RFAO,<EHEI  1, A), 1 = 1,11) 

on  i i = L,n 

1 EHFf  I ,2)=l<J3-»05.-£HEC  1,1) 

00  2 1=12,15 

J=< I-8)**2 

5HF1  ! ,1)*19®M5.-I  109475./ J) 

FHF  ( !,2)*19C4-ra./J 
SHC(I ,31=1-8 

2 EMC  ( I ,4)=4*  J 
00  3 J=  l,  35 
00  3 1=1,35 

3 SP0NA1 I , J 1=0.0 

R FAD* ( SPPNA (1,1), 1=1, 20) 

CO  4 J=2  ,5 

4 P.  E A 0 » (SPCNAll, J)  ,1  = 1,  IB) 

00  5 .1=6,11 

5 READ, CSP0NA(!,J1, 1=1,141 
00  6 J=12,27 

K = J*1 

6 REAO,(SPCNA(L, J) ,L=K, 281 
RFAO,C 

00  7 1=21,35 
J=  I-B 

7 SPCNA ( I , 1 ) = C (1 ) * EHE 1 1 *4) *1 EHtl I , l ) *=2 )/ l J+*5) 

00  8 K = 2 » 5 

CO  8 1=19,35 
J=l-8 

8 SPCNAl  I , K)  = C (K  ) *EHE  1 K ,41*11  EHEI  1 , 1 ) -EHE  (K , 1 ) ) **2  ) / l J*«=5  ) 
00  9 K=6,ll 

00  9 1=15,35 
J*I-B 

9 SP0NA1I ,KJ  = C1!<)*EHF<K,4)*1  1FHE1 I, l)-EME(K,l)  )**2)/I J**5) 
CO  10  1 = 29,35 

J*  1-8 

CC  10  *.=  12 » 29 

IF( I.EO.K)  00  TO  10 

L«K-8 

SPONAl I ,Kl=l.29*iL**3)*J*t C EHFC I ,1 )-EHE IK, l ) )**2)/ 

11 1 IJ**2)-ll**2  ) »**3» 

10  CONTINUE 

DO  11  1=1,34 
!!■ !♦! 

00  11  J= 1 1 , 35 

11  PR  1 N T 1 2 , J , I ,SP3NA1J,I  ) 

12  FORMAT ( * • , 14X, • A( • , ! 2,* , • , 12, • ) * • , E 16. 7) 

PRJNT15 

CO  13  1=1,35 

13  PRINT14, ( F IF  J I , J ) , J= 1 ,4 ) 

14  FOR  M A T 1 • • ,4C16.2) 

15  F0OM AT ( • 1 • J 
PP INT , EH6 
PUNCH, EhE  • 

PUNCH, SPONA 
STOP 

END 
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Picosecond  Dynamics  of  Photoelectric  Emission 

Detailed  study  of  the  time  evolution  of  the  photoelectric  emission 
has  never  been  investigated  experimentally  with  adequate  time  resolution. 
Without  even  taking  into  account  the  difficult  and  long  standing  problem 
of  the  time  dependence  of  the  probability  of  photoelectric  emission  on 
time  scales  on  which  the  oscillations  of  the  optical  field  are  resolved, 
as  considered  for  example  by  Handel  and  Meltzer , * one  can  state  that  the 
response  time  of  a finite  thickness  photoactive  surface  can  be  signifi- 
cantly affected  by  lattice  scattering,  carrier  lifetimes,  pair  production, 

and  formation  and  destruction  of  exclton  and  impurity  states.  Therefore, 

.15  -9 

the  evolution  of  emission  can  conceivably  range  from  10  to  10  sec. 

In  the  literature,  estimates  for  the  response  time  of  pbotocathodes  include 

-12 

such  values  as  "instantaneous,''  *10  sec  for  conventional  materials  in 
which  the  threshold  for  pair  production  is  large  compared  with  the  electron 
affinity, ^ 10  ^ to  3 x 10  ^ sec  for  similar  materials  in  different  circum- 
stances,^ *10  ^ sec  for  some  negative-electron-affinity  materials,^  and, 
surprisingly,  entirely  outside  of  the  microwave  modulation  region  for  the 


specific  case  of  cesium-antimony  deposited  on  beryllium-oxide.  Finally, 

6 

a result  of  a different  nature  but  not  unrelated  was  reported  by  Gabriel 
who  showed  that,  in  photon  counting,  given  the  emission  of  a first  photo- 


electron at  t ■ o,  the  probability  of  emission  of  a second  photoelectron 


Is  significantly  enhanced  at  t 9 15  x 10  sec;  the  work  was  carrle^  out 
-9 

with  *10  sec  resolution. 


It  is  true  that  studying  photoemission  on  picosecond  time  scales  is 


extremely  difficult  and  requires  the  most  modern  high-speed  instruments 
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UNIVERSITY  OF  ILLINOIS 

DISCLOSURE  OF  INVENTION  AND 
LETTER  OF  TRANSMITTAL 


To  the  University  of  Illinois  Date:  February  16,  1973 

Patent  Committee: 

Entered  herein  and  attached  hereto  is  Information  concerning  a 
potentially  patentable  invention,  for  your  consideration. 


1 . Description  Title  of  Invention 

2.  Inventor 

Name:  Henry  Merkelo  Office  Address  and  Telephone  Number: 

155  Electrical  Engineering  Bldg.,  333-2482 

Rank  and  Department:  Permanent  Address: 

Associate  Professor  of  Electrical  Engineering 

3.  Name  and  Address  of  Co-inventor  (if  any) 

4.  Description  of  Invention  (In  as  simple  and  brief  language  as  possible 
describe  invention,  clearly  pointing  out  what  parts  thereof  are  novel 

and  which  are  old  features,  and  which  features  were  invented  or  suggested 
by  colleagues.  Attach  two  sets  of  drawings , sketches,  photos  and  reports 
of  invention,  and  any  parts  which  will  be  used  in  the  description, 
referring  to  by  number,  if  necessary.  Attach  additional  sheets,  if 
needed  ) 

(a)  Description: 

High  speed  sampling  of  photoelectrons  is  demonstrated  and  applied  to 
the  construction  of  a high  resolution  "sampling  photomultiplier". 

The  outstanding  feature  of  devices  built  on  the  sampling  principle 
lies  in  the  attainment  of  extremely  high  time  resolution  of  periodic 
signals.  Opto-electronic  sampling  as  developed  here  is  analogous  to 
that  used  in  conventional  sampling  electronic  instrumentation  (oscillo- 
scopes, counters,  etc.);  however,  the  sampling  of  photoelectrons  is 
distinct  and  technologically  different  from  conventional  electronic 
sampling.  Signal  producing  (sampled)  electrons  are  physically  separated 
with  high  time  resolution  in  a specially  designed,  electrically  broadband, 
photoemissive  structure. 

A device  has  been  built  and  tested  on  light  pulses  of  ~300  psec  rise 
(10  to  9QZ)  fully  resolving  this  time  variation;  the  actual  speed  limit 
of  the  device  is  not  known  yet.  Work  continues  which  is  aimed  at 
establishing  experimentally  the  limiting  resolution  which,  on  the  basis 
of  quantum  mechanical  arguments,  should  lie  in  the  10  - 10“^  sec  range. 

The  remarkable  feature  of  this  high  resolution  sampling  photo- 
multiplier is  that  high  time  resolution  of  periodic  signals  can  be  obtained 
with  slow  auxiliary  instrumentation  (oscilloscopes,  counters,  etc.) 

*See  attached  technical  description  design  drawings  and  performance  tests. 
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Earliest  date  and  place  Invention  was  conceived.  (Brief  outline 
of  circumstances. ) 

Summer  1972;  proposed  M.S.  thesis  to  Sid  Kaiser;  Department  of  Electrical 
Engineering  (Undergraduate  Laboratory  of  Quantum  Electronics) 

Date  and  place  of  first  sketch,  drawing  or  photo. 


January  1973,  Electi'ical  Engineering  Publications  Office,  M.S. 
S.  Kaiser  in  preparation. 

Date  and  place  of  first  written  description,  if  available. 


Thesis 


January  1973,  Electrical  Engineering  Publications  Office,  M.S.  Thesis, 
Piscl£fuE#i3er>  in  Preparation. 

(a)  Disclosure  of  Invention  to  Others 


Name,  Title 
and  Address 


Jamie  Wiczer,  Res. 
Asst.,  Univ.  of  ill. 


Form  of  Disclosure 
(including  theses 
and  other  publications) 


. Verbal 


E.  H.  Eberhardt*, 
Physicist 


Verbal  during  visit  to 
University  of  Illinois. 


Date  and 
Place  of 
Disclosure 


Was  Signature 
Obtained 
(Yes  or  No) 


No;  Mr.  Wiczer 
new  M.S.  thesis 


■ »#:  • • i • ■ r.  ■ ■ 


No;  only  concepts 


Jan.  19,  1973, 
Seminar 


remainder  of  Disclosure  of  Invention  to  Others 


(b)  Date  and  Place  of  Completion  of  first  operating  model  or  full  size 
device. 

January  11,  1973;  267  Electrical  Engineering  Building. 

(c)  Present  Location  of  Model. 

267  Electrical  Engineering  Building. 

^ 6esultI:aC5An<?eli^i5^St:  £lecltl:ngrgf.il^cigt.e,s Wl&r^t¥on:on'S.  Kaiser 

M S Thesis  (in  preparation) 

Support  6f  the  University  of  Illinois  (Complete  explanation,  with  dates, 
of  University  of  Illinois  facilities  and  support  utilized,  in  whole  or 
in  part,  in  development  of  this  invention. ) 

Department  of  Electrical  Engineering 
Undergraduate  Quantum  Electronics  Laboratory 
Special  Processes  Laboratory 
Gaseous  Electronics  Laboratory 

Waiver  of  indirect  costs  on  NQ0014-67-A-0305-0017 


*ITT,  Electron  Tube  Division,  Fort  Wayne,  Indiana 


Disclosure  of  Invention  to  Others  Continued 


0.  L.  Gaddy,  Professor  & 
Associate  Head,  Dept,  of 
Elect.  Engrg.,  University 
of  Illinois,  Urbana 

Verbal 

January  1973  & during 
construction 

No 

S.  Kaiser,  Research  Assistant 
University  of  Illinois 

Assisted  in 
construction 

Summer  1972-January  1973 

No 

J.  Merkel o.  Manager 
Hewlett  Packard  Co. 

Participated  in 
design 

Fall  1972 

No 

Santa  Clara,  California 


This  invention  was  | X | was  not  □ made  while  working  on  research 
sponsored  by  an  agency  outside  the  University.  If  so, 

| x | Name  of  Sponsoring  Agency:  Office  of  Naval  Research 

| X | Description  of  contract:  Research  in  Infrared  and  Optical  Quantum 

Electronics 

1 x|  Number:  N00014-67-A-0305-0017 

If  decision  of  Patent  Committee  is  to  recommend  release  of  interest  of 
the  University  in  this  invention,  the  following  is  recommended: 

n Outright  release  to  sponsoring  agency,  retaining  no  rights 

for  inventor. 


Release  to  inventor,  with  recognition  of  contractual 
obligations  of  the  University  to  sponsoring  agency. 

Other. 


Respectfully  submitted, 

r . t /\  , 

' A f t 1 / /' 

Inventor  Date 


Inventor 


It  is  certified  that  the  statements  made  herein  are  correct,  to  the  best 
of  my  knowledge  and  belief. 


^AJ 


Department  Head,'  J 
or  Other  Administrative 
Officer  J 
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Who  Have  Knowledge 
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High  Speed  Opto-Electronic  Sampler  and  High  Resolution 
Sampling  Photomultiplier 

Design  Criteria  and  Results 

The  high  resolution  sampling  photomultiplier  consists  of  a high  speed 
head,  designed  to  sample  photoelectrons  at  its  own  photocathode,  and  of  a 
conventional  electron  multiplier  section  with  no  special  requirements.  The 
sampling  head  thus  converts  any  conventionally  designed  electron  multiplier 
into  a photomultiplier  capable  of  high  time  resolution  of  periodic,  low-level 
optical  signals.  The  first  complete  system  was  built  and  tested  on  light 
signals  of  “300  psec  rise  (3  x 10  sec,  10  to  90%  amplitude)  fully  resolving 
this  time  dependence;  sufficiently  short  pulses  were  not  available  to  test  the 
device  to  its  limit  of  time  resolution.  Whereas  the  sampling  structure  was 
driven  up  to  10  GHz  bandwidth  for  a corresponding  pulse  resolution  of  “35  psec, 
the  maximum  bandwidth  requirement  on  the  auxiliary  instruments,  including  the 
electron-multiplier  structure  was  only  100  KHz  (“3.5  x 10~^  sec  rise  time). 

Figure  1 shows  the  assembly  of  the  sampling  head  and  the  modified  electron 
multiplier.  The  electron  multiplier  is  that  of  an  RCA  7102  photo-tube.  The 
multiplier  is  installed  in  a new  glass  press  and  housing  which  is  flanged  in 
order  to  accommodate  the  sampling  head.  Figure  2 shows  the  detailed  components 
of  the  sampling  head. 

The  flange  which  is  labeled  with  0-ring  groove  //I  accommodates  a glass 
window;  the  housing,  #2,  is  made  of  cold  rolled  steel;  high  frequency,  500 
feedthroughs,  #3,  connect  to  a perforated  strip,  #6,  which,  in  combination  with 
insert  #5,  forms  a 500  microstrip  transmission  section  terminated  in  a 500  load; 
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insert  if 5 is  gridded  with  a mesh.  It 4,  for  photoelectron  ejection;  moreover, 
accurate  positioning  of  insert  #5  with  respect  to  strip  #6  is  accomplished 
with  a gauged  spacer,  #7.  The  assembly  of  these  components  mounts  on  the 
flange  of  the  multiplier  section  as  shown  in  Figure  1.  Although  many 
construction  details  are  important  in  producing  this  device,  the  most 
significant  feature  of  this  design  is  that  the  opto-electronic  sampling 
section  (photocathode,  strip,  ground  plane,  feedthroughs)  form  a broadband, 
non-dispersive  electrical  system;  this  was  accomplished  with  an  opto- 
electronically  adapted  microstrip  design  suitable  for  ultra  high  vacuum 
construction.  Scaling  of  the  microstrip  section  for  electromagnetic  matching 
was  carried  out  following  considerations  applicable  to  broadband  microstrip 
design  for  reflectionless  and  dispersionless  fundamental  mode  propagation; 
a 300  ym  strip-ground  plane  spacing  and  a 1.58  mm  strip  width  were  chosen  as 
suitable  geometric  parameters;  a cesium  antimony  photocathode  was  used  in 
all  experiments. 

When  properly  biased  and  driven  with  ultrashort  pulses.  Figure  3,  the 
strip-grid  section  (//6  and  #4)  is  capable  of  sampling  photoelectrons  with  a 
time  resolution  which  depends  both  on  sampling  pulse  amplitude  and  pulse 
shape.  Conventional  photoelectron  ballistics  analysis  predicts  that  currently 
available  technology  of  sampling  electronics  is  capable  of  achieving  pico- 
second resolution  and  possibly  reach  the  limiting  resolution  which,  on  the 

-11  -12 

basis  of  quantum  mechanical  arguments,  should  lie  in  the  10  - 10  sec 

range. 
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Figure  1.  High  speed  opto-electronlc  sampling  photomultiplier.  Assembly 
of  sampling  head  and  electron  multiplier. 


6 

Optical  pulses  of  ~300  psec  risetime  (tested  at  high  intensity  with 

a fast  photodiode  and  a sampling  oscilloscope)  were  derived  from  a mode- 

locked  He-Ne  laser  and  detected  with  the  reported  device.  Only  ~7  V, 

250  psec  wide  sampling  pulses  with  a 2.2  V back-bias  were  necessary  to 

achieve  full  time  resolution  of  the  mode- locked  laser  pulses.  Figure  4 

shows  opto-electronically  sampled  light  pulses  as  displayed  on  an  ordinary 

oscilloscope.  The  equivalent  horizontal  scale  is  200  psec/div 
-12 

(200  x 10  sec/div).  The  risetime  (10%  to  90%)  of  the  pulse  shown  is 
*300  psec,  as  measured  independently,  and  the  fall  time  is  identical, 
demonstrating  for  the  first  time  (to  our  knowledge)  the  symmetry  of  mode- 
locked  laser  pulses.  Shorter  optical  pulses  will  be  needed  to  establish 
the  true  limit  of  the  device. 
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Figure  4.  Optoelectronically  sampled 
mode-locked  laser  pulses  displayed  on 
an  ordinary  oscilloscope.  Horizontal 
equivalent  scale:  200psec/div  (200 

•x  10“l2sec/div)  . 
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A HIGH  SPEED  SAMPLING  PHOTOMULTIPLIER* 


Henry  Merkelo,  Janes  J.  Wiczer 

Quantum  Electronics  Research  Laboratory,  University  of  Illinois, 

Urbana,  IL,  61801 


ABSTRACT 

A sampling  approach  to  low-level  opti- 
cal signal  detection  is  reported  and  the 
advantages  of  sampler  based  instruments 
for  the  analysis  of  complex,  high-harmonic 
content  signals  are  re-examined  for  opto- 
electronic considerations.  A sampling 
head,  constructed  for  high  speed  sampling 
and  processed  for  an  S-4  photosensitivity 
is  described.  The  sampling  head,  adapted 
to  an  electron  multiplier  section,  is 
tested  for  sequential  sampling  on  HeNe, 
mode-locked  laser  pulses  showing  good  re- 
solution and  signal  to  noise  ratio.  The 
influence  of  parameter  variation  on  the" 
sampling  function  is  assessed  on  the  basis 
of  photoelectron  ballistics;  discussion 
includes  optimum  and  limiting  resolutions 
which  extend  into  the  picosecond  regime. 


The  development  of  photodetectors 
which  are  capable  of  microwave  modulation- 
frequency  response  and  which  operate  in 
the  visible  or  near-infrared  portion  of 
the  optical  spectrum  has  received  a con- 
siderable amount  of  attention.  The  need 
for  detectors  in  this  category,  particu- 
larly for  low-level  signal  monitoring, 
exists  in  both  communication  and  analyt- 
ical instrumentation  areas.  However,  in 
spite  of  the  significant  progress  that 
has  been  made,  low-level  optical  phenom- 
ena of  picosecond  (10" 12  sec)  and  even 
tens  of  picoseconds  duration  cannot  be 
monitored  routinely  today  as  it  is  pos- 
sible in  some  cases  in  microwave  signal 
processing  (1,2).  In  relation  to  analyt- 
ical instrumentation,  where  not  only  good 
signal-to-noise  ratio  and  high  bandwidth 
are  required  for  satisfactory  detection, 
pulse  fidelity  requirements  limit  time 
resolution  to  several  hundred  picoseconds 
if  signal  deconvolution  is  to  be  avoided; 
the  shortcomings  of  the  current  detector 
status  are  particularly  acute  in  view  of 
the  existence  of  tunable  optical  excita- 
tion sources  capable  of  continuously  e- 
mitting  pulses  of  picosecond  duration, 
suitable  for  many  studies  and  applica- 


tions. Moreover,  high  time  resolution 
studies  that  are  currently  being  carried 
out  at  high  intensities,  by  the  necessity 
of  indirect  detection  methods,  do  not  al- 
ways establish  low  intensity  behavior. 

The  application  of  sampling  princi- 
ples have  long  been  used  on  periodic  wave- 
forms to  achieve  wide  bandwidths  in  oscil- 
loscopes, phase-locked  loops,  random  sam- 
pling voltage  detectors  and  other  rf  sys- 
tems. Resolution  of  10“l°  sec  is  routine- 
ly achieved  and,  with  thin-filin  samplers, 
bandwidths  better  than  18GHz  can  be 
reached  (2). 

A sampling  approach  to  low-level 
optical  signal  processing  is  reported  and 
the  advantages  of  sampler  based  instru- 
ments for  complex,  high-frequency  signal 
analysis  are  re-examined  for  optoelectron- 
ic consideration. 

A sampling  system  generally  consists 
of  a broadband  sampler,  an  integrating 
section,  an  if-amplifier  and,  for  display 
purposes,  some  clock  network  correlating 
high  speed  events  to  the  if  representa- 
tion of  the  original  signal  such  that  the 
output  spectrum,  S0(w),  is  the  convolu- 
tion of  the  signal  spectrum,  S(w),  with 
the  sampling  function  spectrum,  A(u>).  mul- 
tiplied by  the  frequency  response  of  the 
integrator-amplifier  block.  The  low-pass 
portion  of  this  spectrum  is  that  of  a 
periodic  signal  similar  in  spectrum  to 
that  of  the  input  signal  but  with  modula- 
tion harmonics  niu  = 2irnf0  substituted 
for  nui  = 2unf  where  1'0,  t are  fundamental 
frequencies  of  the  output  spectrum  S„(w) 
and  the  original  signal  S(u)  respectively. 
Since  spectral  components  also  appear  cen- 
tered about  each  harmonic  of  the  sampled 
output,  the  harmonic  sampling  extension 
of  the  classical  sampling  theorem  re- 
quires that  mfn<^/2  where  m is  the  high- 
est harmonic  component  in  the  high- 
frequency  waveform.  For  example,  sequen- 
tial sampling  will  take  place  with  high 
signal  fidelity  reproduction  if  the  sam- 
pling function  consists  of  equal  harmonic 
components  and  unity  efficiency.  Then, 
a signal 
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S(id)  * ? (Ai,sinkMt  + B,,co. skint ) 

k»i  k K 

when  sampled  at  a fundamental  frequency , 
f^,  and  low-pass  filtered,  becomes 

SQ ( cu ) ■ ? ( A^sinkAuit  + BkcoskAut) 

where  Au>»  2Tt(T\  + f).  If  the  sampler  devi- 
ates from  the  ideal  situation  of  equal 
harmonic  components,  the  reconstructed 
component  amplitudes  will  be  also  modi- 
fied by  the  rf  frequency  response  of  the 
sampler.  It  is  interesting  to  note  that 
the  requirement  for  equal  harmonic  ampli- 
tudes in  the  sampling  function  on'y  ap- 
plies to  the  terms  that  exist  in  the  sig- 
nal, i.e.  the  first  m terms.  Any  addi- 
tional components  in  the  sampling  func- 
tion do  not  serve  a useful  function. 


Adaptation  of  these  fundamental  con- 
cepts to  photoelectric  devices  requires 
application  of  broadband  design  consid- 
erations of  components  which  are,  at 
least  partly,  assembled  and/or  processed 
under  vacuum.  It  is  shown,  in  partic- 
ular that  tolerances  significantly  affect 
both  the  sampling  efficiency  and  the 
form  of  the  sampling  function.  In  this 
sense,  optoelectronic  sampling  is  ana- 
logous in  principle  to  the  sampling 
used  in  conventional  sampling  electronic 
instruments,  but  implementation  of  the 
sampling  of  photoelectrons  is  techno- 
logically different  from  electronic 
sampling . 


Given  that  the  sampling  process  for 
photoelectrons  consists  of  electrons  tra- 
versing a drift  space  of  dimension  Ax 
under  the  influence  of  a sampling  voltage, 
V(t),  the  sampling  function  A ( t ) is 
arrived  at  on  the  basis  of  the  ballistics 
for  photoelectrons  which  are  distributed 
in  initial  velocity  according  to  a func- 
tion, g\(vQ),  which  is  dependent  on  the 
type  of  photocathode  material  and  on  the 
wavelength,  X,  of  the  optical  signal.  A 
photoelectron  is  considered  sampled  if, 
as  a result  of  the  drift  space  potential, 
V(t)-Vr.  where  Vr  is  a constant  retarding 
potential  (at  least  equal  to  the  most 
energetic  photoelectrons),  the  electron 
traverses  the  drift  space,  Ax.  The  dis- 
placement of  photoelectrons  is  given  by 


x(t)-  U V(t)dt2 
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determining  the  duration  of  the  sampling 
function,  A(t).  In  fact,  due  to  nonzero 
initial  velocities  ol'  photoelectrons,  the 
potential  V(t)  is  affecting  the  position 
of  electrons  well  before  V(t)-V_=0.  More- 
over, the  sampling  function  A(t),  bears 
no  resemblance  to  V(t)  since  application 
of  V(t)  at  t=0  can  effect  the  transversal 
of  Ax  by  photoelectrons  emitted  at  t<0. 

The  sampling  function  A(t),  the  sampling 
voltage  V(t),  and  the  actual  waveform  of 
sampled  photoelectrons  are  all  shifted  in 
time  with  respect  to  each  other  but  retain 
the  enumerated  chronological  order.  For 
the  simple  case  of  Ka(vq)=  fi(vo-voo) 
where  6 is  the  Dirac  delta  function  (i.e. 
monoenerget ic  photoelectrons  correspon- 
ding to  initial  velocity  v00), 

A(-t,  _ 0 if  xt(t’)<  Ax  for  all  t’ 

1 ' A if  xt(t')^  Ax  for  some  t' 

where  x^(t')  designates  displacement  of 
those  photoelectrons  which  are  emitted 
from  the  photocathode  at  time  t;  note 
that  it  is  irrelevant  at  what  time,  t', 
the  electron  traverses  the  drift  space  Ax. 
That  is,  if  V(t)  is  sufficiently  large  to 
effect  displacements  such  that  xt(t')>Ax 
for  some  t’  and  V(t)  is  monotonical ly 
increasing  and  then  monotonically  decreas- 
ing, then  A(t)  is  a square  pulse;  the 
width,  x,  of  A(t)  is,  however,  a function 
of  the  form  of  V(t).  When  the  initial 
velocities  of  photoelectrons  obey  a rel- 
atively narrow  distribution  function,  A(t) 
will  be  still  flat-topped  but  with 
rounded  corners. 

Extensive  analysis  aimed  at  deter- 
mining the  sampling  function,  A(t),  and 
the  sampling  efficiency  for  various  geo- 
metrical conf  igurat ions  and  velocity  dis- 
tributions, gx(vG),  show  that  a practical 
optoelectronic  sampling  device  effecting 
a sampling  time  T=80psec  can  be  routinely 
achieved.  For  illustration,  adapting  a 
specially  designed  sampling  head  on  a 
conventional  electron  multiplier  cage 
(RCA  7102),  HeNe,  mode-locked  laser 
pulses  were  monitored  in  the  optoelectro- 
nic sequential  sampling  mode  and  dis- 
played on  an  ordinary  (real  time)  oscillo- 
scope as  shown  on  the  photograph  of  Fig- 
ure 1 (equivalent  horizontal  scale:  200 

psec/div).  For  a square  wave  sampling 
function,  A(t),  the  magnitudes  of  the 
harmonics  of  order  k are  proportional  to 


where  q,m  are  the  electronic  charge  and 
mass  respectively  and  v0,  xQ  are  the 
initial  velocity  and  displacement  respec- 
tively. Unlike  in  conventional  electronic 
sampling  where  the  duration  of  the  sam- 
pling bridge  forward  bias  and  therefore 
of  the  sampling  function  is  controlled 
by  the  time  the  strobe  pulse  exceeds  a 
fixed  reverse  bias,  the  sampling  voltage 
on  photoelectrons  enters  nontrivially  in 


(sin 


kwx 
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Typical  HeNe,  mode-locked  laser  pulses 
consist  of  only  approximately  15  harmonics 
and  since  the  value  of  the  amplitude  of 
the  fifteenth  term  of  the  sampling  func- 
tion drops  to  only  0.987  of  the  amplitude 
of  the  fundamental,  it  is  consistent  that 
good  sampled  waveforms  exhibiting  the 
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expected  symmetry  observed.  Since  only 
the  sampling  head  is  subjected  to  the 
input  signal  and  sampling  high-frequen- 
cies and  all  the  amplification  takes 
place  at  relatively  low-frequencies,  the 
performance  of  the  sampling  system  is  not 
dependent  on  the  gain-bandwidth  limita- 
tion of  conventional  amplifiers.  The 
signal-to-noise  ratio  is  significantly 
better  than  in  any  fast  photomultiplier/ 
sampling  oscilloscope  combination  and 
ample  gain  is  available  for  detecting 
low-level  radiation. 

Unlike  many  photoelectric  gated  de- 
vices, which  operate  at  high  voltages, 
low-level  signal  design  suitable  for 
broadband  and  high  sampling  rate  opera- 
tion was  stressed  in  this  development. 

The  resulting  disadvantage  is  that  the 
system  is  sensitive  to  initial-  velocity 
distribution,  gx(vQ),  of  photoelectrons 
and  therefore,  the  sampling  function, 
A(t),  is  wavelength  dependent.  For 
many  photocathode  materials,  best  reso- 
lution will  be  obtained  on  the  red  end 
of  the  spectrum.  The  choice  of  cesium 
antimonide  photocatt.ode  and  6328  A op- 
eration, although  fortuitous  here,  is 
particularly  suitable  for  obtaining  good 
resolution.  In  spite  of  the  dispersion 
effects  introduced  by  initial  photoel- 
ectron velocities,  extensive  computations 
show  that  with  high  precision  design  and 
select  use  of  thin  and  thick  film  tech- 
nologies, a sampling  width,  r,  of  less 
than  lOpsec  can  be  realized  with  cur- 
rently available  methods. 


figure  1.  Optoelectronical ly  sam- 
pled laser  pulses  displayed  on  an 
ordinary  oscilloscope.  Horizontal 
equivalent  scale:  200  psec/div. 


Finally,  the  possibility  of  large 
area  photocathodes  in  high  speed  opto- 
electronic sampling  devices  has  also  been 
explored  with  indication  that  significant 
improvements,  with  respect  to  any  other 
known  designs  suitable  for  high  speed 
application,  are  realizable.  Again,  this 
bears  on  the  fact  that  arrival  of  photo- 
electrons  to  the  first  dynodo  is  already 
parr  of  the  slow  structure  operation  and, 
therefore,  does  not  require  stringent 
timing  considerations. 
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OPTOELECTRONIC  SAMPLING— I.  PARAMETERS  FOR  PICOSECOND  RESOLUTION 
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ABSTRACT 

Saapling  of  photoelectrons  by  rapidly  varying  electric  fields  is  for- 
taulated  for  configurations  suitable  .for  high-speed  detection  of  optical 
signals.  The  dynamics  of  photoelectrons  naar  the  surface  of  the  emitter  is 
taken  into  account  in  the  modeling  and  the  numerical  analysis  designed  to 
characterize  the  saapling  process.  The  influence  of  delays  in  the  photoelec 
trie  eaission  on  tha  resolution  and  efficiency  of  optoelectronic  sampling 
is  given  particular  attention.  It  is  shown  that  such  sampling  can  fora  the 
basis  for  high-gain  detectors  capable  of  picosecond  resolution  of  optical 
signals.  Sequential  sampling  is  especially  adaptable  to  studies  involving 
the  periodic  emission  of  continuously  operated,  mode-locked  lasers. 

Quantum  efficiencies  being  comparable  to  those  existing  in  conventional 
pho Cod a tec tors,  the  sampling  approach  is  suitable  for  detection  at  low 
signal  levels;  in  the  limit,  picosecond  resolution  counting  of  photons 
should  be  possible.  Moreover,  it  is  demonstrated  that  high  resolution 
optoelectronic  sampling  is  effective  at  low  sampling  voltages  and,  there- 
fore, can  be  made  compatible  with  conventional  electronic  signal  processing 
instruments  for  applications  in  low  light-lavel  picosecond  spectroscopy 
and  optical  pulse  dispersion  studies. 
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I.  INTRODUCTION 

Optical  studies  requiring  picosecond  resolution  have  taken  a new 
direction  with  the  advent  of  broadband,  mode-locked  lasers.  The  new 
methods  for  generating,  controlling,  and  monitoring  picosecond  optical 
events  have  led  to  significant  results  in  photophysics,  photochemistry, 
photobiology,  and  optical  communication  studies  and  applications.*  Cur- 
rently, however,  the  techniques  for  generating,  amplifying,  frequency 
doubling,  isolating  and  chirping  picosecond  optical  pulses  are  far  more 
advanced  than  the  methods  used  for  accurately  and  conveniently  monitoring 
such  processes.  Whereas  a wide  variety  of  sources  exists  for  generating 
ultrashort  pulses  of  various  durations  and  of  adjustable  intervals,  detec- 
tion methods  suitable  for  picosecond  resolution  of  optical  events  are  not 
yet  sufficiently  routine  for  widespread  usage.  In  particular,  difficulties 
arise  at  low  signal  levels  and  with  incoherent  radiation  when  frequency 
doubling  or  multiphoton  interactions  are  ineffective. 

The  development  of  photodetectors  that  are  capable  of  microwave  band- 
width modulation  response  and  that  operate  in  the  visible  or  near  infrared 
portion  of  the  optical  spectrum  has  received  considerable  attention  as 

evident,  for  example,  from  the  extensive  bibliography  documenting  the 

2 

review  written  by  Anderson  and  McMurtry.  It  is  recognized  that  the.  need 
for  detectors  in  this  category,  particularly  for  low-level  signal  monitor- 
ing, exists  in  both  communication  and  analytical  instrumentation  areas. 
However,  in  spite  of  the  significant  progress  that  has  been  made,  some 
low-level  optical  phenomena  of  picosecond  and  even  tens  of  picoseconds 
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duration  cannot  be  recorded  today  with  adequate  time  resolution.  In  re- 
lation to  analytical  Instrumentation,  where  not  only  good  signal-to-noise 
ratio  and  largp  bandwidth  are  required  for  satisfactory  detection  of  low 
signal  levels,  pulse  fidelity  requirements  limit  time  resolution  of 
optical  waveforms  to  several  hundred  picoseconds  if  signal  deconvolution 
is  to  be  avoided.  The  shortcomings  of  the  current  detector  status  are 
particularly  acute  in  view  of  the  existence  of  tunable  optical  excitation 
sources,  capable  of  continuously  emitting  pulses  of  picosecond  duration, 
suitable  for  many  studies  and  applications.  Moreover,  some  high-time 
resolution  studies  that  are  currently  being  carried  out  at  high  intensities 
do  not  always  determine  low  intensity  parameters.  This  is  particularly 
the  case  with  luminescence  in  semiconductors  and  biological,  optically 

active  materials.  For  instance,  saturation  occurs  at  modest  illumination 
-2 

levels  (<  1 mW  cm  ) in  biological  systems  of  the  important,  photosynthetic 

3 

group  in  which  the  dynamics  of  fluorescence  has  picosecond  characteristics. 

Other  than  indirect  methods  , means  do  not  exist  for  studying  picosecond 

dynamics  of  reactions  taking  place  in  these  materials  at  illumination 

3 

intensities  below  the  luminescence  saturation  limit  . 

The  aim  of  this  study  is  to  examine  the  basic  characteristics  of  photo- 
detectors in  a new  configuration  suitable  for  high-time  resolution  and 
adaptable — in  sampling  operating  modes— to  high-gain  multipliers  and  signal 
averaging  instruments.  Motivation  for  this  work  was  provided  by  the  success- 
ful preliminary  testing  of  the  basic  optoelectronic  sampling  configuration 

> 

which  was  previously  reported^'**;  the  operation  of  the  sampling  prototype 
was  compared  to  the  performance  of  the  most  developed,  high-speed,  commercial 
phototube  suitable  for  low-level  signal  detection,  which  at  that  time  was 
the  Sylvania  #502  Static  Crossed  Field  Photomultiplier  rated  at  beyond  3.0 
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GHz  bandwidth.  Even  without  the  application  of  thin-  and  thick-film 
technologies  in  the  construction  of  the  prototype,  the  advantages  of 
optoelectronic  sampling  were  demonstrated  on  both  signal  quality  and 
signal  fidelity.  Since  electric  ringing  in  the  waveform  trailing  edge 
is  frequently  encountered  in  high-speed  photodetectors — rendering  signal 
comparison  or  deconvolution  analysis  unreliable — signal  fidelity  will 
be  an  Important  consideration  in  this  text. 

It  is  well  known  that  the  evolution  of  sampling  principles  enabled 
significant  advances  to  take  place  in  the  design  of  wide-band  oscillo- 
scopes, phase- locked  loops,  random  sampling  voltage  detectors  and  other 
microwave  systems  suitable  for  periodic  waveform  analysis.  Resolution 
of  a hundred  picoseconds  is  routinely  achieved  in  conventional  electronic 
devices  designed  for  sampling  mode  operation^;  with  thin-film  samplers, 

g 

bandwidths  better  than  20  GHz  have  been  attained  . Sampling  methods 
sidestep  both  the  gain-bandwidth  and  display-size  limitations  of  auxiliary 
devices  by  translating  the  input  signal  to  a much  lower  frequency  domain 
where  conventional  methods  for  amplifying,  recording  or  displaying  are 
more  effective.  For  these  reasons,  a high-speed  sampling  approach  to  low- 
level,  optical  signal  monitoring  has  been  devised  and  the  advantages  of 
sampler  based  instruments  for  complex,  high-frequency  signal  analysis 
are  evaluated  in  optoelectronic  applications.  Although  only  coherent,  se- 
quential sampling  of  periodic  waveforms  will  be  formally  treated  here, 
many  other  aspects  of  signal  processing  by  sampling  methods  will  remain 
within  the  context  of  this  paper. 
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II.  HIGH-SPEED  OPTOELECTRONIC  SAMPLING 
Generally,  a sampling  system  consists  of  a broadband  sampler,  an 
integrating  section  which  stretches  the  samples  into  a continuous  signal, 
an  intermediate  frequency  (IF)  amplifier  which  provides  signal. gain,  and 
some  clock  network  correlating  high-speed  events  to  the  IF  simulation  of 
the  original  signal^.  Thus,  a sampling  system  can  be  viewed  as  a series 
of  separate  processes,  in  which  the  final  output  can  be  determined  by  band- 
width considerations.  Photoelectric  devices,  being  current  sources,  il- 
lustrate particularly  well  the  direct  dependence  of  signal  level  on  load 
resistance  and  the  inverse  relationship  between  signal  amplitude  and  band- 
width. For  these  reasons,  the  advantages  of  down-converting  operating 
frequencies  afforded  by  sampling  are  immediately  apparent. 

The  family  of  optoelectronic  devices  related  to  this  study  includes 
sampling  photomultipliers,  sampling  photon  counters,  adapting  sampling  heads, 
array  detectors,  as  well  as  demultiplexing  and  decoding  samplers,  all 
suitable  for  high-speed  operation.  In  this  text,  the  analysis  will  be 
limited  only  to  the  basic  sampling  configuration  which  is  common  to  all 

designs.  The  sampling  geometry  consists  of  an  electrically  broadband, 

* 9 

vacuum-spaced,  suspended  microstrip  in  which  photoelectrons  are  sampled  by 
rapidly  varying  accelerating  electric  fields  created  by  ultrashort  electrical 
sampling  pulses.  These  pulses,  which  are  generated  in  thin-film,  ultra- 
broadband circuits,  propagate  along  the  photoemisslve  thin-film  microstrip 

% 

line.  The  coordinates  necessary  for  the  analysis  of  optoelectronic  sampling 

*For  achieving  higher  tolerances,  the  vacuum-spaced  configuration  can  be  re- 
placed with  a microchanneled  design.  Precision  fabrication  of  photoelectri- 
cally  active,  mlcrostrlp  samplers,  which  involves  both  thin-  and  thick-film 
technologies,  will  be  described  elsewhere. 
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•re  indicated  in  Figure  1.  Photoelectrons  created  at  the  emitter  E are  swept 
through  the  drift  space  D of  dimension  Ax.  The  surface  and  near-surface  ki- 
netics of  these  charges  which  are  influenced  by  the  energies  of  quantum 
origin  and  by  forces  induced  by  propagating  electric  fields,  significantly 
affect  the  resolution  and  the  efficiency  of  sampling.  For  these  reasons, 
the  analysis  of  high-speed  optoelectronic  sampling  takes  into  account  the 
following  considerations : 


1)  Electrons  originating  at  the  emitter  E are  accelerated  in  the 
drift  space  D by  rapidly  varying,  electric  fields  E(t,z)  of  the  form 
(v(t,z)  - VrJ/Ax  where  is  the  retarding  potential  which  is  transiently 
counteracted  by  the  propagating  potential  v(t,z)  traveling  along  the 


z-directlon  with  time  dependence  t.  Although  constituting  moderate  potentials 
4.1 -11 j — j — - — -i jECt.?2).  Df  500  V cm~*  psec* 


these  fields  are  allowed  to  exceed  rates  of  change 


An  example  of  such  rapidly  varying  potentials  v(t,zo)  (where  zq  is  an  arbi- 
trary point),  generated  for  optoelectronic  sampling, is  illustrated  in 
Figure  2.  These  waveforms  are  represented  quite  well  by  the  following 


relations: 


0 for  t < 0 

v(t,z  ) - V sin2(£t)  for  0 < t < T 
o o r — — 

0 for  t > r 


where  V and  T are  convenient  parameters  for  numerical  analysis.  Superimposed 
o 

on  the  oscilloscope  trace  of  Figure  2,  the  dotted  line  is  a graph  cf  (1)  for 

V - 14.2  V and  r ■ 68  psec.  The  generation  of  such  broadband  signals  is 
o 

essential  for  realizing  effective  high-speed  sampling.  However,  we  will  show 

> 

on  the  basis  of  the  properties  of  surface  emission  of  photoelectrons  that  the 
characterization  of  the  actual  optoelectronic  sampling  function  is  more 
informative  than  a detailed  discussion  of  the  fields  responsible  for  the 
sampling,  since,  for  example,  the  waveform  of  Figure  2 can  give  values  for 
sampling  resolution  ranging  from  a hundred  to  a few  picoseconds  depending  on 


the  choice  of  other  parameters. 
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2)  The  displacement  x of  photoelectrons  in  the  drift  space  depends 


1 

j 
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on  their  x-dlrected,  initial  velocities  v^.  Corresponding  to  each  wavelength 
X of  the  incident  radiation,  the  photoelectrons  emitted  by  each  material  will 
be  distributed  in  velocity,  at  x ■ 0,  according  to  a function  designated  as 
g^(vx).  Experimentally  determined  velocity  distributions  will  be  used  in 
specific  numerical  examples. 

3)  Whereas  a reasonable  amount  of  experimental  data  exists  on  (v^) , 
the  time  of  emission  y of  photoelectrons,  which  establishes  the  response  time 
of  the  emitter,  remains  sparsely  documented.  Without  even  taking  into  account 
the  long  standing  problem  of  the  time  dependence  of  the  probability  of  photo- 
electric emission  on  time  scales  on  which  the  oscillations  of  the  optical 
field  are  resolved,  as  considered  for  example  by  Handel  and  Meltzer^,  the 
response  time  of  a finite  thickness  photoactive  surface  can  be  significantly 
affected  by  lattice  scattering,  pair  production,  formation  and  destruction  of 
exciton  and  impurity  states,  and  free  carrier  lifetimes.  For  Instance,  by 
assuming  that  lattice  scattering  determines  the  escape  depth  in  all  useful 
photoemitters  in  which  the  threshold  for  pair  production  is  large  compared 
with  the  electron  affinity,  it  is  possible  to  construct  arguments  which  give 

an  upper  range  for  respoz  e time  such  that  1 < y < 10  psec  for  surface  thick- 

o 

nesses  on  the  order  of  500A.  Photoemittlng  layers  substantially  thinner 
than  average  escape  depth  will  be  characterized  by  smaller  y's*  For  opto- 
electronic sampling  considerations,  it  is  Important  to  know  more  than  just 
an  average  response  time.  In  particular,  since  ejection  velocities  are 
taken  into  account  in  detail,  ideally,  the  time  evolution  of  photoelectric 
emission  should  be  established  for  each  velocity.  The  time  evolution  of 
emission  should  be  principally  controlled  by  the  escape  mechanism  and  the 
Interactions  concomitant  with  escape.  On  that  basis,  it  can  be  concluded 

from  the  observations  of  Methfessel^  that  the  low  energy  tail  of  8^(vx) 
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will  have  a greater  time  of  emission  dispersion  chan  Che  high  energy  group 
and  thac  this  dispersion  will  be  smaller  for  a fronC-illuminated  Chan  for 
a back- illuminated  (semlcransparenc)  photocathode.  AC  present,  only  an 
upper  limit  of  approximately  30  psec  can  be  deduced  for  y from  the  existing 

literature^. 

4)  Photoelectric  emission  is  enhanced  by  electric  fields  in  the  kV  cm~* 

13,14 . „ 

range  ; of-  course,  only  static  measurements  exist.  Similar  quantum 
efficiency  Improvements  are  expected  with  rapidly  propagating  fields  and, 
moreover,  both  g. (v  ) and  y are  expected  to  be  altered  by  sampling  fields 
which,  for  the  cases  considered,  may  be  as  high  as  3kV  cmT^  It  is  expected 
that  both  resolution  and  sampling  efficiency  would  be  favorably  affected 
by  such  field  enhancement  but,  due  to  the  lack  of  data,  this  consideration 
will  not  be  taken  into  account  quantitatively  at  this  time. 

The  waveforms  shown  in  Figure  3 are  useful  in  discussing  optoelectronic 
sampling.  Only  sequential  sampling  is  considered  in  this  analysis.  The 
intensity  envelope  of  optical  radiation  (incident  on  the  detector)  which  is 
Illustrated  in  Figure  3A,  is  represented  by  a periodic  waveform  f^(t)  of 
period  T.  Electrons  injected  into  the  drift  space  are  tinder  the  Influence 
of  fields  E(t,z),  of  period  T + At,  created  by  propagating  potentials 
(v(t,z)  - V^].  The  sampling  function  fs<t),  which  effects  sampling  of  the 
form  f. (t)f  (t),  is  constructed  on  the  basis  of  the  effectiveness  of  the 
fields  E(t,z)  to  extract  electrons  from  the  drift  space  as  a function  of 
time;  l.e.,  a photoelectron  emitted  as  a result  of  a photon  Interaction  at 
some  time  t is  considered  to  have  been  sampled  if  it  has  completely 
traversed  the  dimension  Ax  of  the  drift  space  within  one  period.  Note 
that  the  time  reference  is  that  of  the  photon  striking  the  emitter  and  not 
the  potentially  delayed  time  t + y,  nor  is  it  the  time  at  which  the  dis- 
placement Ax  is  completed.  In  fact,  unlike  in  a conventional  high-speed 
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photodiode,  the  time  t^  at  which  x - Ax  is  entirely  irrelevant  as  long  as 
tAy-t  < T.  The  duration  of  any  one  pulse  of  — full  width  at  half 

maximum  — is  designated  by  t and  is  shown  in  Figure  3C.  It  is  determined 
by  establishing  the  duration  of  sampling  effected  by  E(t,z)  at  a single  point 
zq.  This  type  of  stationary  analysis  is  valid  for  any  photoactive  length, 

Az,  such  that  Az  « tc  where  c is  the  propagation  velocity  of  E(t,z).  For 
all  such  cases,  voltages  and  fields  are  simply  designated  by  v(t),  as  in 
Figure  3B,  and  E(t),  respectively.  The  effect  of  progressively  delaying 
each  sample  by  an  increment  At  with  respect  to  the  relative  position  of  the 
preceding  sample  on  f^(t)  results  in  a replication  of  the  original  signal 


f.(t)  by  f (t)  such  that  f (t * ) ■ G f . (t)  where  G is  a constant  and  t*  - 
1 o o i 

T 

— t is  a scaling  transformation  illustrated  in  Figure  3D.  Although  over- 
simplified, this  conceptualization  asserts  that  bandwidth  requirements  for 
auxiliary  subsystems  (apart  from  the  principal  sampler)  are  reduced  by  a 
factor  of  — ■ where,  generally.  At  « T. 

It  is  important  to  note  that  the  instantaneous  sampling  is  of  the  form 

f.(t)f  (t)  but  that,  prior  to  amplification  and  low  pass  filtering,  the 
is*' 

amplitude  contribution  to  fo(t),  corresponding  to  any  period  n is  given  by 

/ f (t)  f ft)  dt  . ; W 

(n-l)T  1 

The  intuitive  understanding  is  further  enhanced  in  the  limiting  case  for 

which  every  period  n of  f_(t)  is  a Dirac  delta  function  3(t-t  ).  Sampling 

s n 

could  be  considered  perfect  since  (2)  becomes 


ni 

/ ^(t)  6 (t  - t ) dt  - f (t  ) 

(n-l)T  1 n in 


Clearly,  the  requirement  that  f^(t)  be  a slowly  varying  function  during  the 

time  t - — < t < t + -r,  where  T again  designates  the  duration  of  any  one  pulse 
n 2 n l 


of  f (t),  will  be  an  important  consideration.  Accurately  characterizing  f (t) 
and,  therefore,  establishing  x and  the  efficiency  of  sampling  for  various  para-. 


meters  is  clearly  the  key  factor  to  be  considered  in  appraising  the  effective- 
ness of  optoelectronic  sampling.  The  procedure  for  obtaining  f (t)  involves 

s 

solving  the  photoelectron  displacement  equation,  which  is  of  the  form 


a a' 

x(t,a)  - J da’  / v(a")  da" 


qV  a2 
n r 

2mAx 


+ . 


v a 


(4) 


where  q and  m are  the  electronic  charge  and  mass,  respectively;  t is  the 

time  of  light  interaction  at  the  emitter;  a is  a time  reference  such  that 

the  origin  of  a corresponds  to  the  light  interaction  time  t (i.e.,  x(t,a)  = 0 

for  a ■ 0);  and,  a'  and  a"  are  simply  variables  of  integration.  If  complete 

traversal  of  the  drift  space  is  effected  then  x(t,a)  - Ax  for  a = a^.  This 

takes  place  at  time  t ■ t + a^.  For  optoelectronic  sampling,  it  is  only 

necessary  to  establish  whether  or  not  t^  exists,  and  the  duration  of  time  x 

for  which  it  exists,  for  any  given  v(t)  - V , g.  (v  ) and  Ax.  For  the  potential 

differences  allowed  here,  x is  significantly  smaller  than  a.  — wherein  lies 

AX 

another  advantage  of  optoelectronic  sampling  for  high-speed  detection. 

The  fidelity  with  which  the  sampling  process  replicates  f^(t)  can  be 
analyzed  on  the  basis  of  bandwidth  and  phase  relatioishlp  requirements.  For 
simplicity,  we  continue  assuming  that  the  requirement  for  correlating  the 
high-speed  input  signal  with  the  IF  representation  of  the  signal  is  Implicitly 
Included  in  a sequential  sampling  mode  of  operation.  The  sampling  process  can 
be  stated  as  the  product  of  two  Fourier  series;  one  series  representing  the 
optoelectronic  sampling  function  f^Ct),  and  the  other  representing  the  periodic 
Intensity  of  the  input  optical  signal  f^Ct).  Thus,  in  general. 


K 

f.(t)  ■ l cos(ku .t  + <fr.) 

x k-Q  K 1 * 
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E i 


: 


where,  for  a band-limited  signal,  K is  finite.  Similarly,  the  sampling 

function  f (t)  can  be  written  as 

s 

M 

f (t)  *[  b cos(mu)  t + i^)  . (6) 

8 m-0  “ 8 

The  output  signal  fQ(t)  can  be  represented  as  the  product  of  the  sampling 

function  f (t)  and  the  incident  signal  f (t).  The  resulting  product  can  be 

simplified  greatly  by  considering  the  amplifier  as  a low  pass  filter  which 

rejects  high  frequency  components  such  as  ok  and  ug  as  well  as  their  harmonics, 

passing  only  those  terms  containing  the  difference  u . - u : 

X s 

N 

fQ(t)  - C l anbn  cos[n(u.  - ug)t  + 0 ] (7) 

n-0 

where  0Q  ■ $ - ip,  C is  a constant,  N is  equal  to  K or  M,  whichever  is  smaller. 

Then,  choosing  f (t)  such  that  M > K and,  therefore,  N - K and  requiring  that 
bj  - bfl  for  all  j _<  K,  we  see  that  fQ(t)  fully  replicates  f^(t)  at  a greatly 
reduced  frequency  (for  At  « T ) with  a new  amplitude  and  phase  reference: 

K ‘ 

fQ(t)  - CbQ  [ an  cos [n(ui  - w8)t  + 0^]  - Cb^^t)  . (8) 

n-0 

It  Is  interesting  to  note  that  the  requirement  for  equal  harmonic  amplitudes, 
bj  - bo>  applies  only  to  the  terms  that  correspond  to  the  terms  that  exist  in 
the  incident  signal.  Any  additional  components  in  the  sampling  function,  such 
as  bj  4 0 for  j > K,  do  not  serve  any  useful  purpose.  Therefore,  the  ideal 
sampling  waveform  corresponds,  in  the  frequency  domain,  to  a flat-topped  comb 
function,  band-limited  at  Kwg/2Tr  Hz;  the  phases  of  its  harmonics  are  equal  to 
as  in  (6).  This  Insures  not  only  absolute  signal  fidelity  but  also  minimum 
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sampling  power  requirements.  The  importance  of  accurately  characterizing 


fg(t)  for  picosecond  resolution  at  maximum  efficiency  is  evident  from  these 
conclusions. 


III.  CHARACTERIZATION  OF  THE  SAMPLING  FUNCTION 
Whereas  in  conventional  electronic  sampling  the  width  of  the  sampling 
function  is  determined  simply  by  the  Interval  during  which  the  sampling 
bridge  is  forward  biased  by  a "strobe"  pulse,  in  optoelectronic  sampling 
the  sampling  pulse  v(t)  and  other  parameters  affect  the  value  of  r nontrivially. 
In  fact,  due  to  nonzero  initial  velocities  of  photoelectrons,  the  potential 
v(t)  can  affect  the  position  of  electrons  that  are  emitted  well  before 
v(t)-Vr  >_  0.  Moreover,  the  sampling  function  fg(t)  does  not  always  resemble 
v(t),  since  even  if  v(t)  - 0 for  t < 0 it  can  effect  complete  traversal  of 
Ax  by  photoelectrons  emitted  at  t < 0.  The  sampling  function  fg(t)f  the 
sampling  voltage  v(t),  and  the  actual  waveform  of  current  flowing  at  x ■ Ax 
are  all  shifted  chronologically  with  respect  to  each  other  in  the  above  order. 
Again,  a limiting  case  can  serve  as  a convenient  illustration.  For  the  special 
case  of  y ■ 0 and  g (v  ) ■ 6.  (v  - v ),  where  5 'is  the  Dirac  delta  function 

A X a X XO 

(i.e.,  monoenergetic  photoelectrons  corresponding  to  initial  velocity  v^) , 
we  find  that,  regardless  of  the  exact  waveform  of  v(t), 

0 if  x(t,a)  < Ax  for  all  a 

f (t)  - 

S if  x(t,a)  > Ax  for  some  a < T . (9) 

o 

That  is,  if  v(t)  is  sufficiently  large  to  effect  displacements  such 

> 

that  x(t,a)  >_  Ax  for  some  a and  if  v(t)  monotonically  increases  to  a peak 
which  is  followed  by  a monotonic  decrease  then,  remarkably  f (t)  is  a 
rectangular  waveform.’  The  width  x of  fg(t),  but  not  the  amplitude  S^,  is  a 
function  of  the  form  of  v(t).  When  the  initial  velocities  of  photoelectrons 
obey  a relatively  narrow  distribution  function,  which  Is  generally  true  near 


the  photoelectric  threshold,  f (t)  is  still  flat-topped  but  with  rounded 

s 

corners. 

For  the  exact  determination  of  fg(t) — both  its  waveform  and  amplitude — 

specific  values  of  parameters  are  chosen  and  included  in  a numerical  analysis 

prepared  on  the  basis  of  the  considerations  developed  earlier. 

It  should  be  recognized  that,  since  the  displacement  law  is  an  integral 

equation,  the  quantitative  characterization  of  fg(t)  can  be  given  only  for  ’ 

specific  cases.  A computer  analysis  was  formulated  such  that  all 

important  variables  could  be  suitably  incorporated.  The  obvious  parameters 

to  be  scanned  over  chosen  ranges  of  values  are  the  drift  space  dimension  Ax, 

the  retarding  potential  V , and  the  coefficients  specifying  the  accelerating 

potential  v(t) . Moreover,  for  each  optical  wavelength,  each  time  reference 

must  be  incremented  over  the  range  of  initial  velocities  given  by  g (u  ). 

A X 

Choosing  the  cesium-antimony  photocathode  as  an  example,  parameter 

ranges  commensurate  with  picosecond  resolution  and  corresponding  to  driving 

waveforms  described  by  (1)  were  determined.  Standard  assumptions  were  made 

with  regard  to  spatial  uniformity  of  fields  and  negligible  space  charge 

effects.  Only  the  trends  that  best  illustrate  the  concepts  of  optoelectronic 

sampling  and  the  manifestation  of  parameter  values  are  included  in  this  text. 

Figure  4 shows  a typical  result,  demonstrating  the  time  relationship  between 

v(t)  and  f (t)  for  two  values  of  V ; for  this  example,  the  parameter  values 
S IT 

0 

are:  Vq  ■ 14.2V,  Ax  ■ 10pm,  A = 6200A,  y = 0 and  T = 68  psec;  values  for 

g (v  ) were  taken  from  Khlebnikov  and  Melamid^.  It  should  be  noted  that 
ax  > 

the  electronic  waveforms  illustrated  in  Figs.  2 and  4 could  be  instrument 
limited  since  the  best  commercial  electronic  sampling  components  are  designed 
for  20  psec  rise  time  which  is  accompanied  by  several  picoseconds  of  jitter 
in  display;  when  specially  tuned,  the  performance  can  be  somewhat  better. 

Both  rise  time  and  fall  time  (10Z  to  90Z)  of  the  pulse  shown  in  Figure  2 are 
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12 


close  to  Che  sCandard  racings  of  Che  InsCrumenCs.  Ocher  Chan  complete 

speccral  analysis,  which  requires  accurate  phase  measurements  of  all  harmonics, 

no  ocher  sCandards  were  available  to  settle  the  issue.  If  waveforms  v(t) 

could  not  be  represented  analytically  by  the  form  of  (1) , the  time-resolution 

results  for  optoelectronic  sampling  established  in  this  work  would  represent 

only  conservative  estimates  vis-a-vis  the  values  corresponding  to  a v(t)  with 

a higher  time  rate  of  change  than  that  of  sin^y  t . 

The  example  of  f g Ct ) corresponding  to  = 12.3V,  illustrated  in 

Figure  4,  is  shown  to  approach  a peak  value  Sq  which  is  the  theoretical 

mayimuTn  corresponding  to  the  extraction  of  all  photoelectrons  created  by 

f^(t).  Frequently,  in  real  devices,  as  in  the  above  example  but  unlike  in 

the  example  expressed  by  (9),fg(t)  will  not  reach  the  peak  value  S^.  In 

general,  fg(t)  is  not  a rectangular  waveform  and,  in  addition,  for  any  given 

v(t),  its  shape  changes  with  V , Ax,  y, and  g (v  ).  For  these  reasons,  two 

• IT  AX. 

quality  parameters,  p-  and  n_,  designating  sampling  efficiency  and  sampling 

xl  r 

fidelity,  respectively,  were  defined  and  computed  on  the  basis  of  '‘.he 
following  relations: 

- t +x 

nF  - it-  / ° f (O  dt  (10) 

E sor  \ 8 


t +T 


/ 


fs(t)  dt 


nT 


/ 


(n-l)T 


fs(t)  dt 


(ID 


where,  as  illustrated  in  Figure  5,  Sq  is  the  theoretical  maximum  of  f^Ct);  S, 

the  actual  maximum  of  f (t);  and,  t and  t are  such  that  f (t  ) ■ f (t  + t) 

s o so  so 
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Clearly,  the  sampling  efficiency  parameter  nE  quantifies  the  effectiveness 


of  the  sampling  fields  in  extracting  (from  the  drift  space)  photoelectrons 


created  within  the  time  Interval  T.  Thus,  n„.  relates  to  signal  strength.  As 

b 


stated  earlier,  the  need  for  defining  the  sampling  fidelity  parameter  n. 


arises  from  the  fact  that,  for  any  given  sampling  waveform  v(t),  the  analytical 


form  of  f (t)  changes  for  different  parameters.  For  this  reason,  the 


characterization  of  f (t)  by  T only  is  inadequate,  particularly  when  v(t)-V 

s r 


is  small.  It  is  noted  that  if  g,  (u  ) ■ 6,  (v  - v ),  both  r\„  and  Tip  would  be 

AX  AX  XO  £ r 


unity  and  characterization  of  f s (t)  by  t would  suffice.  It  can  be  said  also 


that  the  proximity  of  rip  to  unity  is  a measure  of  the  squareness  of  f (t)  and, 

* s 


therefore,  the  applicability  of  the  familiar  line  spectrum  relationship 


sinx/x»X  * m j/2,  to  the  estimation  of  the  fidelity  of  the  rendition  of  f (t) 
S I 


by  f (t),  (i.e.,  for  the  verification  of  the  condition  b,  = b , j < K) . 
o jo  — 


Retaining  the  same  values  for  and  T as  before  and  continuing  to  keep 


Y ■ 0,  the  effect  of  changing  the  drift  space  dimension  is  shown  in  Figures 


6 and  7 for  two  portions  of  the  visible  spectrum:  X = 6200A  and  X = 4047A, 


respectively.  Although  both  0-  and  n were  computed,  for  clarity,  only  ri 

£ £ £ 


is  shown;  rip  was  determined  to  remain  always  greater  than  0.75. 

i * 
v 


Typically,  the  sampling  efficiency  drops  off  well  before  the  sampling 


fidelity  indicating  that,  under  conditions  of  marginal  signal  level. 


resolution  may  be  limited  to  approximately  5 psec  in  the  red  and  10  psec 


in  the  violet.  (x  decreases  monotonically  with  increasing  wavelength  for 


cesium-antimony  photocathodes.)  However,  under  conditions  of  ample  signal 


levelj picosecond  resolution  sampling,  exhibiting  good  fidelity,  is 


attainable.  Figures  6 and  7 are  also  useful  for  showing  that,  given  any 


desired  x,  no  optimum  drift  space  dimension  exists.  Electron  extraction 


efficiency  monotonically  increases  with  decreasing  Ax.  A sample  curve. 


labeled  7 psec  in  the  space  corresponding  to  n_  values,  has  been  constructed 

£ 


in  Figure  6,  indicating  the  at  1 ■ 6200  X nearly  70Z  extraction  is  possible 

14 


for  7 psec  resolution.  For  T = 4.5  psec,  rig  is  still  ''<0.30  for 

Ax  - 10  pm.  Finally,  the  effect  of  varying  T is  shown  in  Figure  8. 

As  should  be  expected,  it  is  confirmed  that,  corresponding  to  any  given 

efficiency  tu  < 1.0,  T increases  with  T and  the  minimum  T is  attainable 
£• 

with  minimum  T when  V is  kept  constant. 

o 

The  effect  of  a finite  delay  in  photoelectric  emission,  also  referred 
to  as  photocathode  response  time  y,  can  be  significant  in  optoelectronic 
sampling  under  certain  conditions.  Analytically,  a delay  between  the  time 
the  emitting  surface  is  illuminated  and  the  time  electrons  are  ejected  from 
the  surface  is  equivalent  to  a coordinate  translation  such  that  the  dis- 
placement process  described  by  (4)  is  not  initiated  until  a = y.  Letting 
8 ■ a - y,  d8  - da  , we  can  continue  describing  the  displacement  of  electrons 
with  (4)  if  0 is  substituted  for  a everywhere.  Again,  because  it  is  not 
important  when  t occurs,  if  it  occurs,  (other  than  t^y  - t < T as  before) 
an  equal  delay  on  all  photoelectrons  at  any  given  wavelength  has  neither  an 
effect  on  resolution  nor  on  any  other  parameter  previously  computed,  unless 
the  detector  is  simultaneously  monitoring  signals  at  several  wavelengths. 
However,  even  for  the  monochromatic  case,  the  more  general  condition  will 
be  that  corresponding  to  each  velocity  v*  there  will  be  a different  emission 
delay  and,  moreover,  that  these  delays  will  be  described  by  distribution 

functions  in  time.  Letting  y designate  a time  variable  with  origin  at  t of 

• • , 

f . (t) , we  write  the  said  distribution  function  as  G.  (v  ,y)  such  that 
X A X 

> 

8a(Ux)  ° ^ GX^Vx,Y^  dY  * (12) 

Then,  for  any  given  t,  the  contribution  to  f^(t)  expressed  as 

/ / U1U2Cx(ux,y)  dy  dv/x  (13) 
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where  - U^(y  “ Y^)  and  Uj  ■ Uj(y2  ” are  un*t  step  functions  such  that 
Y2  L and 

0 for  y < Yi 

U1  " 

1 for  y > Y-. 


0 for  y > Y2 

1 for  y lYj 


The  values  for  y^  and  y2  are  found  by  establishing  the  range  of  delays, 

±y  — y2’  ^or  cax  ex*st*  Equation  (13)  is  integrated  over  the 

ranges  y.  < y < y-  and  0 < v < v where  v is  the  maximum  velocity 

in  g,(u  )•  Equation  (4)  is  solved  for  each  velocity.  The  solution  requires 
the  use  of  the  coordinate  translations  8 = a - y,  analogous  to  the  translation 
corresponding  to  the  simple  case  of  a single,  non-zero  delay  y.  Thus,  in  the 
general  case  of  emission,  the  analytical  treatment  in  the  evaluation  of  the 
displacement  of  photoelectrons  for  each  given  and  y can  be  reduced  to  the 
simple  case  of  a single,  non-zero  delay.  And,  although  it  was  noted  that  a 
single  delay  on  all  v did  not  influence  any  sampling  parameters,  the  influence 
of  emission  delays  described  by  G. (v  ,y)  will  be  generally  manifested  on  the 

A X 

sampling  function  f (t)  by  virtue  of  (13).  For  example,  the  following 

s 

general  observation  can  be  made.  When  the  high  energy  side  of  g (v  ) is  delayed 

A X 

with  respect  to  the  low  energy  portion,  the  lagging  edge  of  f (t)  becomes 
more  abrupt  and  the  leading  edge  rises  more  progressively  than  for  y = 0; 
the  converse  is  true  when  low  energy  photoelectrons  are  delayed  with  respect 
to  those  described  by  the  high  energy  portion  of  g (v  ) — which  is  the 
expected  behavior  for  real  materials.  For  illustration  and  for  an  appraisal 
of  the  magnitude  of  the  influence  of  delayed  emission  on  the  sampling 
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The  delay  In  emission  corresponding  to  ■ 0 is  and,  as  before,  the  maximum 
velocity  in  g.  (v  ) is  v . Although  simple  in  form,  this  representation  of 

AX  X »Q<1X 

emission  describes  a real  material  better  than  a single  average  delay  by  spe- 
cifying an  average  delay,  Y(ux)»  for  each  velocity  vx«  The  effect  of  emission 
delays  on  the  characteristics  of  f (t)  for  this  case  was  investigated  on  the 
basis  of  (13)- (16)  by  solving  (4)  in  accordance  with  the  coordinate  translations 
B(V^)  ” a - Y(ux)*  Figure  9 summarizes  the  results  of  the  study  corresponding 

to  v(t)  of  (1),  g, (v  ) of  cesium-antimony,  and  the  following  parameters: 

A X 

0 < < 54  psec,  r * 68  psec,  V **  14.2  V,  Ax  * 10  Urn,  and  A *»  6200A. 

o o 

For  this  form  of  G. (v  ,Y),  even  small  values  of  emission  delays  markedly  affect 
AX 

the  efficiency  of  sampling  without  significantly  altering  resolution.  This 

suggests  that  the  sampling  function  broadens  at  the  base  without  commensurately 

. 

increasing  at  half-amplitude.  The  value  of  the  integral  corresponding  to  the 

denominator  of  (11)  retains  approximately  the  same  value.  At  higher  values 

of  Y0»  both  the  base  of  fg(t)  and  t Increase  but  still  disproportionately  as 

> 

evidenced,  for  example,  by  only  a 10%  increase  in  T for  a 50%  loss  in  rig  when 
Y is  as  large  as  54  psec.  In  addition  to  changes  in  n„  and  T,  emission  delays 

O L 

also  affect  the  waveform  of  f (t).  Of  course,  only  complete  frequency  spectra 

s 

of  fg(t)  with  detailed  phase  information  corresponding  to  (6)  would  determine 
accurately  the  effect  of  waveform  change  on  sampling  resolution. 
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IV.  CONCLUSIONS 

Th«  formalism  required  for  the  analysis  of  optoelectronic  sampling 
is  fully  developed  on  the  basis  of  the  dynamics  of  electrons  produced  as 
a result  of  photoelectric  interaction.  It  is  shown  that  the  sampling  of 
such  electrons  at  the  photocathode  surface  can  be  effected  with  high 
efficiency  within  time  intervals  as  short  as  a few  picoseconds  provided 
certain  conditions  are  satisfied.  For  maximum  time  resolution  at 
maximum  sampling  efficiency,  the  general  guidelines  deduced  from  the 
analysis  of  optoelectronic  sampling  indicate  the  need  for  ultrashort 
sampling  voltage  waveforms,  small  dimensions,  thin  photocathodes — 
preferably  designed  for  front-illumination,  and  operation  at  wavelengths 
near  the  photoelectric  threshold.  Of  course,  only  special  geometries, 
such  as  microstrip  configurations,  are  compatible  with  the  propagation 
of  such  ultrashort  voltage  waveforms.  These  qualitative  observations  are 
supported  by  numerical  examples  which  illustrate  the  potential  performance 
of  optoelectronic  samplers  both  under  ideal  conditions  and  under  limita- 
tions existing  with  real  materials.  In  particular,  characteristics 
describing  optoelectronic  sampling  are  fully  specified  for  some  select 
dimensions,  material  parameters,  and  driving  waveforms.  Sampling  widths 
of  a few  picoseconds  are  obtained  for  several  combinations  of  parameter 
values.  Given  that  sufficient  intensity  exists  to  overcome  considerable 
losses  in  sampling  efficiency,  single  picosecond  structures  of  optical 
signals  can  be  resolved  provided  geometrical  constraints  are  observed. 

It  should  be  noted  that  small  fractions  of  a millimeter  in  path  length 
difference  are  important  in  picosecond  detection,  both  when  considering 
the  extent  of  the  optical  source  and  the  extent  of  the  photocathode. 

The  possibility  of  implementing  relatively  large  area  photocathodes 

in  high-speed  optoelectronic  sampling  devices  is  only  implicitly  included 
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in  this  text  with  the  statement  of  the  condition  Az  « Tc.  For  example, 

whereas  a commercial  high-speed  crossed  field  photodetector  is  limited  in 

2 16 

photocathode  area  to  4 mm  and  in  risetime  to  ^ ISO  psec  , a detector  equiv- 
alent in  speed  but  designed  on  the  optoelectronic  sampling  principle  could  have 
a photoactive  area  ten  times  that  value.  Geometrical  comparisons  are  not 
possible  for  detectors  suitable  for  picosecond  resolution  since  no  fully 
electronic  devices  exist  to  the  best  knowledge  of  the  authors.  However,  it 
is  clear  that  even  for  10  psec  resolution,  geometrical  considerations  will 
play  an  important  role. 

Finally,  it  is  important  to  emphasize  that  the  dynamics  of  photo- 
electric emission  is  inadequately  documented  at  present  for  the  purpose 
of  accurate  optoelectronic  sampling  design.  A complete  specification  of  emis- 
sion parameters,  expressed  here  by  G. (v  ,Y)»  is  required  for  each  material 
for  a complete  quantitative  analysis  of  the  sampling  of  photoelectrons. 
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FIGURE  CAPTIONS 
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. Illustration  of  the  sampling  geometry. 

. Superposition  of  the  waveform  expressed  by  (1),  for  Vq  ° 14.2V 

and  r - 68  psec,  on  an  oscilloscope  trace  of  a driving  potential 

v(t,z  ).  Vertical  scale:  2V/div;  horizontal  scale:  20  psec/div. 
o 

. Illustration  of  waveforms  encountered  in  optoelectronic  sampling. 

A.  Waveform  of  the  envelope  of  optical  radiation  incident  on  the 

detector.  B.  The  sampling  potentials.  C.  The  sampling  function 

effected  by  v(t)  - V . D.  The  sampler  output  shown  on  a time 
T 

axis  scaled  by 

. Example  of  the  time  relationship  between  v(t)  and  ffi(t)  shown  on 
10  psec/div  horizontal  scale  for  two  values  of  V . The  values 
of  T are  9.6  and  4.9  psec  for  V^.  of  12.3  and  12.7V,  respectively. 

Other  parameters  are  Vq  ■ 14.2V,  Ax  ■ 10pm,  y - 0,  T = 68  psec; 

O 

g^(v^)  corresponds  to  cesium-antimony  at  X = 6200A. 

. Graphical  representation  of  quantities  defining  the  quality 
parameters  Tig  and  rip. 

Dependence  of  T and  q on  V at  X » 6200&,  for  five  values  of  drift 

space  dimension  Ax^  2 3 ^ 5 ■ 10,20,30,40,  and  50  pm,  respectively. 

The  monotonic  increase  of  q_  with  decreasing  Ax  is  shown  by  the  constant 

£• 

T line,  labeled  7 psec. 

. Dependence  of  T and  qc  on  V at  X * 4047a,  for  five  values  of  drift 

b r 

space  dimension  Ax^  2 3 4 5 **  ^-0,20,30,40,  and  50  pm,  respectively. 
Dependence  of  T and  q_,  on  V for  r,  " 100,80,  and  60  psec, 
respectively. 

i.  Illustration  of  the  effect  of  emission  delays,  described  by 

(15)  and  (16),  on  T and  q£  for  YQ  1 2 3 4 " 0,18,36,  and  54  psec. 


respectively. 
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Figure  9. 


Adaptation  of  Optoelectronic  Sampling  to  the 
Study  of  the  Dynamics  of  Photoemission 


Spectral  response,  quantum  yield,  and  the  electron  velocity  distri- 
bution of  photoelectric  emission  have  been  the  standard  parameters  used 
in  interpreting  the  mechanisms  of  photoemission.  Even  in  the  study  of 
Garbe*’  on  the  factors  affecting  photoemission  from  cesium  oxide  covered 
GaAs,  a negative-electron-affinity  material,  no  consideration  is  given 
to  the  time  evolution  of  photoemission;  average  escape  depths  are  com- 
puted on  the  basis  of  quantum  yield  and  optical  absorption  coefficients. 
We  can  estimate  only  indirectly  (and  with  a large  margin  of  error)  the 
time  response  of  such  materials.  Clearly,  the  considerations  discussed 
at  the  end  of  Section  B-V  and  the  questions  they  raise  are  not  resolved 
by  simple  escape  depth  estimates. 

The  tools  for  studying  the  entire  evolution  of  photoemission  on 
picosecond  time  scales  can  now  be  designed  and  fabricated.  We  propose 
to  use  our  existing,  cw-operated,  mode-locked  dye  laser  for  the  source 
of  tunable,  picosecond  pulses.  The  time  resolved  detection  of  photoelec- 
trons will  be  carried  out  with  an  apparatus  based  on  an  adaptation  of  the 
high-speed  sampler  described  in  Section  B-V.  Since  all  studies  on  photo- 
emission require  electron  energy  separation,  the  high-speed  sampler  must, 
fortuitously,  operate  under  idealized  conditions,  i.e.,  if  the  photoelec- 
trons are  passed  through  an  energy  filter  before  entering  the  drift  space 
D of  the  sampler,  the  analysis  of  the  sampling  process  is  reduced  to  the 
limiting  case  — given  as  example  on  p.  14,  B-V  — corresponding  to 
g. (V  ) - 5.(v  - V ).  Under  these  conditions  of  monoenergetic  photoelec- 

trons,  the  resolution  of  the  sampler  is  highest.  Picosecond  sampling 
should  be  achievable  for  materials  with  a good  quantum  yield. 
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A schematic  of  the  high-speed,  microstrip  sampler  adapted  to  a 
velocity-selective,  time-evolution  measuring  apparatus  Is  shown  In 


Figure  .1.  A thin  strip  photocmittcr  E is  placed  at  a distance  from  the 
microstrip  sampler  such  that  photoelectrons  arc  first  accelerated  into  a 
magnetic  field  B oriented  perpendicularly  with  respect  to  their  velocity. 
For  any  one  value  of  B,  only  one  velocity  group  of  photoelectrons  will  be 
incident  on  the  strip  of  the  sampling  configuration.  This  strip  is 
now  designed  as  a meshed  slit.  All  other  components  and  configurations 
are  the  same  including  the  conventional  electron  multiplier  for  the 


amplification  of  the  sampled  current.  The  analysis  of  Section  B-V 

simplified  since  g (v  ) is  always  6 (v  - v ).  If  within  any  period  n, 

f^(t)  cC.6(t  - tn),  then,  for  any  given  magnetic  field  B ^ , j = 1,2,..., 

the  apparatus  measures  G, (v  ,y)  „ . • A magnetic  field  scan  yields 

A xu  = v (B, ) 

X x j7 

the  entire  ^(v^.y). 

Since  no  such  measurements  have  been  made  in  the  past,  there  appears 
to  be  no  need  at  this  time  to  reconsider  our  understanding  of  the  time 
evolution  of  photoemission.  However,  it  is  anticipated  that  the  processes 
controlling  the  evolution  of  photoemission,  diffusion  and  tunneling,  will 
be  manifested  quantitatively  in  these  time  studies,  allowing  for  revision. 
Adding  to  the  basic  understanding  of  electron  emission  from  surfaces,  the 
results  of  these  efforts  should  also  contribute  to  our  ability  to  characterize 
photoactive  materials  for  high-speed  optoelectronic  applications. 


S.  Garbe,  "Factors  Affecting  the  Photoemission  from  Caesium  Oxide  Covered 
GaAs,"  Solid-State  Electronics,  vol.  12,  pp.  893-901,  1969. 


Fig.  1-  Illustration  of  the  geometry  of  an  electron  sampler  including 
a photoclectron  energy  separating  provision:  li,  photocm Liter; 
D,  drift  space;  Mj,  meshed  strip;  M^,  meshed  ground  plane; 

Ax,  drift  space  dimension. 
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[57]  ABSTRACT 

An  electronic  sampling  head  comprising  a housing 
having  a path  for  impinging  radiation,  a photoemissive 
element  positioned  across  said  path,  photoelcctron 
ejection  means  positioned  in  proximity  to  said  photo- 
emissive  element,  means  for  applying  ultrashort  sam- 
pling potential  differences  between  said  photoemissive 
element  and  said  photoelcctron  ejection  means  for  ac- 
celerating photoelcctrons  emitted  by  said  photoemis- 
sive element  in  response  to  said  incident  radiation, 
and  means  for  detecting  said  photoelectrons,  all  result- 
ing in  the  capability  of  detecting  high-speed  radiation 
phenomena  with  waveform  fidelity. 

17  Claims,  9 Drawing  Figures 
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Figure  8.  Sampling  function  width,  T,  and  sampling  efficiencies 
and  nA  as  a function  of  mlcrostrlp  spacing  for  sam- 
pling pulses  of  10 V peak,  100  psec  basewldth. 
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I 

OPTOELECTRONIC  SAMPLING  HEAD 

The  invention  herein  described  was  made  under  a 
contract  with  the  Department  of  the  Navy. 

This  invention  relates  to  optical  detection  devices 
and  specifically  to  optical  sampling  heads  usable  in 
conjunction  with  photoclectron  amplifying  devices 
suitable  for  low  power  radiation  detection 

A need  exists  for  a high  resolution  sampling  detector 
for  detecting  radiation  levels,  particularly  optical,  at 
relatively  high  sampling  speeds.  Since  in  many  applica- 
tions the  photon  flux  level  at  high  sampling  speed  will 
be  low,  a sensitive  device  for  detection  is  required. 

A conventional  form  of  detection  device  is  the  photo- 
multiplier. In  its  usual  form,  the  single  stage  dev-  e 
employs  a photocathodc,  emitting  electrons  in  re- 
sponse to  impinging  radiation,  a secondary  emissive 
electrode  or  dynode,  and  a collector  electrode  or  an- 
ode. In  a multistage  electron  multiplier,  a cascade  of 
dynodes  arc  employed,  with  successively  higher  biasing 
potentials,  such  that  impinging  primary  electrons  cause 
emission  of  secondary  electrons  in  a multiplying  pro- 
gression from  successive  dynodes,  and  are  ultimately 
collected  by  the  anode  which  in  turn  provides  an  out- 
put signal. 

It  is  desirable  to  utilize  an  electron  multiplier  for 
providing  high  resolution,  high  speed  sampling  of  opti- 
cal radiation;  however,  conventional  photocathodic 
detecting  structures  do  not  possess  adequate  resolution 
nor  speed  to  enable  rapid  sampling  of  optical  signals, 
particularly  in  the  measuring  or  detecting  of  a rapidly 
varying  light  level.  Such  light  levels,  which  may  be  laser 
generated,  may  require  sampling  rates  approaching  the 
picosecond  range.  Additionally,  the  use  of  high  sam- 
pling rates  is  best  controlled  electronically,  at  the  ac- 
tual sampling  point.  This  latter  feature  is  difficult  to 
achieve  in  conventionally  gridded  photomultiplier 
structures. 

It  is  therefore  the  primary  object  of  the  present  in- 
vention to  provide  a sampling  device  which  is  operable 
with  high  resolution  over  a short  sampling  time. 

It  is  another  object  of  the  present  invention  to  pro- 
vide a sampling  device  adapted  to  be  employed  with  a 
conventional  electron  multiplier  for  improving  the  res- 
olution thereof. 

It  is  a further  object  of  the  present  invention  to  pro- 
vide a sampling  device  which  may  be  electronically 
controlled  to  produce  sampling  intervals  of  shorter 
duration  than  heretofore  possible  while  maintaining 
high  resolving  characteristics. 

The  foregoing  objects  are  realized  in  a device  con- 
structed of  a housing  having  a path  for  impinging  radia- 
tion and  a photoemissive  element  positioned  across 
said  path.  A photoclectron  ejection  electrode  is  posi- 
tioned adjacent  the  photoemissive  element  and  a sam- 
pling potential  difference  can  be  applied  across  the 
element  and  electrode.  Any  electrons  emitted  by  the 
element  due  to  the  incident  radiation  thereon  during 
the  sampling  pulse  duration  is  accelerated  through  the 
ejection  electrode  to  a suitable  detection  device  such  as 
a conventional  electron  multiplier.  The  device  struc- 
ture includes  an  elongated  photoemissive  strip  located 
in  close  proximity  to  the  ejection  electrode,  the  latter 
taking  the  form  of  a circular  wire  mesh.  The  width  of 
the  strip  and  its  distance  from  the  mesh  are  set  with 
predetermined  relationship  for  optimum,  high  speed 
results. 


2 

The  foregoing  objects  and  brief  description  of  the 
present  invention  will  become  more  apparent  from  the 
following  detailed  description  and  appended  drawings 
wherein 

5 FIG.  1 shows  the  head  of  the  present  invention  uti- 
lized with  an  electron  multiplier; 

FIG.  2 is  a cross  section  along  the  optical  path  of  the 
head  showing  the  housing  and  photoemissive  strip; 

FIG.  3 illustrates  the  insert  and  ground  plane  cle- 
10  ment; 

FIG.  4 shows  the  insert  in  position  in  the  housing; 

FIG.  5 is  a top  view  showing  the  relationship  of  the 
strip  to  the  underlying  ground  plane  forming  a micro- 
strip  configuration. 

13  FIG.  6 is  an  illustration  of  the  circuit  utilizing  the 
head  of  the  present  invention;  and, 

FIGS.  7,  8 and  9 illustrate  as  examples  the  graphical 
relationships  of  the  effect  various  dimensional  parame- 
ters and  biasing  have  on  resolution  and  efficiency. 

20  Referring  to  FIG.  1,  sampling  head  10  constructed  in 
accordance  with  the  present  invention  is  shown  in 
mounted  position  on  a conventional  electron  multiplier 
12.  For  purposes  of  illustration,  the  electron  multiplier 
may  be  that  of  an  RCA  Model  7102  photo-tube  includ- 
25  ing  a dynode/anodc  structure  14  including  a plurality  of 
voltage  leads  16  for  anplying  a potential  to  the 
dynode/anode  electrode  structure  14.  Other  conven- 
tional electron  multipliers  may  also  be  employed  with 
the  sampling  head.  The  photo-tube  12  has  been  modi- 
30  fied  to  the  extent  that  an  external  source  of  electrons, 
supplied  by  the  sampling  head  10,  may  be  provided.  To 
this  end,  a flange  18  including  mounting  bolts  20  pro- 
vides a support  frame  for  mounting  the  head  10  to  the 
tube  12.  An  O-ring  seal  22  may  be  provided  to  vacuum 
35  and  optically  seal  the  head  outlet  to  the  tube  inlet. 

On  the  left  side  of  the  head  10,  referring  to  FIG.  1,  a 
lens  or  window  unit  24  is  bolted  to  the  optical  inlet  side 
of  the  head  10.  A further  O-ring  seal  26  may  be  used  to 
vacuum  seal  the  optical  element  through  the  head.  The 
40  head  is  bolted  by  means  of  bolts  28  to  the  lens  structure 
24. 

incident  radiation  30  impinges,  through  the  lens  24, 
upon  the  photosensitive  surface  creating  photoelec- 
trons which  are  sampled  by  the  head  10  in  a manner  to 
45  be  set  forth  in  further  detail  below,  and  the  resulting 
electron  flow  multiplied  in  the  electrode  structure  14  of 
the  electron  multiplier  12. 

The  structure  of  the  sampling  head  10  is  shown  in 
greater  detail  in  F'lGS.  2-5.  The  cross  section  FIG  2 
30  shows  that  the  sampling  head  includes  a centra!  optical 
path  32  passing  through  a housing  indicated  generally 
as  34.  The  upper  flange  36  of  the  head  structure  34 
includes  bolt  holes  38  for  accommodating  the  bolts  28, 
and  a groove  40  for  accommodating  the  O-ring  26, 
53  shown  in  FIG.  I.  The  lower  flange  42  of  the  head  struc- 
ture 34  includes  bolt  holes  44  for  accommodating  the 
bolts  20,  and  a groove  46  for  accommodating  the  ta- 
ring 22. 

t he  housing  34  further  includes  a side  wall  48  which 
W1  may  be  square  or  circular  in  shape.  The  internal  pas- 
sage way  32,  axially  defining  the  optical  path,  is  prefer- 
ably round  in  shape  and  include,  a plurality  of  threads 
50  formed  in  the  interior  wall  for  accommodating  an 
insert,  shown  in  FIG.  3.  The  housing  34,  including 
*>3  flanges  and  side  walls,  as  preferably  constructed  of  a 
cold  rolled  steel  or  like  sturdy  conductive  material. 

Incorporated  into  the  side  walls  48  are  a pair  of  feed- 
through leads  52  which  include  a conductive  cap  54 
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having  provision  for  attaching  or  incorporating  a con- 
ductive wire  56.  The  cap  54  is  supported  on  the  hous- 
ing by  a support  flange  58.  The  feedthrough  is  formed 
as  a plug  fitting  into  the  wall  48,  and  includes  an  outer 
conductive  shell  60  and  an  inner  glass  feedthrough  5 
member  62.  Through  the  feedthrough  is  connected  a 
conductive  photoemissive  strip  64,  composed  of  a suit- 
able photocmissivc  material  such  as  cesium  antimony 
or  the  like.  The  strip  may  be  finely  perforated  or  semi- 
transparent. or  formed  as  a mesh,  and  has  tapered  ends.  10 
The  component  material  structure  of  the  feedthrough 
is  designed  to  permit  the  coupling  of  the  broadbound 
high  frequency  sampling  signals  to  the  strip  64. 

The  strip  64  is  positioned  across  the  optical  path  32. 

The  insert  66,  FIG.  3,  is  a conductive,  externally  15 
threaded,  round  plug-like  component  including  a cen- 
tral electron  path  68  and  a groove  70  for  accommodat- 
ing a spacer  72.  The  upper  part  of  the  insert  66  includes 
a fine  wire  mesh  74  which  can  be  circular  and  which  is 
affixed  to  the  insert  66  as  by  spot  welding  or  the  like. 
The  wire  mesh  74  performs  the  function  of  a photoe- 
lectron ejection  electrode  in  cooperation  with  the  pho- 
toemissive strip  64,  and  requires  a certain  spacing  with 
respect  to  the  strip  64  for  electronic  matching  in  a 25 
microstrip  configuration  The  insert  66  is  positioned 
into  the  body  of  the  structure  34  by  rotating  the  mating 
threaded  portions  of  the  insert  66  and  the  interior  50  of 
the  housing  34  until  the  spacer  72  contacts  the  lower 
projecting  portion  of  the  shell  60,  shown  in  FIG.  4.  30 
Accurate  control  of  the  spacer  dimension  will  control 
the  spacing  distance  d,  between  the  strip  64  and  the 
microstrip  ground  plane  w hich  acts  as  an  ejection  elec- 
trode 74.  The  electrode  74  is  termed  an  ejection  elec- 
trode in  that  it  performs  the  function  of  causing  elec-  35 
trons  from  the  strip  64  to  be  ejected  into  the  path  32. 

The  top  of  the  assembly,  viewed  in  the  direction  of 
the  optical  path  32,  is  shown  in  FIG.  5.  The  strip  64  is 
designed  to  have  a predetermined  width,  d2,  across  the 
optical  path,  which  also  provides  maximum  operational  40 
characteristics. 

By  way  of  example,  the  wire  mesh  may  be  con- 
structed of  a woven  or  etched  nickel  screen.  Other 
conductive  metallic  materials  compatible  with  photo- 
tube processing  may  also  be  employed.  45 

Due  to  the  high  frequency  of  operation,  the  compo- 
nents forming  the  housing,  insert,  strip  and  feedthrough 
are,  as  shown  in  FIGS.  2-5,  designed  to  be  assembled 
employing  techniques  which  will  result  in  a broadband, 
nondispersivc  electrical  system  suitable  for  high  vac-  50 
uum  construction.  The  relationship  of  the  dimensions 
</,  and  d2  arc  critical  with  respect  to  the  device  opera- 
tion. These  parameters  are  mutually  interdependent  in 
determining  the  mode  propagation  in  accordance  with 
broadband  design  criteria,  and  should  be  scaled  to  55 
produce  a reflcctionless  and  dispcrsionlcss  fundamen- 
tal mode  propagation.  An  example  of  such  dimension 
selection  will  be  found  in  the  IRE  Transactions  of  the 
Professional  Group  on  Microwave  Theory  and  Tech- 
nique entitled,  "Symposium  on  Microwave  Strip  Cir-  60 
cuits",  March,  1955. 

According  to  the  embodiment  disclosed  herein, 
where,  as  illustrated  in  FIGS.  7-9,  resolution  of  optical 
pulse  durations  in  the  picosecond  (10  '■  seconds) 
range  can  be  achieved.  I he  spacing  d,  was  chosen  at  (,5 
300  microns  to  which  corresponds  a strip  width,  of 
1.5  millimeters,  all  designed  in  accordance  with  micro- 
strip design  criteria  and  matched  for  a 50  ohm  termma- 
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tion.  The  optical  path  68,  by  way  of  reference,  is  about 
1.5  inches  in  diameter. 

The  electron  multiplier,  which  serves  as  the  integrat- 
ing-amplifier  section,  along  with  the  output  circuitry,  is 
mounted  in  a large  glass  to  metal  seal.  This  is  then 
brazed  to  an  O-ring  flange  to  allow  the  sampler  portion 
to  be  bolted  to  the  multiplier  section.  Before  the  two 
sections  are  joined  together,  antimony  is  deposited  on 
the  silver  coated  stainless  steel  strip  by  a vacuum- 
evaporation  process.  The  strip  is  soldered  to  the  center 
pins  of  the  feedthrough  connectors,  and  the  sampling 
portion  is  bolted  to  the  integrating  amplifier  section. 
Next,  cesium  gas  is  introduced  into  the  previously  evac- 
uated system  while  the  tube  is  baked  in  an  oven  at  1 50° 
C to  form  a cesium  antimonide  photocathode  (Cs3Sb). 

Operation  of  the  sampling  head  is  set  forth  in  FIG  6. 
The  head.  10,  in  position  on  the  electron  multiplier  12, 
is  biased  by  applying  a retarding  potential  V0  to  the 
strip  64.  This  potential  V„  is  applied  through  a conven- 
tional bias  tee  76,  available  from  the  General  Radio 
Corp.,  for  keeping  the  d.c.  level  free  of  sampling  pulses. 
The  potential  V0  is  supplied  to  the  electrode  64  by 
contacting  lead  78  to  feedthrough  lead  56A.  thereby 
rendering  the  potential  of  the  electrode  74  negative 
with  respect  to  the  strip  64.  An  optical  signal  30  is 
directed  to  the  sampling  head,  as  described  above.  A 
sampling  pulse  is  supplied  by  the  sampling  pulse  gener- 
ator 80  along  the  line  82  through  the  feedthrough  lead 
56B  to  the  strip  64.  The  sampling  pulse  — Vs,  applied  to 
strip  64,  is  negative  with  respect  to  the  ground  elec- 
trode 74  and  serves  to  both  cancel  and  exceed  in  mag- 
nitude the  retarding  potential  V0.  Should  radiation  be 
incident  upon  the  strip  during  the  sampling  pulse  inter- 
val, electrons  emitted  from  the  strip  64  in  response  to 
incident  photons  will  be  accelerated  by  the  potential 
difference  toward  the  photoclectron  ejecting  element 
(the  ground  plane  wire  mesh  74)  and  pass  through  the 
electron  path  68  into  the  electron  multiplier  12.  jince 
the  electron  multiplier,  under  the  bias  -t-V„,  serves  to 
amplify  the  electron  flow,  an  output  signal  Vs  will  be 
provided  along  line  84. 

By  way  of  example,  and  with  the  spatial  dimensioning 
set  forth  above,  the  sampling  head  was  employed  with 
an  optical  pulse  of  about  300  pscc  rise  time  derived 
from  a mode-locked  HeNe  laser.  Using  a —7  volt,  250 
psec  wide  sampling  pulse,  with  Vo  set  at  —2.2  volts,  full 
time  resolution  of  the  mode  locked  laser  pulses  were 
achieved. 

Referring  to  FIGS.  7,  8 and  9 graphical  relationships 
showing  the  criticality  of  the  dimensional  relationships, 
as  well  as  the  effect  of  varying  potential  levels,  and 
sampling  pulse  widthes  are  demonstrated. 

The  parameters  found  useful  in  characterizing  the 
optoelectronic  sampler  are  the  sampling  efficiencies 
t}t  , rjB  and  the  pulse  width  r.  Thus.  tjt  represents 
the  ratio  of  the  relative  number  of  photoelcctrons  of 
the  derived  sampling  function  to  the  maximum  possible 
number,  A„t.  of  photoelectrons  that  could  be  sampled 
in  the  same  interval  for  the  same  set  of  conditions,  and 
r)t  represents  the  ratio  of  the  sampled  photoelectrons 
of  the  sampling  function  to  the  total  number  of  sampled 
photoelcctrons  in  the  entire  sampling  function  Clearly, 
t)t  and  7] , are  maximally  unity  and  are  specified  to 
designate  amplitude  sampling  ctficicncy  and  sampling 
quality  respectively. 

As  a result  of  this  analysis,  changes  in  tjt  . tj,  and  r 
due  to  variations  in  the  input  conditions  can  be  used  to 
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optimize  the  output  signal  with  respect  to  resolution 
and  signal  magnitude. 

The  retarding  potential  represents  the  constant  volt- 
age applied  to  the  sampler  structure  which  prevents 
photoelcctrons  from  reaching  the  sampler  anode  at  all 
times  that  the  sampling  voltage  pulse  is  not  present. 
When  the  sampling  voltage  pulse  does  occur  with  suffi- 
cient amplitude,  the  net  voltage  difference  between  the 
photocathode  and  anode  can  result  in  photoelcctrons 
traversing  the  sampler  drift  space.  As  the  retarding 
potential  is  increased,  the  effective  time  during  which 
the  photocathode  is  forward  biased  decreases,  resulting 
in  improved  resolution,  as  shown  in  FIG.  7. 

Although  improved  resolution  is  desirable,  one  must 
consider  the  effect  of  increased  retarding  potential  on 
sampling  efficiencies.  tjt  and  tj  , tend  to  fall  when  the 
retarding  potential  is  increased  to  a potential  for  which 
the  effective  duration  of  the  sampling  pulse  is  too  short 
to  allow  those  photoelcctrons  with  comparatively  small 
initial  velocities  to  reach  the  anode. 

FIG.  8 illustrates  the  effects  of  different  spacings,  d,, 
between  the  photocathode  and  the  sampler  anode  for  a 
fixed  set  of  parameters.  From  these  numerical  results,  it 
can  be  shown  that  if  the  duration  of  the  accelerating 
field  is  approximately  equal  to  or  smaller  than  the  time 
necessary  for  all  the  emitted  photoelcctrons  to  traverse 
the  drift  space  then  only  those  photoclectrons.  if  any, 
with  large  initial  energies  will  be  sampled.  This  results 
in  a decrease  in  time  dispersion  due  to  the  initial  veloc- 
ity spread  of  the  photoelectrons  as  is  illustrated  in  FIG. 
8 for  d,  greater  than  70  microns.  Also  note  that  the 
sampling  efficiency  tjt  rapidly  deteriorates  as  fewer 
photoelcctrons  arc  collected. 

FIG.  9 illustrates  the  dependence  of  the  sampling 
function  width  and  the  sampling  efficiency  on  the  inci- 
dent wavelength  for  a cesium  antimonidc  photocath- 
ode. It  should  be  noted  that,  in  all  cases  considered, 
wherever  the  sampling  efficiencies  dropped,  the  order 
i)A  > rjr  was  preserved,  indicating  that  though  signal 
magnitudes  could  become  small,  signal  fidelity  would 
be  retained.  For  the  exemplary  cesium  antimony  pho- 
tocathode considered,  it  appears  that  the  highest  reso- 
lution and  greatest  sampling  efficiencies  occur  in  the 
red  part  of  the  visible  spectrum.  Those  incident  wave- 
lengths which  result  in  the  narrowest  range  of  energies 
for  the  emitted  photoclectrons  will  be  capable  of  the 
highest  resolution  since  the  transit  time  dispersion  ef- 
fect which  degrades  the  time  response  to  a sampling 
photomultiplier  tube  is  due  to  the  spread  of  initial  ve- 
locities of  the  emitted  photoclectrons. 

Although  the  invention  has  been  described  in  con- 
junction with  certain  preferred  embodiments,  it  will  be 
understood  that  other  variations,  modifications,  substi- 
tutions, additions  and  deletions  may  be  made  within  the 
framework  and  spirit  of  the  invention. 

What  is  claimed  is: 

1.  A high  speed  sampling  head  for  photoclectrons 
comprising  a housing  having  a path  for  impinging  radia- 
tion, a strip-shaped  photoemissive  element  positioned 
across  said  path,  planar  photoelectron  ejection  means 
positioned  in  proximity  to  said  photo-emissive  element 
in  said  path,  means  for  propagating  a high  speed  sam- 
pling potential  difference  between  said  photoemissive 
element  and  said  photoelectron  ejection  means,  means 
for  accelerating  photoelec  Irons  emitted  by  said  photo- 
emissive  element  in  response  to  said  incident  radiation, 
and  means  for  detecting  said  photoclectrons,  whereby 
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said  ejection  means  constitute  a microstrip  ground 
plane  for  said  photoemissive  element. 

2.  The  head  of  claim  1 wherein  said  photoclcctron 
ejection  means  is  a circular  wire  mesh. 

3.  The  head  of  claim  I wherein  said  photoemissive 
element  is  an  elongated  shaped  strip  across  said  optical 
path  having  a predetermined  narrow  dimension  width 
relative  to  the  spacing  from  said  photoclcctron  ejection 
means. 

4.  The  head  of  claim  3 wherein  said  photoemissive 
element  is  perforated. 

5.  The  head  of  claim  1 w herein  said  means  for  detect- 
ing is  an  electron  multiplier. 

6.  rhe  head  of  claim  I wherein  a first  bias  potential 
is  applied  to  said  photoclcctron  ejection  means 
whereby  said  photoclcctron  ejection  means  becomes 
negative  with  respect  to  said  photoemissive  element 
and  wherein  said  sampling  potential  has  a magnitude 
greater  than  said  first  bias  potential,  thereby  rendering 
said  photoclectron  ejection  means  momentarily  posi- 
tive with  respect  to  said  photoemissive  element 

7.  An  optoelectronic  sampling  head  comprising  a 
housing  having  an  optical  path  for  impinging  light,  a 
photoemissive  elongated  strip  positioned  across  said 
optical  path  and  supported  in  said  housing  by  broad- 
band insulated  electrical  feedthrough  elements,  an 
insert  positionablc  in  said  housing,  said  insert  including 
a photoelectron  ejection  ground  plane  and  an  electron 
path,  and  spacer  means  mounted  to  said  insert  and 
positioned  with  said  insert  against  said  housing  for 
positioning  said  photoelectron  ejection  ground  plane 
with  predetermined  spatial  distance  from  said  photo- 
emissive strip. 

8.  The  head  of  claim  7 wherein  said  photoelcctron 
ejection  ground  plane  is  a circular  wire  mesh. 

9.  The  head  of  claim  8 wherein  said  photoemissive 
clement  is  perforated  or  semitransparent. 

10.  The  optoelectronic  sampling  head  of  claim  7 
further  including  means  for  applying  a high  frequency 
content  sampling  potential  difference  along  said  broad- 
band feedthrough  elements  relative  to  said  ground 
plane  for  accelerating  photoclectrons  emitted  by  said 
photoemissive  strip  toward  said  ground  plane  and  along 
said  electron  path,  and  means  coupled  to  said  housing 
for  receiving  and  detecting  electrons  passing  along  said 
electron  path. 

1 1.  The  head  of  claim  10  wherein  said  photoclcctron 
ejection  ground  plane  is  a metallic  mesh. 

12.  The  head  of  claim  10  wherein  said  photoemissive 
element  is  perforated  or  semitransparent. 

13.  The  head  of  claim  10  wherein  said  means  coupled 
to  said  housing  is  an  electron  multiplier. 

14.  T he  head  of  claim  10  wherein  a first  electron 
retarding  potential  is  applied  to  said  photoclcctron 
ejection  ground  plane  and  wherein  said  sampling  po- 
tential exceeds  said  first  retarding  bias  potential  and  is 
negative  w ith  respect  to  the  potential  of  said  photoelec- 
tron ejection  ground  plane. 

15.  The  sampling  head  of  claim  7,  wherein  said  pho- 
toelectron ejection  ground  plane  is  spaced  from  said 
strip  a distance  up  to  300  microns. 

16.  An  optoelectronic  sampling  head  comprising  a 
conductive  housing  defining  an  optical  path  in  the  axial 
direction  of  said  housing  for  impinging  light,  a pair  of 
diametrically  opposite  apertures  in  said  housing,  insula- 
tion means  scaled  in  said  apertures,  electrical  leads 
extending  through  said  insulating  means,  an  elongated 
photoemissive  strip  supported  on  said  leads  inside  said 
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housing,  a conductive  insert  positioned  within  said 
housing  and  having  affixed  thereto  a perforated  photo- 
electron  ejection  ground  plane  normal  to  the  axis  of 
said  housing,  said  ground  plane  being  spaced  from  said 
elongated  strip  a distance  up  to  300  microns  in  the  axial  5 
direction  of  said  housing,  and  electron  multiplier 
means  affixed  to  said  housing  in  the  optical  path  on  the 
end  thereof  toward  said  ground  plane,  whereby  said 
ground  plane  constitutes  a microstrip  ground  plane  for 
said  strip.  10 
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17.  The  sampling  head  of  claim  16  wherein  said  hous- 
ing is  provided  with  internal  threads,  said  insert  having 
external  threads  engaging  said  internal  threads 
whereby  the  axial  position  of  said  ground  plane  is  ad- 
justable, whereby  said  ground  plane  is  electrically  con- 
nected to  said  housing,  and  further  comprising  a source 
of  a sampling  potential  connected  to  said  leads,  and  a 
source  of  a bias  potential  connected  to  said  leads. 
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